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ABSTRACT 


Various  methods  are  presented  for  computing  heavy  or 
light  fluid  loading  (i.e.,  added  mass)  of  thin  finite 
rectangular  plates.  Based  on  the  results,  preferred  methods 
of  computation  are  recommended.  These  methods  and  a 
corresponding  computer  program --Opt ion  3- -are  of  particular 
value  in  extending  previously  formulated  digital  computer 
programs  for  obtaining  the  vibroacoustic  response  to  turbu¬ 
lence  excitation  of  a  plate.  Computer  results  are  given  for 
a  particular  case  which  includes  the  effect  of  fluid  load¬ 
ing  on  the  vibratory  response  of  a  plate  subject  to  turbu¬ 
lence  excitation. 


ADMINISTRATIVE  INFORMATION 

This  study  was  conducted  at  the  Naval  Ship  Research  and  Development 
Center  (NSRDC)  and  supported  by  the  Naval  Ship  Systems  Command  (NAVSHIPS) 
Code  037.  Funding  was  provided  by  NAVSHIPS  037  under  Subproject  S-4628008 
Task  14919. 

INTRODUCTION 

\ 

Reference  1  documents  four  available  computer  programs  for  determin-, 

ing  the  vibratory  response  and  the  associated  acoustic  radiation  of  a 

* 

finite  rectangular  plate  to  fully  developed  turbulence.  Several  compu¬ 
tational  frameworks  are  provided  which  can  be  modified  and  extended  through 
additional  research  to  furnish  more  accurate  and  realistic  programs  to  meet 
naval  needs.  The  chief  objective  of  the  original  study  was  to  furnish  a 
base  for  future  development. 

Extension  of  these  computations  are  treated  in  Reference  2  (Option  1) 
and  3  (Option  2).  Reference  2  includes  a  correction  in  the  computer  pro¬ 
grams  for  the  effects  on  vibroacoustic  response  of  the  boundary  layer 
thickness  and  pressure  pickup  dimensions.  Reference  3  includes  in  the  pro¬ 
grams  the  vibration  modes  and  natural  frequencies  of  thin  rectangular 
plates  with  clamped  and  rotational  supports  and  cylindrical  curvature. 

The  overall  program  now  includes  the  response  of  simple  and  clamped 
plates  in  air  and  in  water.  However,  the  fluid  loading  on  the  plate  has 

* 

References  are  listed  on  page  172. 
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hitherto  been  determined  for  the  case  of  an  infinite  plate.  Hence  there 
is  a  need  to  incorporate  into  the  program  the  results  of  recent  investi¬ 
gations  of  the  loading  effects  of  a  heavy  or  light  fluid  on  a  thin  finite 
rectangular  plate.  This  modification  is  required  to  improve  the  accuracy 
of  the  computations  and  to  extend  the  applicability  of  the  program. 

Accordingly,  the  present  report  presents  a  modification  (Option  3) 

of  coi  lodaied  version  of  oris  of  the  ovigivjxi  urograms — that  of  MaestrslLo — 

** 

to  include  the  effects  of  fluid  loading  for  finite  plates.  The  modified 
program  is  based  on  results  obtained  by  the  use  of  various  analytical 
methods.  The  following  titles  identify  the  methods  treated  and  their 
location  in  the  report;  notations  relevant  to  each  method  are  also  included 
in  the  Appendixes. 

Appendix  A  -  The  Feit-Junger  Method 

Appendix  B  -  The  Davies  Method 

Appendix  C  -  The  Leibowitz  Method  I 

Appendix  D  -  The  Leibowitz  Method  II 

Appendix  E  -  The  Bolt  Bei'anek  and  Newman  Method 

Appendix  F  -  The  Greenspon  Method 

Appendix  G  -  The  Leibowitz  Method  III 

For  the  convenience  of  the  reader,  the  Appendixes  include  an  adequate 
amount  of  the  mathematical  development  underlying  these  mathods.  An  under¬ 
standing  of  the  development  will  assist  the  reader  to  appreciate  the  merits, 
shortcomings,  subtleties,  and  complexities  of  a  particular  method  and  to 
apply  the  various  methods.  Certain  figures  are  adapted  from  the  basic 
references . 


* 

See  the  Dyer  representation,  includable  in  all  programs,  for  added  mass 
in  Appendix  A  of  Reference  1. 

** 

In  Appendix  H,  the  original  Maestrello  program  is  updated.  This  up¬ 
dated  version  is  then  modified  to  include  fluid  loading  (Option  3).  Other 
programs  presented  in  Reference  1  may  be  modified  in  a  similar  manner. 
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extend  the  applicability  of  the  results  for  Appendix  A  to  even-odd  and  odd-odd  nodes. 
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1.  Results  are  also  approxisately  correct  for  cylindrical  shells  of  interest  in  References  11  and  12. 

2.  Results  are  used  in  Reference  12  as  an  approxiaation  for  nodes  having  fairly  large  node  njsbers 
(khaik^a  large  where  a  is  the  shell  radius.)  As  kfca  becomes  large  the  nodal  vibration  acre  nearly 

approxicaces  beading  vibrations  on  a  flat  plate  for  which  tangential  cct.^ribut ions  to  kinetic 
energy  are  unisportant  compared  with  the  radial. 

Results  are  inapplicable 
vanish  for  these  nodes. 


■TED  AND 


m  ;„b  ,  fii  * 

",  '  -2-  '  V  _Si 

2*  R 

_  P  ‘i 


and  B^  depend  on  the  bean  functions  used  to  represent  the  node  shapes  and  therefore  depend 
cn  the  boundary  conditions  of  the  plate.  The  values  of  A^  and  for  several  of  the  lover  nodes 
of  plates  which  are  clasped  or  sieply  supported  on  all  edges  are  given  in  Table  A. 


,  .  (*  ) 
i  c  U  )  an  vacuun 

en  an  water  o 
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/  s  b  ..a  A* 
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V  2“  Bu 


gfT.O  AND 

t 


M  .  0.48s  (,b) *  A?  . 

*  — « —  _a  f 

=p  *u 

,  .  <,  > 
x  U  )  -  an  vacuua 

on  am  water  on 


1.  For  low  frequencies  the  present  procedure  does  not  require  that  f(~)  be  expllcitely 
deterained  fron  a  curve. 


<,"<■  .Me  i  o.<8i,,(lb)i  A^ 

V  1  °R  B1 I 


TABLE 


?  2 

Input  Data  for  Computing  Normalized  Modal  Mean  Square  Displacement  Y~(t^)/p 

Simply  Supported  Aluminum  Plate  in  Water 


for 


TABLE  2a 


Computed  Natural  Frequencies  and  Total  Damping  Ratio  for 
Simply  Supported  Plate  in  Air 


node 

San&er 

(R.ft) 

f(») 

KarturtOfl 

,  («*} 
cn'scc 

K>rtx:rtcr. 

i  *  -|2. 
f=n 

M 

83.82 

553.10 

0.16S33 

1.2 

345.83 

2179.54 

0.04324 

1.3 

776.93 

4831.94 

0.01930 

l.< 

1379.13 

8665.30 

0.01037 

1.5 

2153.31 

13529.62 

0.006% 

1,6 

3099.53 

19474.89 

0.004% 

2.1 

610.96 

0.15425 

2.2 

2232.40 

0.04221 

2.3 

735.40 

4934.80 

0.01909 

2,4 

1337.54 

8718.16 

0.01081 

2.5 

2161.72 

13532.48 

0.00695 

3.1 

111.26 

699.06 

0.13431 

3,2 

369.32 

2320.50 

0.04C61 

3,3 

799.42 

5022.90 

0.01876 

3.4 

1401.56 

8806.26 

0.01070 

3,5 

2175.74 

13670.58 

0.00690 

4,1 

130.89 

822.40 

0.11460 

4,2 

2443.84 

0.03356 

4,3 

819.05 

5146.24 

0.01831 

HI 

1421.19 

8929.60 

0.01055 

El 

2195.36 

13793.92 

0.00684 

5,1 

156.13 

980.93 

0.09607 

5.2 

414.19 

2602.42 

0.03621 

5,3 

844.29 

5304.22 

0.01776 

5,4 

1446.43 

9088.17 

0.01037 

5,5 

2220.61 

13952.49 

0.00681 

6,1 

186.93 

1174.80 

0.03022 

6,2 

445.03 

2796.23 

0.03370 

6,3 

875.13 

5498.63 

0.01714 

6,4 

1477.28 

9281-99 

0.01015 

6,5 

2251.46 

14146.31 

0.00666 

node 

UzB&er 

(«t.n) 

f(Ht) 

Ujrturton 

— ) 

TO  SW 

Uirburtoa 

4m  1  r~ 

ro 

7.1 

223.43 

1403.% 

0.05713 

7.2 

481.49 

3025.29 

0.03115 

7.3 

911.59 

5727.69 

0.01645 

7.4 

1513.73 

9511.(6 

0.00990 

7.5 

2237.91 

14375.37 

0.00655 

8.1 

265.49 

1663.15 

0.(6649 

3.2 

523.55 

3289.59 

0.02865 

3.3 

953.65 

5991.59 

0.01572 

8.4 

1555.79 

9775.35 

0.00964 

3,5 

2329.93 

14639.67 

0.00644 

9,1 

1957.63 

0.04789 

9,2 

Sh  9 

3539.13 

0.02625 

9.3 

1001.33 

6291.52 

0.01493 

9,4 

1603.47 

10074.83 

0.00935 

9.5 

2377.65 

14939.20 

0.00631 

10,1 

2302.46 

0.04093 

10,2 

3923.90 

0.02401 

10,3 

1054.61 

6626.30 

2.9K22 

10.4 

1656.75 

10409.66 

0.00905 

10,5 

2430.93 

15273.93 

0.00616 

a  =  3.0  ft 

b  =  0.541666  ft 

h  =  0.04  in. 

E  =  10  x  106  lb/in.2 

c  =  0.33 

weight  density  (o )  =  0.1075  lb/in.3  =  185.8  lb/ft3 

of  plate  "  alininua 

g  *  32.2  ft/sec2  *  384.6  in./sec2 
<  =  9.64  *  10'4  ft 

c.  =  17,000  ft/sec 


TABLE  2b 


Y2(<o) 


Input  Data  for  Computing  . . . L  for  Simply  Supported  Aluninua  Plate  in  Kater- 

2 

P 

Frequency  Computation  by  Method  I;  q  =  0 


neat 

immttr 

Kptl!Vft?) 

w^OVft7) 

{«•!«■) 

P  ^ 
(IWft7) 

CirfoMftV 

n 

KOI 

7.345* 

7.965? 

63.4* 

215.97 

1356.33 

3.55 

0.00573 

m 

5.52*3 

6.12*3 

37.56 

333.34 

2093.49 

4.61 

0.0014 

4.42*3 

5.0243 

25.24 

572.21 

3631.21 

5.63 

0.0031 

H 

iH 

1 

3.(974 

4.2876 

18.33 

918.20 

5766.33 

6.59 

0.0C21 

2.3 

■ 

I 

7.3257 

7.5257 

62.4? 

218.87 

1374.4? 

3.57 

0.0055 

2.* 

5.5115 

6.1115 

37.35 

4J7.53 

2747.73 

4.(3 

O.OOJ3 

2.i 

I 

<  4155 

S.0156 

25.15 

752.16 

472V60 

5.64 

O.CC23 

3.3 

9 

7.261/ 

7.2612 

61  33 

213.23 

Wt.t2 

3.60 

o.oca?? 

3.* 

V 

5.4*33 

6.0839 

37.01 

i  411.15 

2770.43 

4.65 

0.00335 

■ 

I 

4.4014 

5.0014 

25.01 

756.17 

4743.73 

5.65 

0.00237 

pH 

iH 

m 

7.1737 

7.7737 

60.43 

224.40 

1409.26 

3.64 

0.00516 

H® 

mt 

■ 

5.4*59 

6.9*59 

36.55 

4*6.01 

2200.96 

4.68 

0.00234 

pH 

h® 

9 

4.3817 

4.9817 

24.8? 

762.58 

4729.01 

5.67 

O.0C236 

■ 

K 

9 

10.0883 

10.(23^ 

114.23 

59.17 

*22.24 

2.65 

0.00850 

7.0(56 

7.665$ 

5?. 76 

231.52 

1453.98 

3.69 

0.00507 

CT® 

■ 

■ 

5.39=7 

5.9982 

JS.93 

451.05 

2251.4* 

4.71 

0.00330 

SpH 

■ 

4.3567 

4.9567 

24.57 

771.3? 

434*. 3* 

5.70 

0.00235 

6.2 

■ 

■ 

9.7321 

10.3321 

106.75 

1C8.38 

620  64 

2.74 

0.00X5 

6.3 

■ 

6.9432 

7.5130 

56.85 

2*0.63 

1511.20 

3.75 

0.CC496 

6.4 

■ 

5.3415 

5.9*15 

31-23 

*64.30 

2915.85 

4.76 

0.00326 

6.5  j 

4.3263 

4. 5263 

24.27 

722.62 

4914.86 

5.74 

0.CC233 

1.2 

9.2565 

9.C565 

SJ.K 

119.4* 

750.4* 

2.93 

0.00731 

1.3 

6.7999 

7.3993 

51.76 

2St.11 

1581. IS 

3.8? 

0.00*33 

7.4 

t 

5.2762 

5.870 

*75.83  j 

2934.3* 

4.81 

0.00221 

3.5 

<.7977 

4.8922 

93-93 

796.28  ; 

5oco.se 

S.J3 

C.CG?31 

3.2 

8.9722 

9.5773 

91.61 

132.45 

831.8* 

r-ss 

0.00709 

8.3 

6.6422 

7.2452 

M.« 

264.ce 

1664.12 

3.90  i 

O.COK3 

0.4 

s.resa 

5.8050 

33.30 

451.57 

3C87.J1 

4.87 

0.00315 

8.5 

4.2533 

4.8533 

W.55 

812.40 

5101.90 

5.2 2  ; 

0.00272 

9.? 

8.5903 

9.1933 

61.46 

147.*? 

926.25 

309 

0.00665 

9.3 

6.4331 

7.0031 

50.24 

220.32 

1760.44 

3.59 

0.00*53 

9.4 

5.1271 

5.7271 

32.20 

503.65 

3194.37 

4.91 

0.00309 

9.5 

4.210* 

4.8101 

23.14 

831.01 

5212.75 

5.22 

0.00225 

10.2 

e.2152 

3.2158 

77.72 

161.63 

1033.90 

3.21 

0.00520 

10.3 

6.3221 

6.9721 

47.92 

297.81 

1870.43 

4.03 

0.00*36 

10.< 

5.0*40 

5.64*0 

31.85 

522.08 

3316.3* 

5.01 

0.0030? 

1(1.5 

( 

>.6 

4.1640 

4.76*0 

22.70 

257.12 

5351.3? 

5.93 

0.00221 

3.0  ft 
0.5*1666  ft 
0.0«  in. 

_ _ 

17.000  ft/S« 

**  *  1.5  ft 
y’  *  0.270033  ft 


Co.  1  '  0.1075  lb/1*. 

jlistirxn 

9  *  38*. fc  in./s«? 

2  7 

V  -  1.0  (Ib/ftM 
4*  *  0.033  ft 


18S.8  lb/ft 


Mter 

'I  »  2.0 


0.223  lb/ft 

•  61.2  lb/ft3  ■  0.027  lb/i-..3 


i 

** 

r. 

i 

1 

1.6 

0.47 

2 

7.7 

3.0 

3 

12.0 

14.0 

Uc(ft/5«) 

‘(Wc)  |  | 

8 

2.25 

*  !0'3 

16 

3? 

_ j 

6- 

2.25 

»  10*  * 

2  ^ 


GwUttw 

cot 

p«rfc**<£ 
fcr  tl«e 
•Off 


Note:  In  computing  f  or  (see  Program  C'  in  Section  3  of  Appendix  H)  the  program 

input  in  in  lb/in. J  units  for  (o  )  In  computing  f  or  u  (see  Program  B'  in  Section 

w  aluainum  nn  1:115 


5  of  Appendix  H)  the  program  input  is  in  lb/ft-5  units  for  both  (p,. J a ^ umi jiinn  and  (P,.)w 


water 
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TABLE  2c 

onputation  of  Normalized  Modal  Mean  Square  Displacenent 
~  T 

Y~(w)/p“  for  Simply  Supported  Aluminum  Plate  in  h'ater 

AliWI.MN  PLATE 

■  •  o  (oxir) 


°c 

(fUsrz) 


B 

—  ,  -  •* 

p7 

T7?-) 

3.1  *  10'* 

3.1  x  10'S 

3.3  x  IO'10 

4.4  x  10'9 

4.5  x  10'*° 

4.8  i  io:" 

5.0  •  10'!0 

5.0  x  IO'*1 

5.5  .  IO'17 

6.5  *  IO'" 

6.6  x  IO'17 

7.3  .  IO'15 

1.9  *  10'* 

1.9  x  IO*9 

7.0  .  IO'10 

8.9  x  ID'10 

9.0  x  IO'11 

9.8  x  IO'17 

1.0  •  10'10 

1.0  x  10*" 

1.1  .  IO'17 

1.6  «  10'* 

1.6  x  IO'9 

1.7  .  IO'10 

1.2  .  ID'10 

7.3  .  IO11 

7.9  x  IO'17 

7.5  x  10'" 

7.6  .  IQ'17 

8.5  x  IO'13 

1.3  »  IO*8 

1.3  x  IO'9 

1.3  ,  IO'10 

6.*  x  10'10 

6.4  x  10'” 

7.0  x  IO'17 

6.7  <  10'" 

6.3  x  IO''7 

7.0  x  IO'13 

7.5  x  IO'7 

7.5  x  10'8 

7.6  .  IO'9 

9.9  x  10'9 

1.0  x  10'9 

1.1  «  IO'10 

5.7  ,  IO'10 

5.7  .  .0'” 

8.7  x  IO'17 

5.7  >  IO'11 

5.3  x  IO*17 

6.0  x  IO'13 

3.4  *  10'* 

3.4  x  10'9 

3.6  x  ID'10 

5.7  *  ID'9 

5.7  x  IO'10 

6.0  x  10'" 

7.1  x  10*'° 

7.1  x  10‘" 

7.6  x  IO'17 

8.7  x  10*" 

8.8  •  IO*17 

9.6  x  IO'13 

7.1  x  10'* 

7.1  .  IO'9 

7.7  x  IO'10 

1.7  x  IO'9 

1.7  .  IO''0 

1.3  .  10'" 

1.3  ,  IO''0 

1.3.10'" 

1.4  x  IO'17 

1.6  x  10'S 

1.6  x  IO'9 

1.7  x  IO'10 

8.7  x  10'*° 

8.7  x  10'" 

9.4  x  10',? 

9.1  x  10'" 

9.7  x  IO'17 

1.0  X  IO'17 

I.«  x  10'* 

1.4  x  IO'9 

1.5  x  IO'10 

8.«  x  IO'10 

8.5  x  10"” 

8.9  x  IO*17 

7.0  x  10'” 

7.0  x  IO'17 

7.8  x  IO"'3 

7.5  x  IO'7 

7.5  x  10'8 

7.6  x  10'9 

a 

1.0  x  10  0 

1.0  x  ID"9 

1.1  x  IO'10 

5.9  x  IO'10 

5.9  x  10'" 

6.4  x  IO'*7 

6.5  x  IO'17 

6.5  x  IO'17 

7.7  x  IO'13 

3.8  x  10'8 

3.8  x  IO'9 

4.0  x  IO'10 

7.1  x  IO'9 

7.1  ,  IO'10 

7.3  x  10'" 

note:  The  extensive  results  for  the  response 
:btained  by  Hr.  lucio  Heestrello  and  firs. 
Christine  Brown  at  Langley  Research  Center, 
HASA,  using  the  author’s  program,  are  tabu¬ 
lated  in  Table  Zc. 


) 


(LB/FT7) 


COMPUTED  DAMPING  VALUES  GIVEN  IN  TABLE  2b 
ONE-TENTH  COMPUTED  DAMPING  VALUES  GIVEN  IN  TABLE  2b 
TEN  TIKES  COMPUTED  DAMPING  VALUES  GIVEN  IN  TABLE  2b 


Figure  1  -  Normalized  Modal  Mean  Square  Displacement  of 
Simply  Supported  Aluminum  Plate  with  Fluid  Loading 
Effects  Included 


The  response  results  computed  at  NSRDC  and  plotted  in 
Figure  1  were  duplicated  by  Mr.  Lucio  Maestrello  and  Mrs. 
Christene  Brown  who  assisted  the  author  by  performing 
similar  computations  on  their  computer  at  Langley  Research 
Center  NASA,  using  the  author's  program.  Their  more  ex¬ 
tensive  results  are  tabulated  in  Table  2c. 
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TABLE  5 


Input  Data  for  Computing  Normalized  Modal  Mean  Square  Displacement  Y~(u)/p 
for  Simply  Supported  Steel  Plate  in  fc’ater 

TABLE  3a 

Computed  Natural  Frequencies  and  Total  Damping  for 
Simply  Supported  Steel  Plate  in  Air 


|  Mode 
Nunber 
(m,n) 


4.1 

5.1 

6.1 

7.1 

8.1 

9,1 

10,1 

11,1 

12,1 

13.1 

14.1 

15.1 

16.1 

17.1 

18.1 

19.1 

20.1 


f(Hz) 

Warburton 

umn(rad/sec) 

Warburton 

6  =  0.5/f 

mn  mn 

3.37 

21.16 

0.1483 

13.37 

83.98 

0.0373 

30.03 

188.68 

0.0166 

53.36 

335.25 

0.0093 

83.35 

523.71 

0.0059 

120.01 

754.05 

0.0041 

163.33 

1026.26 

0.0030 

213.32 

1340.35 

0.0023 

269.98 

1596.33 

0.0018 

333.30 

2094.18 

0.0015 

403.28 

2533.91 

0.0012 

479.93 

3015.52 

0.0010 

563,25 

2539.01 

0.0008 

653.23 

4104.38 

0.0007 

749.88 

4711.63 

0.0006 

853.19 

5360.75 

0.00058 

963.17 

6051.76 

0.00051 

1079.81 

6734.65 

0.00046 

1203.12 

7559.41 

0.00041 

1333.09 

8376.05 

0.00037 

a  =o  w 
mn  mn  mn 


a  =  10  ft 

b  =  97.39  ft 

h  =  0.5  in. 

E  =  30  x  106  lb/in.2 

a  =  0.30 

(p)  =  0.283  lb/in.3 
w  steel 

g  =  32.2  ft/sec2  =  384.6  in. /sec2 
k  -  1.204  x  10'2  ft 
c0  =  17,000  ft/sec 


TABLE  3b 


* 


Input  Data  for  Computing 


Y2(w) 


for  Simply  Supported  Steel  Plate  in  Water- - 


2 

P 

Frequency  Computation  by  Method  I;  q  =  0 


Mode 

Ntcber 

(o.n) 

Mp(lb/ft2) 

wCT(ib/ft2) 

p  cn 
(lb/ft2) 

2 

(M‘ }2(lb/ft2) 

?cn(h2) 

w  (rad/sec) 
cn 

S^O/sec) 

1.1 

2C 

.4 

406.57 

426.97 

182303.38 

0.74 

4.62 

0.0748 

0.03238 

2.1 

204.09 

224.49 

50355.76 

4.03 

25.33 

0.1425 

0.01125 

3,1 

136.16 

156.56 

24511.03 

10.85 

68.13 

0.2045 

0.00600 

4.1 

102.14 

122.54 

15016.05 

21.79 

136.83 

0.2612 

0.00381 

5.1 

81.72 

102.12 

10428.49 

37.28 

234.14 

0.3135 

0.00267 

6.1 

68.11 

38.51 

7834.02 

57.66 

352.12 

0.3617 

0.00199 

7.1 

58.38 

78.78 

6206.28 

83.18 

522.39 

0.4064 

0.00155 

8.1 

51.08 

71.48 

5109.39 

114.05 

716.23 

0.4479 

0.00125 

9.1 

45.41 

65.81 

4330.95 

150.43 

944.72 

0.4865 

0.00102 

10.1 

40.87 

61.27 

3754.01 

192.47 

1208.72 

0.5226 

0.00086 

11.1 

37.15 

57.55 

3312.00 

240.28 

1508.97 

0.5564 

0.00073 

12.1 

34.06 

54.46 

2975.70 

293.97 

1846.10 

0.5879 

0.00063 

13.1 

31.44 

51.84 

2687.38 

353.60 

2220.64 

0.6177 

0.00055 

14.1 

29-19 

49.59 

2459.16 

419.27 

2633.03 

0.6457 

0. 00049 

15.1 

27.25 

47.65 

2270.52 

491.03 

3083.68 

0.6721 

0.00043 

16.1 

25.54 

45.94 

2110.48 

568.94 

3572.93 

0.6966 

0.00038 

17,1 

24.04 

44.54 

1983.81 

653.04 

4101.07 

0.7190 

0.00035 

18.1 

22.71 

43.21 

1867.10 

743.37 

4668.38 

0.7411 

0.00031 

19,1 

21.51 

41.91 

1756.44 

839.98 

5275.09 

0.7641 

0.00028 

20.1 

20.4 

20.43 

40.93 

1675.26 

942.90 

5921.40 

0.7824 

0.00026 

—  2 
p2  =  1-0  (lb/ft2) 

rw  =  0.323  lb/ft2 

(p  )  =  64.2  lb/ft3  =  0.037  lb/in.3 
water 

a  =  2.0 


Kote:  5nn  ‘  3an  «V'K")  =  l°-5*}  (^0^V) 

m' 


Uc(ft/sec) 

e(sec) 

6*(ft) 

8 

B29ESI 

0.044 

16 

0.041 

32 

0.038 

64 

BEEHT8I 

0.035 

a  =  10  ft 
b  «  97.39  ft 
h  =  0.5  in. 

<  =  A  =  1.204  x  10‘2  ft 

ZV 3 

cf  =  17,000  ft/sec 

x  =  x‘  =  5.0  ft 
y  =  y*  =  48.695  ft 


V (lb/ft2)  (see  Colton  E) 
o  =  0.30 

E  -  30  x  106  lb/in.2 

(ow)  =  0.283  lb/in.3 
steel 

g  =  384.6  in./sec2 


1 

A. 

K. 

1 

1 

1 

1.6 

0.47 

2 

7.2 

3.0 

3 

12.0 

14.0 

results 

for 

inodes 

of 

interest 


Note:  In  computing  f  or  (see  Program  C'  in  Section  5  of  Appendix  H)  the  program 

input  is  in  lb/in.0  units  for  (P..)steej*  In  computing  f  or  w  (see  Program  B'  in 

Section  3  of  Appendix  H)  the  program  input  is  in  lb/ft°  units  for  both  (Pw)steej  and 

(o  )  „  . 

w'water 


'"t 


10 


q  =  o.3 


TABLE  3c  (Continued) 


p  mn 


(lb/ft2) 


?  5  2 

(H‘)^(1b/ft2) 


WHz>  “mn(rad/sec> 


*»(,/see)  *« 


199451.56 

54915.24 
26611.40 
16251.15 
11250.84 

8427.24 
6658.56 


5468.60 

4624.00 

3999.30 

3522.42 

3147.21 

2847.29 

2601.00 

2397.08 

2225.95 

2079.36 

1953.64 

1843.84 

1748.91 


0.00085  results 

0.00072  f°r 
modes 

0.00062  of  interest 


q  =  0.6 


TABLE  3d  (Continued) 


Mode 

Number 

(m,n) 

HpClb/ft2) 

wmn(W) 

M'=W+W 
p  mn 

(lb/ft2) 

?  p  2 

(W)  (lb/ft2) 

imn(rad/sec) 

amn(,/sec> 

6 

mn 

1  >1 

2,1 

3.1 

4.1 

5.1 

6.1 

7,1 

20.4 

508.21 

255.11 

528.61 

275.51 

279428.53 

0.66 

4.16 

0.0604 

0.02903 

75905.76 

3.64 

22.86 

0.1161 

0.01015 

170.20 

190.60 

36328.36 

9.83 

61.75 

0.1679 

0.00543 

127.68 

148.08 

21927.69 

19.82 

124.48 

0.2161 

0.00347 

102.16 

122.56 

15020.95 

34.04 

213.74 

0.2612 

0.00244 

85.14 

105.54 

11138.69 

52.80 

331.63 

0.3033 

0.00182 

72.98 

93.38 

8719.82 

76.41 

479.83 

0.3428 

0.00142 

8,1 

9,1 

10,1 

11,1 

12,1 

13.1 

14.1 

15.1 

16.1 

17  1 

63.86 

84.26 

7099.75 

105.05 

659.73 

0.3800 

0.00115 

56.76 

77.16 

5953.67 

138.93 

872.47 

0.4149 

0.00095 

51.09 

71.49 

5110.82 

178.19 

1119.03 

0.4479 

0.00080 

46.44 

66.84 

4467.59 

222.97 

1400.24 

0.4791 

0.00068 

42.57 

62.97 

3965.22 

273.38 

1716.80 

0.5085 

0.00059 

39.30 

59.70 

3564.09 

329.51 

2069.34 

0.5364 

0.00051 

36.49 

56.89 

3236.47 

391.46 

2458.40 

0.5628 

0.00045 

34.06 

54.46 

2965.89 

459.31 

2884.45 

0.6104 

0.00042 

31.93 

52.33 

2738.43 

533.11 

3347.90 

0.6119 

0.00036 

18,1 

19.1 

20.1 

30.05 

50.45 

2545.20 

612.92 

3849.13 

0.6347 

0.00032 

28.38 

49.28 

2428.52 

698.80 

4388.49 

0.6498 

0.00029 

26.89 

25.54 

47.29 

45.94 

2236.34 

2110.48 

790.80 

888.95 

4966.20 

5582.60 

0.6771 

0.6970 

0.00027 

0.00024 

results 

for 

modes 

of  interest 


TABLE  3e  (Continued) 


rt2) 

p  mn 

(lb/ft2) 

in 

(lb/ft2)* 

WH2> 

-  <r*J> 
unn's*c 

Smn<1/sec) 

6 

m,n 

932.7 

4 

953.14 

908475.86 

0.49 

3.10 

0.0335 

0.02161 

468.2 

0 

488.60 

238729.96 

2.734 

17.17 

0.0654 

0.00761 

312.3 

6 

332.76 

111394.74 

7.44 

46.93 

0.0962 

0.00415 

234.3 

3 

254.73 

64887.37 

15.11 

94.91 

0.1256 

0.00264 

187.4 

9 

207.89 

43218.25 

26.13 

164.11 

0.1540 

0.00187 

156.2 

r 

176.65 

31205.22 

40.82 

256.33 

0.1811 

0.00141 

133.9 

3 

154.33 

23817.75 

59.43 

373.24 

0.2073 

0.00111 

117.2 

0 

137.60 

18933.76 

82.21 

516.27 

0.2326 

0.00090 

104.1 

8 

124.58 

15520.18 

109.34 

686.68 

0.2570 

0.00074 

93.7 

6 

115.16 

13261.83 

141.01 

885.56 

0.2780 

0.00062 

85.2 

4 

105.64 

11159.81 

177.37 

1113.88 

0.3031 

0.00054 

78.1 

3 

98.53 

9708.16 

218.56 

1372.53 

0.3249 

0.00047 

72.1 

2 

92.52 

8559.95 

264.70 

1662.28 

0.3460 

0.00041 

66.9 

7 

87.37 

7633.52 

315.90 

1983.85 

0.3664 

0.00036 

62.5 

1 

82.91 

6874.07 

372.27 

2337.86 

0.3862 

0.00033 

58.6 

0 

79.00 

6241.00 

433.90 

2724.91 

0.4053 

0.00029 

55.1 

5 

75.55 

5707.80 

500.88 

3145.53 

0.4239 

0.00026 

52.0 

9 

72.49 

5254.80 

573.28 

3600.20 

0.4417 

0.00024 

49.3 

5 

69.85 

4879.02 

651.18 

4089.38 

0.4584 

0.00022 

46.8 

8 

67.28 

4526.90 

734.63 

4613.47 

0.4760 

0.00020 

fo 
mod 
of 
interest 


TABLE  3f  (Continued) 


Mode 

Number 

1,1 

0.237 

2,1 

1.325 

3,1 

3.638 

4,1 

7.445 

5,1 

12.972 

6,1 

20.410 

7,1 

29.331 

8,1 

41 .691 

9,1 

55.831 

10,1 

72.482 

11,1 

91.765 

12,1 

113.795 

13,1 

138.678 

14,1 

166.517 

15,1 

197.405 

16,1 

231.433 

17,1 

268.69 

18,1 

309.256 

19,1 

353.209 

20,1 

400.626 

results 

for 

modes 

of 

interest 


2 

Note:'  No  computations  of  -  were  made  for 


TABLE  3g 

Computation  of  Normalized  Modal  Mean  Square 

~2  ~2 

Displacement  Y  (w)/p  for  Simply  Supported 
Steel  Plate  in  Water 

Hote:  The  extensive  results  for  the  response  obtained  by  Hr.  Lucio  Kaestrello 
and  Hrs.  Christine  Brown  at  Langley  Research  Center,  NASA,  using  the  authors 
program,  are  tabulated  in  Tables  3g-3i. 
n  _  n  STEEL  PLATE 


Uc 

(ft/sec) 

Mode 

Humber 

(m.n) 

* 

** 

*** 

v2(u) 

P2 

Y2{u) 

P2 

Y2(u) 

P2 

8 

8,1 

4.4  x  10~9 

4.4  x  10"10 

4.4  x  10'11 

9,1 

1.7  x  10‘9 

1.7  x  10'10 

1.7  x  10'11 

10,1 

7.1  x  10'10 

7.1  x  10"11 

7.1  x  10-12 

11,1 

3.0  x  10'10 

3.0  x  10'11 

3.0  x  10'12 

12,1 

1.4  x  IQ*10 

1.4  x  10"11 

1.4  x  ID'12 

13,1 

6.4  x  10*11 

-12 

6.4  x  10 

6.4  x  ID’13 

14,1 

3.0  x  10'11 

-12 

3.0  x  10  xc 

3.1  x  10'13 

15,1 

1.6  x  10*11 

1.6  x  10'12 

1.7  x  10'13 

16,1 

9.3  x  10'12 

9.3  x  10‘13 

9.4  x  10'14 

17,1 

5.1  x  10'12 

3.1  x  19‘13 

5.2  x  10‘14 

18,1 

3.1  x  10'12 

2.1  x  10~13 

3.1  x  10'14 

19,1 

1.8  x  10'12 

1.8  x  ID'13 

1.9  x  10'14 

20,1 

1.0  x  10'12 

1.0  x  10'13 

1.1  x  10'14 

16 

8,1 

4.2  x  10'9 

4.2  x  10’10 

4.2  x  10'11 

9,1 

1.7  x  10'9 

1.7  x  10_1° 

1.7  x  10'11 

10,1 

7.5  x  ID'10 

7.5  x  10'11 

7.6  x  10'12 

11,1 

3.1  x  10'10 

3.1  x  10'11 

3.1  x  10'12 

12,1 

1.4  x  10'10 

1.4  x  10'11 

1.4  x  10"12 

13,1 

6.8  x  10'11 

6.8  x  10~12 

6.8  x  10'13 

14,1 

3.9  x  10'11 

3.9  x  10"12 

3.9  x  10'13 

15,1 

2.3  x  10'11 

2.3  x  10'12 

2.3  x  10'13 

16,1 

1.2  x  10'11 

1.2  y  10*12 

1.2  x  10'13 

17,1 

5.7  x  10~12 

5.7  x  ID'13 

5.7  x  10'14 

18,1 

3.1  x  10'12 

3.1  x  10'13 

3.2  x  10‘14 

19,1 

1.9  x  10'12 

1.9  x  10'13 

2.0  x  10‘14 

20,1 

1.2  x  10'12 

1.2  x  10'13 

1.2  x  10'14 

32 

8,1 

3.2  x  10'9 

3.2  x  10'10 

3.2  x  10*11 

9,1 

1.1  x  10‘9 

1.1  x  10'10 

1.1  x  10'11 

10,1 

4.6  x  10'10 

4.6  x  10'11 

4.7  x  10'12 

13 


TABU:  5!i 


STEEL  PLVTE 


q  =  0.3 

U  Mode 

<ft'sec>  "ST 


4.5  x  10* 
1.7  x  10' 

7.3  x  10' 

3.1  x  10' 

1.4  x  10' 

6.5  x  10' 

3.2  x  10' 

1.7  x  10" 

9.5  x  10* 
5.4  x  10' 

3.2  x  10" 
1.9  x  10' 

1.1  x  10' 

4.3  x  10" 

1.8  x  10' 

7.7  x  10' 

3.2  x  10' 

1.4  x  10' 

6.8  x  10' 

3.9  x  10' 

2.4  x  10" 

1.2  x  10" 

6.1  x  10' 

3.2  x  10' 
2.0  x  10" 

1.2  x  10' 

3.3  x  10" 
1.2  x  10' 
4.7  x  10' 


2.2  x  10" 


4.5  x  10' 
1.7  x  10" 

7.3  x  10' 

3.1  x  10' 

1 .4  x  10' 

6.5  x  10' 

3.2  x  10" 

1.7  x  10" 

9.5  x  10' 
5.4  x  10' 

3.2  x  10' 
1.9  x  10' 

1.1  x  10' 

4.3  x  10" 

1.8  x  10* 

7.7  x  10* 

3.2  x  10' 

1.4  x  10' 

6.8  x  10' 

3.9  x  10" 

2.4  x  10' 

1.2  x  10' 

6.1  x  10* 

3.2  x  10'1 
2.0  x  10" 

1.2  x  10" 

3.3  x  10' 
1.2  x  10' 1 
4.7  x  iO"1 


2.2  x  10' 


4.6  x  10 


1.8  x  10' 


7.3  x  10'“ 

3.1  x  10'12 

1.5  x  10'12 

6.6  x  !0'13 

3.2  x  10'13 

1.7  x  10'13 

9.6  x  10'14 
5.5  x  10'14 

3.2  x  10~14 
1.9  x  10'14 

1.1  x  10'14 

4.3  x  10'11 

1.8  x  10'11 

7.7  x  10'12 

3.2  x  10'12 

1.4  x  10'12 

6.9  x  10'13 
4.0  x  10'13 

2.4  x  10"13 

1.2  x  10'13 

6.1  x  10'14 

2.3  x  10"14 
2.0  x  10'14 

1.3  x  10'14 

3.3  x  10'11  c2 

1.2  x  10'11  c 

4.7  x  10'12  a,uminum  243,000 

_lz  stee1  19,400 

2.3  x  10  lc 


Table  3h  (Continued) 

STEEL  PLATE  (Continued) 


q  =  0.5 


Uc 

(ft/sec) 

Kode 

Nieier 

(=,n) 

_  ** 

H 

Y2^) 

P2 

Y2^) 

P2 

12.1 

1.3  X  10"10 

1.3  x  10*11 

1.3  x  10*12 

13.1 

8.8  x  ID’11 

8.8  x  10'12 

8.8  x  lO*13 

14,1 

6.6  x  10'11 

6.6  x  IQ'12 

6.6  x  10'13 

15.1 

4.7  x  10'11 

4.7  x  10'12 

4.7  x  10"13 

16.1 

2.5  x  10'11 

2.5  x  10'12 

2.5  x  10-13 

17,1 

1.0  x  10-11 

1.0  x  10'12 

1.0  x  10*13 

18,1 

4.4  x  10'12 

4.4  x  10'13 

4.5  x  10"K 

19,1 

1.8  x  10'12 

1.8  x  10'13 

1.9  x  10~14 

20,1 

6.0  x  10*13 

6.1  x  10‘14 

6.3  x  lO’15 

64 

8,1 

3.6  x  10‘9 

3.6  x  lO*10 

3.6  x  10'11 

9,1 

1.4  x  10'9 

1.4  x  10*10 

1.4  x  10'11 

10,1 

5.2  x  10‘10 

5.2  x  10-11 

5.3  x  10‘12 

11,1 

2.1  x  10'10 

2.1  x  10'11 

2.2  x  10‘12 

12,1 

7.7  x  10"11 

7.7  x  10-12 

7.8  x  10'13 

,3., 

2.2  x  10'11 

2.2  x  10'12 

2.3  x  10'13 

14,1 

t 

f 

15,1 

t 

t 

A 

16,1 

t 

r 

t 

17,1 

t 

■ 

4. 

18,1 

t 

t 

A. 

19,1 

t 

t 

f 

20,1 

t 

t 

t 

For  one-tenth  computed  damping  value  given  in  Table  3b. 

** 

For  computed  damping  value  given  in  Table  3b. 

*** 

For  ten  times  computed  damping  value  given  in  Table  3b. 
‘Computer  error  {recomputation  not  made). 


TABLE  3i 


°C 

Hote 

l^(-) 

T2^) 

r2^) 

(ft/sec) 

Kunbsr 

(o,n) 

P2 

P2 

p2 

8 

8.1 

4.9  x  10'9 

4.9  x  10"10 

4.9  =  10"11 

9.1 

1.9  x  10"9 

1.9  x  10‘10 

1.9  x  tO'11 

10.1 

7.9  x  10‘10 

7.9  x  10*11 

8.0  x  10’12 

11.1 

3.5  x  IQ’10 

3.5  x  10-" 

3.5  x  10'12 

12,1 

1.6  x  10-10 

1.6  x  10"11 

1.6  x  10'12 

13.1 

7.3  x  IS'11 

JX'J 

7.3  x  10 

7.4  x  ID*13 

14.1 

3.7  x  10-" 

3.8  x  10* iV 

3.8  x  ID*13 

15.1 

1.9  x  10-11 

1.9  x  10'12 

1.9  x  10'13 

16,1 

1.0  x  10'11 

1.0  x  10'12 

1.0  x  10'13 

17,1 

5.6  x  10*12 

5.6  x  10‘13 

5.7  x  10'U 

18,1 

3.2  x  10*12 

3.2  x  10"13 

3.2  x  10‘K 

19,1 

1.9  x  10'12 

1.9  x  10"13 

1.9  x  10'14 

20,1 

1.3  x  10'12 

1.3  x  ID*13 

1.3  x  10'14 

16 

8,1 

4.5  x  10*9 

4.5  x  10*10 

4.5  x  10*11 

9,1 

1.9  x  10'9 

1.9  x  10‘10 

1.9  x  10'11 

10,1 

8.4  x  10_1C 

8.4  x  10'10 

8.5  x  10‘12 

11,1 

3.7  x  10'10 

3.7  x  10'11 

3.8  x  10~12 

12,1 

1.6  x  10*10 

1.6  x  10"11 

1.6  x  10*12 

13,1 

7.3  x  10*11 

7.3  x  10‘12 

7.4  x  10'13 

14,1 

4.1  x  10*11 

4.1  x  10'12 

4.2  x  10"13 

15,1 

2.4  x  10'11 

2.4  x  10*12 

2.4  x  10'13 

16,1 

1.4  x  10'11 

1.4  x  10-12 

1.4  x  10'13 

17,1 

7.1  x  10'12 

7.1  x  10'13 

7.1  x  10'14 

18,1 

3.5  x  10'12 

3.5  x  10‘13 

3.5  x  10'14 

19,1 

2.0  x  1C'12 

2.0  x  i0'13 

2.0  x  10'14 

20,1 

1.3  x  10'12 

1.3  x  10*13 

1.4  x  10‘14 

32 

8,1 

3.6  x  10'9 

3.6  x  10'10 

3.6  x  10'11 

9,1 

1.3  x  10'9 

1.3  x  IQ*10 

1.3  x  10'11 

10,1 

5.0  x  10'10 

5.0  x  10'11 

5.1  x  10'12 

<3  i  0.6 


Table  3i  (Continued) 

STEEL  PLATE  (Continued) 


Node 

X-jiber 

(«=,n) 

,Jc 

(ft/sec) 

^(u) 

P2 

T2^) 

7 

r2^) 

p2 

11.1 

2.3  x  10'10 

2.3x10-" 

2.3  x  10"12 

12,1 

1.3  x  ID’10 

1.3  x  10'" 

1.3  x  10"12 

13,1 

8.1  x  10'11 

8.1  x  10"12 

8.2  x  10"13 

14.1 

6.0  x  ID'" 

6.0  x  10"12 

6.0  x  10"13 

15,1 

4.4  x  10‘" 

4.4  x  10"12 

4.4  x  10"13 

16,1 

3.0  x  10"" 

3.0  x  10"12 

3.0  x  10"13 

17.1 

1.5  x  10‘" 

1.5  x  10"12 

1.5  x  10‘13 

18,1 

6.1  x  10"12 

6.1  x  10"13 

5.2  x  10"K 

19,1 

2.7  x  10"12 

2.7  x  10"13 

2.8  x  10"14 

20,1 

1.1  x  10"12 

1.1  x  10"13 

1.2  x  10K 

64 

8,1 

3.9  x  10"9 

3.9  x  10"10 

4.0  x  10"" 

9,1 

1.5  x  10"S 

1.5  x  10*10 

1.5  x  10“" 

10,1 

6.0  x  10"10 

6.0  x  10"" 

6.0  x  10"12 

11,1 

2.5  x  10-10 

2.5  x  10"" 

2.5  x  10"12 

12,1 

9.9  x  10"" 

9.9  x  10"12 

1.0  x  10"12 

13,1 

3.6  r  10"11 

3.6  x  10"12 

3.7  x  10"13 

14.1 

4.8  x  10"12 

4.9  x  10*13 

5.3  x  10"K 

15,1 

r 

* 

16,1 

t 

f 

17,1 

f 

f 

18,1 

f 

t 

19,1 

f 

f 

20,1 

f 

f 

f 

* 

For  one-tenth  computed  damping  value  given  in  Table  3b. 

** 

For  computed  damping  value  given  in  Table  3b. 

***  . 

For  ten  tires  computed  damping  value  given  in  Table  3b. 

' Computer  error  (recomputation  not  made). 
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Uc 

(ft/sec 

Mode 

Kucher 

(n,n) 

** 

^(o) 

PT 

T2^) 

PZ 

Y2^) 

7 

11.1 

3.2  x  10-W 

3.2  x  ID"” 

3.2  x  10"12 

12,1 

1.4  x  10'10 

1.4  x  10"11 

1.4  x  10"12 

13,1 

7.4  x  10'11 

7.4  x  10"12 

7.5  x  10"13 

14,1 

4.6  x  1C'11 

4.6  x  ID"12 

4.6  x  10‘13 

15,1 

3.0  x  10'11 

3.0  x  10"12 

3.1  x  10"13 

16,1 

2.4  x  10"11 

2.4  x  10‘12 

2.4  x  10"13 

17,1 

2.1  x  10'11 

2.1  x  10"12 

2.1  x  10"13 

18,1 

1.6  x  10'11 

1.6  x  10"12 

1.6  t  10"13 

19.1 

9.4  x  10'12 

9.4  x  10"13 

9.4  x  10'14 

2G.1 

4.4  x  ID"12 

4.4  x  10"13 

4.4  x  10'14 

64 

Not  Computed 

9,1 

2.3  x  10'9 

2.3  x  10*10 

2.3  x  10'11 

10,1 

9.5  x  10'10 

9.5  x  10'11 

9.6  x  10"12 

11,1 

4.1  x  10"10 

4.1  x  10"11 

4.2  x  10"12 

12,1 

1.8  x  10‘10 

1.3  x  10"11 

1.8  x  10"12 

13,1 

8.5  x  10"11 

8.5  x  10"12 

8.6  x  10"13 

14,1 

3.9  x  lO"11 

3.9  x  10"12 

3.9  x  10'13 

15,1 

1.4  x  10'11 

1.4  x  10'12 

1.4  x  10'13 

16,1 

2.8  x  10"12 

2.9  x  10'13 

3.0  x  10'14 

17,1 

t 

f 

f 

18,1 

t 

f 

f 

19,1 

t 

f 

f 

20,1 

t 

T 

t 

* 

For  one-tenth  computed  damping  value  given  in  Table  3b. 

** 

For  computed  damping  value  given  in  Table  3b. 

*** 

For  ten  times  computed  damping  value  given  in  Table  3b. 

'"'Computation  not  made  for  8,1  mode  because  for  q  =  0.9, 
fg  j  <  100  Hz  (see  Table  3e). 


'Computer  error  (recomputation  not  made). 


<LB/F 


Figure  2  -  Normalized  Modal  Mean  Square  Displacement  of  Simply 
Supported  Steel  Plate  Kith  Fluid  Loading  Effects  Included 

The  response  results  computed  at  .VSSDC  and  plotted  in 
Figure  2  were  duplicated  by  Mr.  Lucio  Maestrello  and  Mrs. 

Christens  Brown  who  assisted  the  author  by  performing 
similar  computations  on  their  computer  at  Langley  Research 
Center  .VASA,  using  the  author's  prog  ran.  Their  Bore  ca¬ 


ll,  WFT/SEC 

Figure  2a 


Figure  2b 
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CALCULATIONS  AND  RESULTS 

The  analytical  results  for  a  fluid-loaded  plate  are  presented  in  a 

series  of  Appendixes  (A-G).  The  salient  features  of  this  study  are 

summarized  in  Table  1  for  the  convenience  of  the  reader. 

In  addition,  computer  results  were  obtained  for  the  vibratory 

response  of  a  water- loaded,  simply  supported,  rectangular-aluminum  isotropic 

plate  and  for  a  water- loaded,  simply  supported,  rectangular-steel  isotropic 

plate  subject  to  turbulence  excitation.  The  computer  results  were  obtained 

by  modifying  the  updated  version  of  the  original  formulation  devised  by 

Maestrello *  with  the  fluid  loading  program  devised  here.  The  corresponding 

mathematical  analysis,  methods  for  determining  the  imput  data,  and  computer 

* 

program  documentation  are  given  in  Appendix  H.  Results  are  now  presented 
for  the  aluminum  and  steel  plates. 

ALUMINUM  PLATE 

Table  2a  presents  the  natural  frequencies  and  total  damping  ratio  of 
a  simply  supported  aluminum  plate  in  air  for  several  modes  of  vibration. 

The  natural  frequencies  were  computed  by  use  of  the  Warburton 
method  (see  Reference  3  for  the  associated  computer  program)  and  were 

checked  by  use  of  the  simple  frequency  expression  u>mn  =  kc^  (T)2  *  (r)2] 
for  simply  supported  plates.  The  corresponding  natural  frequencies  in  water, 
as  well  as  other  relevant  input  data  required  for  the  computation  of  the 

Y2(o)) 

normalized  modal  mean  square  displacement — == — .  are  tabulated  in 

p2  -T~ 

**  Y  (o>) 

Table  2b.  The  normalized  modal  mean  square  displacements  - —  for  the 

2 

P 

simply  supported  aluminum  plate  in  water  are  tabulated  in  Table  2c.  All  of 
these  data  are  computed  in  accordance  with  the  detailed  procedure  given  in 
Appendix  H. 


Both  the  updated  program  (designated  MTURAD)  which  supercedes  the  ori¬ 
ginal  Maestrello  program  (designated  TURAD)  and  the  fluid  loading  program 
(Option  3)  are  presented  in  Appendix  H. 

**  12 
The  Maestrello  method  for  obtaining  p  is  given  in  Appendix  H.  Dyer 

2  -31  22 

(page  32  of  Reference  1)  uses  p  =  [6  x  10  •  r  p  UB  ]  and  Jacobs  (page  301 

2  2  ^ 

of  Reference  1)  used  p  =3.1  x  .  Here  p  is  the  mean  square  turbulence 
pressure,  p  is  the  fluid  density,  U^  is  the  free- stream  velocity,  and  Tw  is 
the  local  wall  shear  stress  (see  Appendix  H  for  further  discussion). 
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Figure  1  is  a  plot  of  some  of  the  computed  results  for  the  normalized 

~  ~2 

modal  mean  square  displacements  Y“(u)/p  of  a  simply  supported,  water-loaded 

aluminum  plate  subject  to  turbulence  excitation  over  a  range  of  convection 

velocities  U  for  various  values  of  damping  and  q  =  k/k  =  0  only  (see 
c  *  s 

notation  for  Appendix  H) .  More  extensive  results  are  tabulated  in  Table 
2c.  The  computer  procedure  used  for  the  calculations  is  given  in  Appendix  H. 

STEEL  PLATE 

Table  5a  tabulates  the  natural  frequencies  of  a  simply  supported 
steel  plate  in  air  for  several  modes  of  vibration  as  computed  by  use  of  the 
Ivarburton  method  (see  Reference  3  for  the  associated  computer  program) . 

The  corresponding  natural  frequencies  in  water  as  well  as  other  relevant 
input  data  required  for  the  computation  of  the  normalized  modal  mean  square 
~2  ~ 

displacement  Y~(co)/p“  are  tabulated  in  Tables  3b-3f.  The  normalized  modal 

2  2 

mean  square  displacements  Y  (w)/p  for  the  simply  supported  steel  plate  in 

water  are  tabulated  in  Tables  3g-3j . 

Figures  2a-2d  are  plots  of  some  of  the  computed  results  for  the  nor- 

2  2 

malized  modal  mean  square  displacement  Y  (<a)/p  of  a  simply  supported,  water- 

loaded  steel  plate  subject  to  turbulence  excitation  over  a  range  of  con- 

* 

vection  velocities  U£  for  various  values  of  damping  and  q  =  k/ks  <  1  (see 
notation  for  Appendix  H) .  More  extensive  results  are  tabulated  in  Tables  3g- 
3j.  The  computer  procedure  used  for  the  calculations  is  given  in  Appendix  H. 

DISCUSSION  AND  EVALUATION 

This  section  discusses  (1)  the  analytical  results,  (2)  the  compu¬ 
tational  results,  and  (3)  the  turbulence-vibroacoustic  relationships. 

it 

For  the  present  problem,  the  contribution  of  radiation  damping  but  not 
added  mass  is  excluded  in  the  (acoustically  slow)  region  for  which  the  com¬ 
putations  are  applicable;  see  Appendix  H  for  a  more  detailed  discussion  of 
this  point.  Simple  equations  for  determining  the  radiation  damping  con¬ 
tributions  for  inclusion  in  the  computer  program  will  appear  in  a  companion 
report . 


ANALYTICAL  RESULTS 


Table  1  identifies  and  compares  the  various  methods  of  computation. 
Notation  pertinent  to  each  method  is  found  in  Appendixes  A-G  inclusive. 

The  results  have  been  included  in  the  overall  computer  program  (see 
Appendix  H) . 

From  the  summary  for  Appendix  C  given  in  Table  1,  we  observe  that 
the  solution  to  the  Feit-Junger  basic  equations  has  been  extended  to  in¬ 
clude  the  added  mass  results  for  coupled  modes  in  addition  to  their  results 
for  uncoupled  modes  presented  in  Appendix  A;  thv,  extended  results  are  then 
identical  to  corresponding  results  obtained  by  Davies  in  Appendix  B.  The 
physicomathematical  basis  clarifying  the  plate-fluid  coupling  mechanism  is 
discussed  in  Appendix  A.  In  general,  results  can  be  classified  both  in 
terms  of  modes  in  wave  number  space  according  to  their  radiation 
characteristics  and  in  terms  of  frequency  band. 

It  is  interesting  to  note  that  despite  the  variety  of  analytical 
methods  used,  the  results  obtained  in  Appendixes  A-E  are  identical  for  the 
uncoupled  modes.  As  discussed  in  Appendixes  A-C,  the  results  for  these 
modes  are  considered  to  be  dominant  with  respect  to  the  results  for  coupled 
modes.  Moreover,  we  perceive  from  the  discussion  in  Appendix  D  that  for  the 
uncoupled  (dominant)  modes,  the  methods  of  Appendixes  A-D  are  applicable, 
with  minor  modification,  to  low  wave  numbers  (m,n=(l,l)  (1,2)  (2,1)...)  as 
well  as  to  high  wave  numbers. 

The  results  presented  in  Table  1  for  Appendixes  A-E  are  based  on 
analyses  which  assume  a  plate  with  simple  supports  for  the  boundary  con¬ 
dition.  However,  with  proper  modification  and  interpretation,  these 
results  can  be  used  to  yield  corresponding  results  for  fluid-loaded  clamped 
plates,  at  least  for  the  dominant  modes.  Following  the  method  and  using 
the  computer  program  (Option  2)  described  in  Reference  3,  we  first  obtain 

the  in  vacuo  natural  circular  frequency  a>  ■*  (w  )  ,  ,  for  a  clamped- 

*  n  mn  mn  clamped  r 

clamped  plate.  The  frequency  for  the  fluid-loaded  elamped-clamped  plate  is 


w  (w  )  i  ,  means  w  =  (w  )  ,  ,  and  similarly  for 

mn  mn  clamped  mn  mn  clamped  1 


^°mn^  damned' 
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then  computed  from  the  equation  for  to  ■*  (to  )  ,  .  given  in  Table  1 

r  mn  mn  clamped 

where  now,  in  this  equation,  to^  (tonn)cj ampe(j  the  previously  computed 

in-vacuo  frequency  for  a  elamped-clamped  plate.  Finally  to  obtain  the 
response  of  a  fluid-loaded  clamped-clamped  plate  subject  to  turbulence  ex¬ 
citation,  we  use  (“jnjjicianjpgjj  hut  mode  shapes  corresponding  to  simply 
supported  end  conditions.  In  extending  the  method  of  Reference  3  to  the 
case  of  fluid- loaded  clamped-clamped  plates }  we  assume  that  the  sensitivity 
of  the  fluid  loading  to  the  change  in  boundary  conditions  is  sufficiently 
small  so  that  the  preceding  procedure  will  yield  approximately  correct 
results  for  the  response  of  clamped-clamped  plates. 

Using  an  entirely  different  theoretical  approach.  Appendixes  F  and 
G  (see  Table  1)  give  directly  obtainable  results  for  both  simply  supported 
and  clamped  plates.  Again  the  results  depend  on  the  mode  numbers. 

Appendix  G  presents  a  relatively  simple  formulation  relevant  for  low  fre¬ 
quencies. 

If  we  treat  the  radiation  of  boundary  layer  noise  into  a  closed 

* 

rectangular  cavity  as  well  as  into  free  space  (see  Dyer  Model,  Appendix  A 
of  Reference  1),  then  the  approximate  value  of  the  added  mass  due  to  the 
enclosed  fluid  is  given  as 


p  tanhlk  L  I 
1  mn  z' 


cavity 


if  a)  >io 
mn  c 


n  =  4- —  tan  k  L  if  to  <  to 

cavity  k  mn  z  mn  c 

J  mn 


The  total  added  mass  is  then  the  sum  of  the  added  mass  of  the  enclosed 
fluid  and  the  free  (half)  space  added  mass  given  in  Table  1. 


In  the  Dyer  Model  all  interior  surfaces  except  the  plate  are  assumed  to 
be  pressure  release  surfaces. 
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COMPUTATIONAL  RESULTS 


For  the  results  shown  in  Figures  1  and  2,  which  include  the  effects 

of  fluid  loading,  we  make  the  following  observations;  of  course  evaluation 

of  the  theoretical  results  requires  a  comparison  between  theory  and  experi- 
* 

ment  in  water.  _  _ 

?  ? 

To  compute  the  normalized  modal  mean  sequare  displacement  Y“(a))/p“ 

of  a  water-loaded  plate  subject  to  turbulence  excitation,  we  require  the 

corresponding  modal  frequencies.  Appendix  H  presents  two  methods, 

designated  Methods  1  and  2  for  computing  the  natural  frequencies  of  a 

plate  in  water.  Method  1  (Equation  (H3))  is  based  on  the  analyses  in 

Appendixes  A-E  whereas  Method  2  (Equation  (H6))  is  based  on  the  analyses  in 

Appendixes  F  and  G.  Computations  show  that  the  frequencies  computed  by 

Method  2  are  significantly  greater  than  those  computed  by  Method  1  and  this 

** 

difference  increases  at  higher  modes.  The  discrepancy  is  attributed  to 

the  less  sophisticated  assumptions  involved  in  deriving  Equation  (H6) ;  see 

Appendix  F.  Consequently,  we  consider  the  results  obtained  by  Method  2  to 

be  less  valid  (i.e.,  more  inaccurate)  than  those  obtained  by  Method  1,  and 

~~2  2 

use  only  Method  1  for  the  computations  of  Y  (w)/p  now  presented. 

For  the  fluid-loaded  aluminum  plate.  Figure  1  shows  the  contributions 

2 

of  the  lowest  modes  to  the  normalized  modal  mean  square  displacement  Y  (w)/p 
for  both  the  range  of  convection  velocities  considered  (0  <  Uc  <  64  ft/sec) 
and  the  range  of  modal  frequencies  of  interest  (100  <  f  <  1000  Hz) ;  see 
Table  2.  It  is  evident  from  the  figure  that  for  any  damping  value  con¬ 
sidered,  the  contribution  of  the  modes  to  the  total  normalized  mean  square 
displacement  (Y^  +  Y^^  +  Y^  +  Y^4  )+..-)/p  decreases  with  successive 
mode  orders  (1,2),  (2,3),  (1,4),  and  (2,4).  The  major  contributors  are  the 
(1,3)  and  (2,3)  modes;  the  relative  contribution  of  the  latter  mode 


* 

The  fluid  loading  does  not  include  the  influence  of  hydrostatic 
pressure.  The  effect  of  hydrostatic  pressure  on  the  natural  frequencies  is 
discussed  in  Appendix  I. 

** 

As  shown  in  Appendixes  F,  G,  and  H,  the  data  given  in  Table  4  are  used 
in  computing  the  frequencies  by  Method  2. 
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TABLE  4 


Modal  Values  of  A. .  and  B- -  for  Clamped  and 
ij  ij  * 

Simply  Supported  Plates 


Mode 

Clamped  Plate 

Simply  Supported  Plate 

A.. 

ij 

B-  - 
ij 

A.. 

ij 

B-  - 
ij 

i  =  1  j  =  1 

0.6904 

1 

■a 

i  =  1  j  =  2 

0 

1 

mm 

i  =  1  j  =  3 

0.3023 

1 

0.25 

i  =  1  j  =  5 

0.1924 

1 

0.25 

i  =  3  j  =  1 

0.3023 

1 

0.1351 

0.25 

i  =  3  j  =  2 

0 

1 

0 

0.25 

i  =  3  j  =  3 

0.1324 

1 

0.0450 

0.25 
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decreases  for  values  of  approaching  64  ft/sec.  However,  we  also  ob¬ 
serve  that  for  given  values  of  Uc  and  damping,  the  root  mean  square  dis¬ 
placement  of  each  mode  treated  is  of  the  same  order  of  magnitude.  Ke  con¬ 
clude  therefore  that  each  of  the  several  modes  considered  makes  a  signi¬ 
ficant  contribution  to  the  total  vibrator}'  displacement. 

A  comparison  of  the  curves  of  Figure  1  also  indicates  that  the 
character  of  the  curves  is  essentially  independent  of  the  damping  values 
and  that  the  values  for  the  modal  amplitudes  decrease  with  increasing 

~ 7 

values  of  the  damping.  If  the  damping  is  sufficiently  small  and/or  p~  is 
sufficiently  large,  we  may  expect  the  turbulence  excitations  to  produce 
undesirably  large  magnitudes  of  plate  vibration. 

For  the  steel  plate  (Figure  2),  similar  results  were  obtained  for 

the  lowest  modes  for  various  values  of  q  =  k/kg;  see  Table  3.  Again,  we 

observe  that  for  any  value  of  q,  each  of  the  several  modes  considered  makes 

a  significant  contribution  to  the  total  vibratoiy  displacement,  the  modal 

amplitudes  decrease  with  increasing  values  of  damping,  and  the  character 

of  the  curves  is  essentially  independent  of  the  damping  values.  In 

addition,  we  observe  that  the  response  increases  with  increasing  values  of 

q.  For  q  =  0.995,  Table  3b  shows  a  significant  increase  in  the  mode 

numbers  corresponding  to  the  frequency  range  of  interest  (100  <  < 

1000  Hz) .  Once  more  we  conclude  that  if  the  damping  is  sufficiently  small 
2 

and/or  p  is  sufficiently  large,  we  may  expect  the  turbulence  excitations 
to  produce  undesirably  large  magnitudes  of  plate  vibration;  the  magnitudes 
are  enhanced,  but  not  radically,  at  sufficiently  large  values  of  q.  Hence 
in  computations  it  appears  practical  to  use  a  single  representative  value 
for  q  for  the  range  of  turbulence  frequencies  distributed  about  u>  =  Uc  Kj. 

=  k  c,  corresponding  to  0  <  q  =  k/ks  <1.0.  A  practical  alternative  would 
be  to  select  the  average  of  the  mean  square  displacement  responses  computed 
for  small  and  large  values  of  q. 


* 

These  relationships  are  discussed  in  the  following  section. 
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TURBULENCE- VI BROACOUSTIC  RELATIOXSHI PS 

Background  material  is  cited  in  the  bibliography. 

For  the  convenience  of  the  reader,  we  list  the  notation  consaonly 
used  in  this  section. 

a,b  Length  and  width  of  the  plate,  respectively 

3  "> 

B  Plate  bending  stiffness  equal  to  E  h  /12(l-o“) 

C„  Trace  speed  of  the  plate  bending  wave  in  the  direction  of  flow  or 

B  1/2  1/4 

free  flexural  phase  velocity  for  a  thin  plate  equal  to  ta  (B/M) 

c  Velocity  of  sound 

2  ,1/2 

Compress ional  wave  velocity  of  the  plate  equal  to  [E/ps(l-o  )] 

E  Young's  modulus 

f ,fc ,r^  Natural,  acoustic  coincidence  or  critical,  and  hydrodynamic  coinci¬ 
dence  or  critical  frequencies,  respectively 

h  Plate  thickness 


ki'k5 

k  ,k 
ra  n 

k  ,k 
mn  s 


Acoustic  wave  nuaber  equal  to  u/c 

Save  number  exponents  lying  along  the  x-  and  v-axes,  respectively 
Modal  wave  nunbers  equal  to  nm/a  and  nn/b,  respectively 

r~ 7  7 

Save  nunbers  equal  to  ykQ“  +  k^ 

1/2 

Free  plate  bending  wave  nuaber  equal  to  (w/sc  c^) 

Turbulence  wave  nuaber 

Plate  structural  aass  per  unit  area 

Effective  mass  per  unit  area  (i.e.,  mass  per  unit  area  of  fluid- 

loaded  plate)  equal  to  m  +  m 

r  n  p  rnn 

Added  mass  per  unit  area  (or  apparent  mass  or  virtual  mass  per  unit 
area) 

Mode  numbers 
Equal  to  k/kg 

Free-stream  velocity 

Convection  velocity  equal  to  0.8  U 

Mean  convection  speed  along  the  flow  direction 
Hydrodynamic  coincidence  speed 

Equal  to  1  for  fluid  loading  on  one  side  of  plate  only;  equal  to  2 
for  fluid  loading  on  both  sides  of  plate 

Radius  of  gyration  equal  to  h/2f/J 
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p  Mass  density  of  fluid  medium 

Ps  Mass  density  of  plate 

o  Poisson's  ratio 

u  Circular  frequency  equal  to  2nf 

ii>c  Acoustic  coincidence  or  critical  frequency  equal  to  2nfc 

Hydrodynamic  coincidence  or  critical  frequency  equal  to  2»fg 

w  Plate  resonance  frequency  equal  to  2irf 

ran  n  n  mn 

-  Bar  over  quantity  denotes  quantity  for  water-loaded  plate 

The  convection  properties  of  the  turbulence  pressure  field  relate 
the  circular  frequency  w  and  the  turbulence  wave  number  k^.  by  the  equation 

JCj.  U.  =  m  (1) 

The  approximation  obtains  because,  in  actuality,  a  range  of  wave  numbers 
contributes  to  the  frequency  spectral  density  of  the  turbulent  pressures 
(e.g.,  see  Equation  (B60)  of  Reference  1).  Alternatively  stated,  for  a 
turbulent  field,  a  particular  wave  number  component  is  generally  associated 
with  a  distribution  of  frequencies  and/or  convection  speeds. 

Following  Dyer  (see  page  18  of  Reference  1),  the  hydrodynamic  coin¬ 
cidence  speed  vq  is  defined  as  the  speed  at  which  the  magnitude  of  the  mean 
flow  convection  velocity  U  along  the  flow  direction  is  equal  to  the  trace 

c  * 

speed  of  the  bending  wave  Cg  in  the  direction  of  flow,  i.e.,  v  =  vq  =  Cg. 
The  frequency  at  which  this  occurs  is  a  significant  parameter  because  it 
represents  a  value  at  which  we  can  expect  a  large  vibratory  response  due  to 
increased  amplitudes  of  the  modal  forcing  function  for  the  panel.  This 
parameter  is  called  the  hydrodynamic  coincidence  frequency  co^  and  is  given 
by  the  equation 

w,  =  k  CD  =  k  U  (2) 

n  p  B  pc  ' 

Thus  from  Equations  (1)  and  (2) ,  kT  =  k^  =  «h/Uc  at  the  hydrodynamic  coin¬ 
cidence  frequency. 


* 

The  bending  wave  is  also  referred  to  as  the  free  flexural  phase  velocity 
and  free  plate  bending  velocity. 
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Here 


k  = 
P 


because  (see  pages  18  and  28  of  Reference  1) 


ta 


2  /Bv4 


Cg  =  -j— -  =  (a  (— )  ;  w  =  at  hydrodynamic  coincidence 
P 


(3) 


(4) 


and 


(5) 


The  region  within  which  hydrodynamic  coincidence  can  exist  is  defined 

2  2  2 

by  a  semicircle  in  k  space  represented  by  (k^  -  Uc/2ic  c^)  +  =  (Uc/2k  c^). 

For  the  plate  coordinate  system  shown  in  Figure  3,  the  corresponding  classi¬ 
fication  of  plate  modes  in  wave  number  space  is  shown  in  Figure  4.  For  low- 
convection  velocities  U  «  CD,  the  locus  of  the  hydrodynamic  coincidence 
(HC)  curve  collapses  to  the  origin  and  the  response  is  composed  entirely  of 
hydrodynamically  fast  (IIF)  modes. 

From  these  equations  we  obtain  the  significant  relationship  from 
which  we  can  determine  for  a  given  plate  material  and  value  of  Uc, 
namely. 


CR  U 

“h  =  ^y/2  (6) 

The  hydrodynamic  critical  frequency  is  now  defined  as  the  resonance 
frequency  of  the  plate  mode  with  bending  wave  velocity  equal  to  the  con¬ 
vection  velocity.  Table  5a  presents  values  of  the  hydrodynamic  critical 
* 

frequency  f^  =  oo^/2tt  computed  for  the  fluid- loaded,  simply  supported 
aluminum  and  steel  plates  for  values  of  Uc  used  in  the  previous  computations 


* 

The  condition  of  hydrodynamic  coincidence  has  been  used  by  Maestrello 

2  2 

in  obtaining  computer  solutions  Y  (w)/p  (see  Appendix  H). 
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k  {k,,iul  SPACE 

kT  uc  ' 


in  the  nodal  lattice, 

the  projection  of  k  on  k,  is  k  =  a~/l  and 
r  J  >3n  1  ra  1 

the  projection  of  k  on  k_  is  k  =  n;/i.. 

2.  For  an  infinite  plate  the  curves  in 

{kpw}  space  show  that  above  wh  hydrodynamic 

coincidence  does  not  occur  for  any  wavenumber. 

At  very  high  frequencies  the  sonic  curve  intersects 

the  free  wave  curve  because  —  increases  faster 

c 


than 


and  acoustic  coincidence  occurs. 


For  additional  details  on  corresponding 
relationships  in  {Kj ,K.}space  see  Bibliography 

(White,  P.H. ) 


Figure  4  -  Classification  of  Modes  in  Wave  Number  Space 
and  Turbulence  Vibroacoustic  Relationships 
in  Wave  Number-Frequency  Space 
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TABLE  5 


Computation  of  Hydrodynamic  Critical  Frequencies  and  Acoustic 
Critical  Frequencies  for  Fluid-Loaded,  Simply  Supported  Alum¬ 
inum  and  Steel  Plates 
TABLE  5a 

Computation  of  Hydrodynamic  Critical  Frequencies 
for  Aluminum  and  Steel  Plates 


uc 

Uc2 

fh  - 

uc2 

vh  ■ 

(ft/sec) 

alum,  plate 

steel  plate 

8 

0.62 

0.05 

16 

2.48 

0.20 

32 

9.95 

0.79 

64 

39.87 

3.17 

ALUMINUM:  <  =  9.64  x  10  ft 
c„  =  17,000  ft/sec 

A 

Steel  :  k  =  1.204  x  10'2  ft 
c.  =  17,000  ft/sec 

cwater  :  =  £->000  ft/sec 


TABLE  5b 

Computation  of  Acoustic 
Critical  Frequencies 
for  Aluminum  and  Steel 
Plates 


Material 

fc  ■  <H2> 

Aluminum 

243,000 

Steel 

19,400 

of  the  plate  response.  The  results  were  obtained  by  use  of  liquation  (6). 

2  2  2  2  2 

For  a  finite  plate  at  resonance,  k^  =  kg  r  =  [(mTt/a)  +  (rnt/b)  ]. 
For  air  (or  in  vacuo )  loading 


mn 


-  k  c„ 


(7) 


2 

the  equation  w  =  K  c£  (m7I/a)  =  (mu/ a)  uc  represents  the  hydrodynamic 

coincidence  condition  for  the  finite  plate  at  resonance  because 'U  is  the 
convection  velocity  in  the  flow  direction  only  (see  page  132)  and 


w  w 

_  mn  ..  mn  „  „ 

CB  =  IT"  or  Uc  =  —  =  K  C£  kmn 


mn 


K  C, 


(?) 


at  to  = 


wh 


(8) 


For  water  loading  (see  liquation  (113)) 


0)  -»  w  =  w 
mn  mn 


a  p 


I 

-i"2 


1  + 


P„  kh  fl-c.V/2 


s  s 


(9) 


The  radiative  properties  of  the  plate  immersed  in  a  fluid  relate  the 
given  frequency  w  (Equation  (1))  and  the  acoustic  wave  number  k  by  the 
equation: 


kc  =  us  (10) 

From  Equations  (1)  and  (10),  we  obtain: 

II  U 

k  =  ~  kT  =  0.8  kT  •  —  =  0.8  k?  •  Mach  No.  (11) 

These  last  two  relations  are  given  to  add  to  the  theoretical  picture.  How¬ 
ever,  no  computations  for  the  radiation'  of  plates  were  made  for  the  present 
report . 

On  page  24  of  Reference  1,  we  defined  the  sound  coincidence  or  cutoff 
or  critical  frequency  us  =  us^  as  corresponding  to  c  =  C^,  i.e.,  the  fre¬ 
quency  at  which  the  flexural  wavelength  of  a  thin  plate  equals  the  acoustic 
★ 

wavelength  in  water.  At  this  frequency  we  expect  a  greatly  increased 
acoustic  response. 


* 

The  critical  frequency  is  used  for  a  resonant  plate. 
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R 


1 


Then 


Moreover  since  at  this  frequency  (see  page  27  of  Reference  1) 


1 


Table  5b  presents  values  of  f£  =  W./2TT  computed  for  the  fluid-loaded,  simply 
supported  aluminum  and  steel  plates.  The  results  were  obtained  by  use  of 
Equation  (13) .  Obviously  there  is  one  and  only  one  sound  coincidence  (or 
critical  frequency)  for  each  plate  material. 

According  to  Leehey  (see  Bibliography),  the  frequency  range  of 
interest  for  sonar  self-noise  application  is 


U2 
-  c 

“■ft  ”  K'  C„ 


<  to  <  to 


2 

c 

IC 


Hydrodynamic  coincident  effects  are  not  important  for  this  range  (see 
Tables  2  and  3),  and  the  radiated  power  per  unit  area  is  associated  with 
both  the  decay  of  turbulence  and  the  presence  of  plate  boundaries. 

We  observe  that  the  near  field  in  the  vicinity  of  the  panel 
represents  the  predominant  portion  of  the  transduced  pressure  spectrum. 

This  spectrum  is  composed  of  nonpropagating  components  caused  by  hydro- 
dynamic  coincidence  effects.  Thus,  the  wave  numbers  chiefly  in  evidence  are 
those  in  the  immediate  vicinity  of  the  hydrodynamic  coincidence  curve  which 
corresponds  to  an  excitation  field  progressing  in  the  flow  direction  with 
velocity  Uc  and  frequency  w;  see  Figure  4.  The  spectrum  of  boundary  layer 
pressure  is  distributed  about  =  to/Uc  with  most  of  the  energy  at 

rather  iow  values  of  k^.  As  we  proceed  into  the  far  field,  the  near-field 
components  decay  and  only  the  wave  numbers  at  or  below  the  sonic  line, 
represented  in  Figure  4  by  a  circle  k  =  k0  =  w/c  is  present. 

We  now  show  that  at  sufficiently  low  wave  numbers,  the  added  mass  (or 
near  field)  of  a  fluid-loaded  plate  subject  to  turbulence  excitation  may 
contribute  significantly  to  the  vibroacoustic  response  of  a  structure  in  its 
vicinity  and,  in  turn,  may  also  be  significantly  affected  by  that  response. 


I' 


s! 


T 


jS! 


i 
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The  added  mass  or  fluid  inertial  loading  of  a  flexurally  vibrating 
plate  represents  the  imaginary  or  reactive  part  of  the  impedance  associated 
with  the  reaction  of  the  fluid  to  the  vibrating  plate.  The  corresponding 

fluid  pressure  is  considered  to  be  effective  over  a  distance  d-s?-^-  from  the 

A  b 

plate;  the  distance  d  =  -g-  is  called  the  near  field.  Thus,  at  sufficiently 

small  wave  numbers  (d  <  A/6  or  kd  <  1),  the  near  field  will  exert  a  pressure 
on  an  adjacent  structure  lying  within  the  effective  bounds  of  the  field. 

If  the  structure  is  flexible,  we  will  then  have  a  complex  coupled  vibro- 
acoustic  system  involving  the  plate,  the  adjacent  structure,  and  the  inter¬ 
vening  fluid  medium.  The  interaction  would  in  general  couple  the  near-  and 
far- field  (or  radiation)  pressures  on  both  plate  and  adjacent  structure  and 
would  include  the  phenomena  of  reflections,  scattering,  etc.  Wo  explain 
this  by  recognizing  that  each  flexible  body  or  source  works  against  its  own 
sound  pressure  which  represents  the  reaction  of  the  medium  to  its  motion  as 
well  as  against  the  sound  pressure  that  is  generated  by  the  adjacent  source 
(e.g.,  the  opposing  body).  Thus,  two  sound  sources  (i.e.,  flexible  bodies) 
in  close  proximity  (d  <  -|-)  react  with  each  other  and  the  sources  may 
generate  considerably  more  sound  energy  than  if  they  were  further  apart. 

Sound  sources  of  this  type  in  close  proximity  are  said  to  be  dependent,  and 
the  power  that  each  source  generates  individually  cannot  simply  be  added  to 
yield  the  actual  or  total  power  of  the  system  which  is  due  to  the  interaction 
effects  associated  with  the  near  field. 


CONCLUSIONS 


The  chief  conclusions  drawn  from  this  investigation  are: 

1.  The  added  mass  and  corresponding  natural  frequency  of  a  fluid- 
loaded  rectangular  plate  are  more  significant  for  the  uncoupled  modes  and 
are  easily  computed  using  the  results  given  in  Table  1.  For  a  first 
approximation,  only  the  results  for  uncoupled  modes  need  to  be  considered 
in  vibroacoustic  computations. 

2.  The  added  mass  results  for  the  coupled  modes  can  be  computed 
(using  the  results  in  Table  1)  and  added  to  the  results  for  the  uncoupled 
modes  to  refine  the  accuracy  of  the  computation  for  the  natural  frequency  or 
to  determine  the  effect  of  these  modes  on  the  natural  frequency  and  vibro- 
acoustic  response. 

3.  Separate  results  for  the  added  mass  and  natural  frequencies  of 
only  the  uncoupled  modes  can  be  computed  (using  the  results  in  Table  1)  to 
determine  the  absolute  numerical  contribution  of  these  modes  to  the  vibro¬ 
acoustic  response  or  to  identify  a  corresponding  response. 

4.  The  methods  of  analysis  used  in  Appendixes  A-E  yield  the  same 
added  mass  and  natural  frequency  results  for  uncoupled  modes  and  are 
applicable,  in  their  essence,  to  both  high  and  low  wave  numbers  and  for  low 
and  high  frequencies  (see  the  remarks  column  of  Table  1  relative  to 
Appendix  B;  see  also  Appendix  E) . 

5.  The  Feit-Junger  method  of  analysis  used  in  Appendix  A  can  be  ex¬ 
tended  by  the  methods  of  Leibowitz  presented  in  Appendix  C  to  yield  the  same 
added  mass  and  natural  frequency  results  for  the  coupled  modes  as  obtained 
by  Davies  (Appendix  B) . 

6.  The  results  obtained  in  Appendixes  A-E  for  the  added  mass  and 
natural  frequency  of  a  fluid- loaded,  simply  supported  plate  can  be  extended 
to  yield  corresponding  results  for  a  fluid- loaded  clajnped-clamped  plate  by 
using  the  Leibowitz-Wallace  methods  given  in  Reference  3.  The  results  ob¬ 
tained  by  Greenspcn  and  Leibowitz  in  Appendixes  F  and  G  can  be  used 
directly  to  compute  the  added  mass  and  natural  frequency  of  either  a  simply 
supported  or  clajnped-clamped  plate;  the  results  for  Appendix  G  are 
particularly  applicable  to  the  low  frequency  response.  The  Greenspon- 
Leibowitz  results  are  considered  to  be  less  accurate  than  the  former 
results. 
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7.  The  total  mass  and  corresponding  natural  frequency  for  a  plate 
radiating  into  free  (half)  space  as  well  as  into  a  closed  rectangular 
space  can  be  computed  using  the  results  presented  in  Table  1  and  the  results 
obtained  from  Appendix  A  of  Reference  1  given  in  the  Discussion. 

8.  The  computed  results  for  both  the  aluminum  and  steel  plates  show 
that  the  contribution  of  the  higher  modes  to  the  total  vibratory  response 
is  not  negligible  e.g.,  for  the  aluminum  plate,  the  magnitudes  of  the  root 
mean  square  displacement  for  the  (1,4)  and  (2,4)  modes  are  of  the  same 
order  as  that  for  the  (1,3)  and  (2,3)  modes  for  a  given  value  of  con¬ 
vection  velocity  and  damping.  Thus,  determination  of  the  total  vibratory 

-  displacement  requires  that  the  computations  include  the  contributions  of  the 

several  modes  of  vibration  deemed  to  be  significant. 

9.  Turbulence-induced  plate  vibration  may  be  of  significant  magni¬ 
tude  for  small  damping  and/or  sufficiently  large  mean  square  pressure 
fluctuations;  the  magnitude  is  enhanced,  but  not  radically,  at  sufficiently 
large  values  of  q  =  k/ks.  Hence  in  computation  it  appears  practical  to  use 
a  single  representative  value  for  q  for  the  range  of  turbulence  frequencies 
distributed  about  co  =  Uc  kT  =  kc  corresponding  to  0  <  q  =  k/ks  <1.0.  A 
practical  alternative  would  be  to  select  the  average  of  the  mean  square 
displacement  responses  computed  for  small  and  large  values  of  q. 

RECOMMENDATIONS 

To  simplify  the  computational  procedure  (and  the  computer  program) 
and  to  achieve  reasonably  accurate  vibroacoustic  results  for  a  vibrating 
plate  fluid  loaded  on  one  side,  the  following  recommendations  are  made. 

Note,  however,  that  the  user  who  wishes  to  refine  the  accuracy  of  compu- 
•  tations,  determine  the  coupled  mode  contribution,  identify  a  coupled  mode 

response,  or  treat  a  closed  rectangular  cavity  can  incorporate  the 
additional  relevant  results  presented  here  into  the  computer  program. 

1.  It  is  recommended,  as  a  first  approximation,  that  the  equations 
for  the  added  mass  and  corresponding  natural  frequency  of  only  the 
uncoupled  or  self  (dominant)  modes  of  a  fluid- loaded,  simply  supported 
rectangular  plate  be  used  in  making  vibroacoustic  computations.  This  in¬ 
cludes  the  uncoupled  mode  equation  common  to  all  of  the  results  obtained  in 
Appendixes  A-E;  see  Appendix  H  for  the  corresponding  computer  program. 
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2.  It  is  recommended  that  vibroacoustic  computations  for  fluid- 
loaded,  clamped-elamped  rectangular  plates  be  made  (for  uncoupled  modes 
only)  by  extending  the  equation  for  a  fluid- loaded,  simply  supported 
rectangular  plate  to  include  this  case,  in  accordance  with  the  Leibowitz- 
Wallace  methods  of  Reference  3;  see  Appendix  H  for  the  corresponding  com¬ 
puter  program. 
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APPENDIX  A 

THE  FEIT-JUNGER  METHOD 


NOTATION 


A 

c 

c 

P 


D 

E 

e 


F(x,y) 

F 

mn 


h 

I 

mnpq 


i 


k  ,k  ,k  ,k 
nr  n’  p’  q 


k 

s 


L 


y 


p. 


q 


m 

mnpq 

p(x,  y,  z  >  0) 
p(x,  y,  z  =  0) 
P(Yx>  Yy 
r 

mn 

r 

mnpq 


Amplitude  of  p 

Velocity  of  sound  in  fluid  medium 
Compressicnal  wave  velocity  of  the  plate  equal  to 


fcfcf 


Flexural  rigidity  of  plate  equal  to 


2  ,3 

c  p  h 
P  s 

12 


Eh' 


12  (1-V  ) 


Young’s  modulus 

Equal  to  2.718;  base  for  natural  or  Naperian  system  of 
logarithms 

Driving  force  applied  at  coordinates  x,  y,  z  =  0 
Generalized  force  for  the  mn  mode 


Plate  thickness 

Defined  by  Equation  (A17b) 

Equal  to  fCT 

Acoustic  wave  number  equal  to  w/c 

Modal  wave  numbers  defined  by  Equation  (A3) 

2  2  1/2 

Surface  wave  number  equal  to  (k  +  k  ) 

n  m  n 

Half  length  and  half  width  of  plate,  respectively 
Total  plate  structural  mass 
Mode  numbers 

Added  modal  mass  per  unit  area 

Added  mass  of  coupled  modes  mnpq  per  unit  area 

Pressure  in  fluid 

Pressure  on  surface  of  plate 

Double  inverse  Fourier  transform  in  Y  and  Y 

x  y 

Radiation  modal  damping 

Mutual  specific  acoustic  resistance  for  coupled  modes  mnpq 
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Structural  modal  resistance 


w  ,  w 
mn  pq 

w(x,y,  z=0) 


w(Yx,  Yy) 


w  (Y  ,  Y  ) 
pq  X*  yJ 

x,  y,  z 


SOc-Yj) 


O) 

(A) 

c 

w  ,  w 
mn  mn 


Displacement  amplitude  of  vibration  of  a  plate  for  the  mn 
and  pq  modes  of  vibration,  respectively 

Modal  displacement  at  surface  of  plate 

Fluid  particle  displacement  at  surface  of  plate  equal  to 
plate  displacement 

Transform  of  displacement  w(x,  y,  z  =  0);  equal  to  a 
series  of  modal  transforms 

Modal  transform  of  displacement  at  the  surface  of  plate 

Rectangular  coordinates;  x  and  y  are  in  the  plane  of  the 
plate  and  z  is  normal  to  the  plate 

2  2  2 

n  ,  *  3  3  3 

Equal  to  — =•  +  — r  +  — j 

3x  3y  3z 

/ 

Delta  function;  1 6 (k .  -  Y-)d„  = 

J  i  J  Y. 

Wave  numbers  which  are  the  coordinates  in  Fourier 
transform  space 

A  small  quantity 

Structural  loss  factor  equal  to 

Poisson's  ratio 
Mass  density  of  fluid  medium 
Added  mass  density  of  fluid  per 

Sum  of  mass  densities  of  plate  and  fluid  equal  to  p  +p 
Mass  density  of  the  plate 


/  2  , 
r  w/u)  p  h 
s  /  mn  s 

unit  volume  equal  to  m  /h 


1,  y.  =  k. 

J  l 


0,  Y-  ^  k. 

j  l 


Natural  circular  frequency  of  vibration 


Coincidence  frequency  equal  to  c 


In  vaouo  and  submerged  natural  circular  frequency, 
respectively,  for  the  mn  mode  of  vibration  i.e.,  resonance 
frequency  for  the  in  vaouo  and  submerged  plate, 
respectively 
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DERIVATION 


From  Reference  4,  the  even  modes  of  a  simply  supported  rectangular 

* 

plate  vibrating  in  vacuo  have  a  configuration  described  by 


w  (x,y)  =  W  cos  k  x  cos  k  y,  I x I <  L  ,  |y|<  L 
mn  mn  m  nJ  ’  1  1  x’  1/1  y 


=  0 


with  boundary  conditions 


M*  Lx»  lyls  Ly 


(Al) 


cos  k  L  =  cos  k  L  =  01 
m  x  n  y 


m 


sin  k  L  =  (-1) 
m  x  v  3 

sin  k  L  =  (-l)n 
n  y  v  3 

The  boundary  conditions  restrict  the  wave  numbers  k  ,  k  ,  or  k  to  the 
3  m*  n  s 

values 


(A2) 


km  Lx  =  (2m  +  l)~  ;  m,n  =  0,  1,  2* 

n 

1/2 


n  y  n 

1/2 

ks  ■  (C *  K2)1'2  -- t  [(-rf)2  *  fr1)2]  !  »• 2- 


(A3) 


Also  k  ^  >  k^  for  w  <  w  =  c2//h/l/T2)c  (see  Reference  4)  so  that  each  wave 
s  mn  c  \  r  /  p 

is  characterized  by  one  pair  of  wave  numbers  k^  and  k^,  i.e.,  the  normal  modes 
of  a  plate  vibrating  in  vacuo  are  described  in  terms  of  a  discrete  wave 
number  spectrum.  The  modal  configurations  are  defined  to  be  orthogonal  to 
each  other  so  that  each  mode  can  be  excited  independently  by  a  suitable 
distribution  of  the  load. 

For  the  submerged  plate  (exposed  to  water  on  one  side),  however,  each 
of  the  originally  normal  (km,  k  )  modes  generates  an  acoustic  pressure  in 
the  plane  of  the  plate  (which  can  be  represented  as  an  inverse  Fourier 


★ 

This  configuration  is  of  practical  importance  because  it  matches  the 
modes  of  vibration  of  a  simply  supported,  rectangular  plate  driven  at  the 
center. 
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b- 


>1 

a 


a 

JS 


transform)  whose  wave  number  spectrum  is  continuous  and  thereby  encompasses 

the  discrete  wave  numbers  of  the  other  originally  (i.e.,  in  vacuo )  normal 

modes.  The  resulting  pressure  distribution  due  to  a  single  mode  is  not 

orthogonal  to  the  other  modes.  This  causes  the  modal  configurations  to 

couple,  and  thereby  to  lose  their  normal  mode  character.  Thus  if  we 

attempt  to  formulate  the  problem  in  terms  of  in  vacuo  normal  modes,  we  shall 

find  that  these  modes  become  coupled;  the  subsequent  discussion  will  clarify 

these  features.  The  modes  also  lose  their  standing  wave  character  in  the 
2  2 

range  k  >  k  because  the  plate  boundaries  are  energy  sinks. 

S 

For  a  distributed  load,  the  forced  equation  of  motion  for  the  (m,n) 

mode  of  an  undamped  thin  rectangular  plate  extending  from  -L  to  L  and  -L 
+  .  .  4-6  x  x  y 

to  L  is 

y 


■-1 

ft 

I 


L  L  p  h  L2  -  a)2 1  W 
x  y  s  mn  J  i 


=  F 
mn  mn 


(A4) 


Here,  L  L  p  h  =  M  /4  where  M  is  the  mass  of  the  plate, 
x  y  s  P  P 

he  2  2  /  d  2  ? 

oj  =  — E.  fk  +  k  ")  =  (  — r  |  [k  +  k  “]  is  the  in  vacuo  natural  fre- 
mn  .rrz  m  n  \  p_h  /  1  m  n  J 


quency  of  this  particular  mode,  Fmn  is  the  generalized  force  associated 
with  a  concentrated  driving  force  applied  at  x=0,  y=0,  and  the  surface 
pressure  p(x,y,0)  is  represented  by 


F  =  F(0 
mn 


„  L  L 

r  X  r  y 

,0)  -  I  I  p(x,y,0)  cos  km  x  cos  kRydxdy 


(A5) 


-L  -L 
x  y 


To  obtain  an  explicit  expression  for  p(x,y,0)  as  a  double  integral,  the 

surface  pressure,  which  is  spatially  aperiodic  in  x  and  y,  is  written  as  a 

double  inverse  transform  in  Y  and  y  . 

x  'y 


With  inclusion  of  a  structural  loss  factor 


h  c 


0) 


mn 


yn 


r,  2  ,2. 
+  K  1 

m  n 


(-?) 


(see  Reference  4) . 
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p(y. 


n[-i(v  x  +  Yy)] 
p(x,y,z)  e  y 


dx  dy 


(A6) 


-OO  _CO 


P(x,  y;  »)  =  ^2  J  J  PCyx,yy;  *)  e'  CYx  Yyy)1dYx  ^Yy  (*?) 


-OO  -OO 


Now  the  three-dimensional  Helmholtz  wave  equation  is 


L  3x  3y  3z  J 


z)  =  [V2  +  k2]  p(x,y;z)  =  0  (A8) 


The  double  transform  of  the  foregoing  equation  is 


(' 


2  _  y  2  _  y  2  +  3  \  CY  Y  z)  =  o 

x  y  3z2/  x  y 


(A9) 


xctz 

Assuming  a  solution  of  the  form  p  =  Ae  ,  the  solution  of  Equation  (A9)  is 

~  ifk2-Y  2-Y  z 

pCyy’Yv’2^  =  A®11  x  y  J  caio) 

A  y 

2 

where  the  boundary  condition  3p/3z  =  poo  w  at  z=0  yields  the  value  (see 
Chapter  IV  of  Reference  4) 


*  3 

The  Fourier  transform  of  the  term — j  p(x,y;z)  is  obtained  by  inte- 

3x 


limits  x=+  «.  This  yields 


grating  by  parts  and  setting  p  and  its  derivative  equal  to  zero  at  the 

9  -iYx  x  ? 

e  dx  =  -  p(Yx,y;z).  Also 

2  2  2  2 

because  3  /3y  ,  3  /3z  ,  and  k2  are  independent  of  x,  they  can  be  taken  out- 

22  22  2  °°  -  i  y  x 

side  the  integral  sign  so  that  (3  /3y  +  3  /3x  +  k  ) f  p(x,y;z)e  x  dx 


-  OO 

j  32p/3x^ 

V-  OO 


j. 


2  2  2  2  2 
(3  /3y  +  3  /3z  +  k  )  p(Y  ,y;z).  The  procedure  is  now  repeated  for 

2  x 
ri  ~ 

—  p(Y  ,y;z)  to  obtain  Equation  (A9) . 

c.  X 

3y 
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(All) 


Substituting  the  time  integral  of  the  first  of  Equations  (A15)  in 
Equation  (A14) ,  we  get 


Fmn  =  F(M)  +  >^ff  - -  -  i/2~  dYxdYy  (A16) 

n  4tt  W  JJ  2  v  2  v  2.1/2  X  y 

mn  pq  (k  -  Yx  -  Y  ) 


Using  the  time  integral  of  the  second  of  Equations  (A15)  for  both  w  and 
~  mn 

w  ,  we  get 

pq 


F  =  F(0,0)  + 
mn  5 


i4pc2k2k  k  (-l)m+nYl  V  ff 

: - -  2l2l  W  k  k  (-l)P+q 

J  p  q  pq  p  q  J  J 


cos2  Y  L  cos2  Y  L  dY  dY 
xx  y  y  x  y 


(A17) 


7  7  7'??  7  7  7  7  7  7 

(k  -Y  -Y  )  (k  -Y  )(k  -Y  )(k  -Y  (k  -Y  ) 
x  y  m  x^^n  y  p  x-'  q  y 


The  exact  solution  of  Equation  (A17)  requires  a  numerical  integration  of 
the  branch  cut  integral.  Reference  4  avoids  the  determination  of  the  exact 
solution  by  evaluating  the  integral  for  the  high  wave  numbers  (short  wave¬ 
length)  limit  only.  This  simplifies  the  analysis,  and  the  final  expression 
is  considered  to  hold  for  a  finite  rectangular  plate  over  a  particular  wave 
number  and  frequency  range. 


Write 


F  =  F  (0 
mn 


,0)  +  ia)  I 


W 

mnpq  pq 


(A17a) 


where  we  define 


4pwk  k  (-1) 
m  nv 


z 


(-l)p+q  k  k 

N  '  r 


(A17b) 


(equation  continued  on  page  48) 


(Equation  continued  from  page  47) 


c°s2  U  cos2  Y  L  dY  dY 

XX  y  y  x  v 


x  x  y  y  w  x  y 

4 


(A17b) 

Cont’d 


r  -  icon  iL  L 
(  mnpq  mnpq  /  x  y 


Thus  the  complex  quantity  1^^  is  the  product  of  L^Ly  and  the  sum  of  the 

mutual  specific  acoustic  resistance  r  and  the  reactance  -  icom  ,  also 

mnpq  mnpq 

referred  to  as  modal  coupling  coefficients,  where 


r  =  Re  (I  ) 
mnpq  mnpq 


m  - - Imag  (I  ) 

mnpq  to  mnpq 


F  =  F (0,0)  + 
mn  ’  J 


ico  x  ,  L  L  /r  -ico  \  W 
“J  xy\  mnpq  mnpq j  p 


=  F (0,0)  +  i 


icoL  L  ( i 

xy\ 


r  -icom  j  W  + 
,  mn  mn/  mn 


Z 


pq/mn 


icoL  L  (r  -  icom  \  W 
x  y  \  mnpq  mnpq/  pq 


=  F(0,0)  +  L  L  (co  m  +  icor  )W  + 
x  y  \  mn  mn  /  mn 


Z  (w‘ 

pq/mn 


m  +  icor  \  W 
mnpq  mnpq  /  pq 


For  the  simplified  analysis,  we  write  the  Lagrange  equation  for  the 
forced  motion  of  a  mode,  replacing  Equations  (A4)  and  (A5)  by  (note  that 
harmonic  time  dependence  is  in  the  form  e  iart) 


L  L  p  h  [co2  -  co2 1  W  -ico 
x  y  s  L  ran  J  mn 


/L  L  p 

x  y  ■ 

\  “ 


,  2 

hw  n 
s  mn  s 


(A18) 


F(0,0)+L  L  fco2m  +  icor  \  W  +  Zi*  L  L  ( co2in  +  icor  ^  W  =  F 

x  y  V  mn  mn/  mn  ^  x  y  V  mnpq  mnpq/  pq  mn 
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where  the  left  member  now  includes  a  term  for  the  structural  resistance 

9 

r  =  u“  p  hn  /(i)  expressed  in  terns  of  a  structural  loss  factor  n  and  the 
s  ran  s  s  ‘  s 

radiation  loading  portion  of  the  generalized  forces  are  written  as  two 
terms.  The  first  of  these  terms  represents  the  uncoupled  (added  mass  and 
radiation)  damping,  i.e.,  the  self- impedance  due  to  radiation  loading 
whereas  the  second  of  these  terms  represents  the  coupled  added  mass  and 
radiation  damping  i.e.,  the  mutual  impedance  due  to  radiation  loading.  If 
included  in  the  computation,  the  coupled  terms  must  be  evaluated  numerically. 

Transferring  the  radiation- loading  portion  of  the  generalized  forces 
to  the  left  member  of  Equation  (A18) ,  the  Lagrange  equation  becomes 


f  2  /  ,  \  .  /^s^mn  ns  \ 

(psh  ♦  mnm  /  "ltl)  y - - -  *  rm„/ 


+  0)  p  h 
mn  s 


(A19) 


-  Z 

pq^mn 


i  +  iarr 
mnpq  mnpq 


which  is  a  doubly  infinite  set  of  equations  for  the  unknowns  W  coupled  by 

the  terms  m  and  r  .  Thus  I  is  a  coupling  coefficient  linking  an 
mnpq  mnpq  mnpq  r  t>  t> 

(m,n)  mode  to  a  (p,q)  mode. 

For  kL  ,  kL  >3  which  for  w  «  w  ,  is  equivalent  to  k  L  ,  k  L  »  1 
x’  y  c’  n  m  x’  n  y 

(i.e.,  the  criteria  for  large  plates),  see  Equation  (A3),  the  coupling  terms 
are  much  smaller  than  the  self-impedance  components  so  that  they  can  be 
ignored  in  making  an  approximate  evaluation  of  the  far  field. 

Feit  et  al.  now  proceed  to  solve  the  integral  Equation  (A14)  for 
the  case  kfflLx  and  k^L^  approaching  infinity  in  order  to  gain  some  insight 
into  the  above  approximation.  In  this  case  the  orthogonality  of  the 
cosines  in  Equation  (A14)  yields 


However,  even  without  going  to  the  limit  of  infinite  k  L  ,  K  L  ,  the  pq 
’  o  o  m  x  n  y  r  1 

summation  is  small  when  k  L  ,k  L  is  large;  see  Reference  7. 

m  x’  n  y  & 
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L  L 
y  r  x 


lim  f  f  cos  Y  x  cos  k  x  cos  Y  v  cos  k  y  dxdy  =  4n2  6(k  -Y  )S(k  -Y  ) 

II  x  m  y  ir  J  v  m  xJ  y  n  yJ 

kL  +»  J-L  -L 


m  x 

k  L  -*  00 
n  y 


v  x 


(A20) 


Using  these  6  functions  in  Equation  (A14) ,  we  find 


F  =  F(0,0)  +  ipc2k2 
mn 


CO 

Z  JJ 

pq 


w  (Y  ,Y  )6(k  -Y  )6(k  -Y  )dY  dY  (A21) 
pq  x  y  m  x  n  y  x  y 

(k2-Y  2-Y  2)1/2 
x  y 


or 


F  =  F(0,0)  + 
mn  ^  J 


ipc2k2  /  ,  w  (k  ,k 
pqv  m*  n 
pq  ^ 


(k2-k  2-k  Wn 

a  n  ' 


Using  the  second  of  Equations  (A15) ,  Equation  (A22)  becomes :  (In 
Equation  (A15),  let  m-*p,  n-*q,  x*Hn,  y^n.) 

i4pc2  /  ,  W  k  k  (-1)P+C*  cos  k  L  cos  k  L 

— ^  pq  p  q  m  x  n  y 


F  =  F(0,0)  + 
mn  ^  ’  ’ 


(k2-k  2-k  2)1/2(k  2-k  2)(k  2-k  2) 
v  m  n  J  p  m  '  ^  q  n 


(A22) 


(A23) 


The  second  term  in  the  right  member  vanishes  when  kp^km=Yx> 

k  jik  =Y  (see  Appendix  B) .  For  k  =k  =Y  ,  k  =k  =Y  ,  the  denominator  vanishes 
q  n  y  v  p  m  x’  q  n  y 

but  the  second  term  is  finite  because  the  boundary  conditions.  Equation  (A3), 
require  that  the  numerator  also  vanish.  Hence  we  evaluate  the  indeterminate 
quantity 


lim 
k  -k 


p  m 
k  -k 


q  n 


rcos  k  L  — i 

rcos  k  L  ~i 

+  0 

m  x 

n  y 

■*  0 

(k  2-k  2) 
Lv  p  m  -J 
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f 


For  lira  cos  k  L  ?  7 

k  -k  0  - —  let  k  =  k  +  e,  k  “  «  k  ~  +  2ek 

pffi  ■7'’  PQ  P  m  m 

k  “-k  “ 
p  m 


cos  k  L  cos  (R  -e)L 

my  , .  p  > 

lim  - -  - - - - 

V*.  -  0  k  2-k  2  6+0  2ek 

P  1,1  pm  m 


Similarly, 


cos  k  L  cos  eL  +  sin  k  L  sin  eL 
m  x  x  p  x  x 


(-DpeLx  C-DPI.X 


2ek  2k 

m  m 


(A24) 


k  -k  -* 0  k  -k 
<1  n  q  n 


cos  k  L  (-1)C,L 

ny/  y 

,  2  ,  2  *  2k 


(A25) 


Substituting  Equations  (A24)  and  (A25)  in  Equation  (A23) ,  setting 
(-l)"P+2<*  =  1,  letting  W  ->  W  ,  and  dropping  2  since  the  coupling  terms 

Pq  mn  pq 

vanish  in  the  large  k  L  ,  k  L  limit,  we  obtain 
b  m  x  n  y 


F  =  F  (0,0)  + 
mn  ^  J 


ipc2k2  W  L  L 
ran  x  / 

(k2-k  2-k  V/2 

v  m  n 


(A26) 


Hence  in  the  large  k  L  ,  k  L  limit 
b  m  x  n  y 


(A27) 


F  =  F(0,0)  + 
mn 


u  p  LxLy  "mn 

(K  2*k  W)1'2 
m  n  J 


=  F(0,0)  +  u  L  L  m  K 
x  y  mn  mn 


2  2  2 
(k  z  +  k  n  > 

m  n  ' 


where  (see  Equation  (A18)) 


(k  2+k  2-k2)1/2 
v  m  n  J 


_ _  P_ 

(ks2-k2)'/2  ks 


(A27a) 


is  the  added  mass  per  unit  area  which  is  shown  in  Reference  4  to  correspond 
to  an  infinite  train  of  straight-crested  parallel  waves  or  of  straight- 
crested  orthogonal  waves.  Feit  and  Junger  assume  that  this  added  mass  will 
hold  for  a  finite  rectangular  plate,  i.e.,  the  added  mass  of  a  finite  plate 
whose  dynamic  configuration  embodies  many  nodal  lines  is  effectively  that 
of  an  infinite  train  of  standing  waves. 

For  the  low  frequency  range  where  few  modes  contribute  to  the  far 

field,  a  deterministic  approach  can  be  used.  The  resonance  frequency  for 

the  submerged  plate  a>  is  then  determined  as  follows: 

mn 


(kg2-k2)1/2h 


;  k  >  k 


(A28) 


Therefore 


p  =  p  +  p  =  p  +  -  ;  k  >  k 

e  s  Ka  Hs  2  2  1/2  s 

s  ' 


(A29) 


„  2  ,2.1/2, 

ps(ks  -k )  h 


(A30) 


APPENDIX  B 
THE  DAVIES  METHOD 


NOTATION 

Ap  Area  of  plate  equal  to 

Coefficient  of  defined  by  Equations  (B46)  and  (B47) 
cq  Velocity  of  sound  in  fluid  medium 

D  Flexural  rigidity 

e  Equal  to  2.718;  base  for  natural  or  Naperian  system  of 

logarithms 

I  (k Functions  defined  by  Equations  (B23)  and  (B26) , 
nr  respectively 

I^mCkj)  Function  defined  by  Equation  (B24) 

i  Equal  to  yn 

2  2 

k  Wave  number  equal  to  |k|=k^  +k^ 

k  Wave  number  vector  with  components  {k^k^l 

kQ  Acoustic  wave  number  equal  to  co/c^ 

kj  Component  of  k  lying  along  x.^  axis 

k^  Component  of  k  lying  along  x^  axis 

k  Wave  number  equal  to  imr /£., 

m  1 

k^  Wave  number  equal  to  nir/S.^ 

k  Surface  wave  number  equal  to  \J  k  2+k  2 

mn  n  y  m  n 

£^,^2  Length  and  width  of  plate,  respectively 

Mp  Total  mass  per  unit  area  represented  by  Equation  (B49) 

m^  Mass  per  unit  area  of  panel 

m  Added  mass  or  fluid  loading  per  unit  area 

mn  r 

m,  n,  q,  r  Mode  numbers 

P  Acoustic  pressure  generated  by  motion  of  panel 

Pmn  Acoustic  modal  pressure  defined  by  Equation  (B8) 

p  Pressure  field  driving  panel 

p  Driving  modal  pressure  defined  by  Equation  (B7) 


t 
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R  («) 
mnqr 

S  (k) 
mn 


Coupling  coefficient  connecting  the  m,n  mode  with  the 
q,r  mode 

Shape  function 

Modal  radiation  coupling  term 
Modal  mass  loading  coupling  term 


V  X2’V*'X2 


y(x,t) 


Panel  normal  velocity  displacement 
Modal  velocity  amplitude 

Rectangular  coordinates;  x is  normal  to  the  panel  and 
the  origin  is  at  one  comer  of  the  panel;  x  =  {xj,x^} 
Panel  normal  displacement 


Z(k,w) 
sgn  to 


«(k.-k.) 
1  y 


»  (x) 

mn  - 


Radiation  impedance 

Equal  to  -1  for  to  <  0,  +  ]  for  to  >  0 

Denotes  complex  conjugate 

Coefficient  accounting  for  mechanical  damping  of  panel 


Equal  to 


fa  y  23  y  3  y^ 

7  +  ~~z — \  +  — T  1  for  isotropic  plate 
^3x  3x  3y  3y4  / 


oo 

Delta  function:  (  I  6(k.-k.)dk.  =  i  f0r  k.=k. 

JJ  1  1  i  il 

=  0  for  k.#c. 

1  J 

Kronecker  delta  equal  to  1  for  i=j ,  equal  to  0  otherwise 

Modal  structural  loss  factor 

Mass  density  of  fluid  medium 

Normalized  characteristic  functions 

Natural  circular  frequency  of  vibration 
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DESCRIPTION 


Reference  7  treats  a  simply  supported,  thin,  rectangular  plate  in¬ 
serted  in  an  infinite  rigid  baffle  and  loaded  with  a  dense  fluid  on  one 

side.  The  normal  vibration  velocity  field  of  the  plate  is  expanded  in  a 

* 

series  of  in  vaeuo  normal  modes.  The  effect  of  structure- fluid  interaction 
leads  to  the  coupling  of  in  vaeuo  modes  represented  by  an  infinite  set  of 
simultaneous  linear  equations  to  be  solved  for  the  infinite  number  of 
unknown  modal  response  amplitudes.  Fluid  loading  terms  or  coefficients  in 
these  equations  are  defined  by  integrals  which  are  evaluated  approximately 
for  various  regimes  of  frequency.  Coupled  and  uncoupled  plate  modes  are 
included.  The  imaginary  part  of  the  coefficients  associated  with  these 
modes  leads  to  a  virtual  mass  which  is  added  to  the  plate  mass.  This 
causes  a  decrease  in  the  modal  resonance  frequencies. 


DERIVATION 

Assume  that  neither  the  panel  vibration  nor  the  acoustic  field 
affects  the  applied  external  force  acting  on  the  thin  panel  (Figure  3) . 
The  equation  of  motion  representing  the  normal  displacement  of  the  panel 
driven  by  a  pressure  field  is  then 


DV^y  +  m  8  —  +  m  =  p(x,t)  -  P(x,x0  =  0,t) 

P  3t  P  2  "  ~  1 


As  discussed  in  Appendix  A,  the  wave  number  spectrum  of  the  structure 
is  discrete  and  that  of  the  acoustic  field  is  continuous.  Hence  for  the  sub¬ 
merged  plate,  in  vaeuo  normal  modes  do  not  exist.  However,  the  expansion 
of  the  velocity  response  of  the  structure  in  terms  of  its  in  vaeuo  modes  is 
still  valid.  For  convenience,  we  refer  to  these  functions  as  modes  and  also 
refer  to  the  resonance  frequencies  of  these  modes.  Thus,  we  do  not  refer 
to  a  frequency  associated  with  some  natural  mode  of  vibration  but  rather  to 
a  frequency  corresponding  to  a  maximum  value  of  the  amplitude  response  of  a 
mode.  The  coupling  together  of  the  in  vaeuo  modes  by  the  structure-field 
interaction  is  a  significant  aspect  of  this  problem.  The  effective 
coupling  depends  on  both  wave  number  matching  and  resonance  frequency 
proximity  and,  therefore,  on  the  relative  magnitudes  of  the  widths  of  the 
resonance  peaks  and  the  frequency  spacing  of  the  resonances. 
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The  nodal  equation  for  the  frequency  Fourier  transform  of  panel  velocity 
is  then 


4  7  2 

[Dk  -  u~m  -  ia  a  i) _  sen  wlv  fw)  =  -iup  +  iajP 

nn  p  p  «n  '  mn  *nn  i 


nn 


(82) 


obtained  by  use  of  the  following  relationships 


where 


O 

r(x,u)  =  J  v(x,t)e”1Ut  dt 


(83) 
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sin  k  x.  sin  k  x, 
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mw  nn 

is  a  nornal  mode  of  a  siaply  supported  panel .  Here  k  =  ,  k  =  , 

n  Xj  n  x3 


A  =  £,  and  k  2  =  k  2  +  k  2.  Also 
p  1  o*  Bn  n  n 


B  =  ns  |w|  =  ns  oj  sgn  u 
and  (see  Chapter  V  of  Reference  8) 


(86) 


W*0  =  j  P(?’w)  ¥nn(^  d? 
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P 

p«n(w)  =  J  d5 

A 
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where  dx  =  dx,dx„ 

1  o 

The  boundary  condition  in  the  plane  of  the  plate  relating  the 
acoustic  pressure  and  panel  velocity  v(x,t)  =  3y/3t  is 


3P(i,x2,t) 


avfe,t) 

3t 


(B9) 


so  that 


3P(x,x2,u) 


iupov(x,u) 


(BIO) 


Using  Equation  (BIO)  together  with  the  wave  equation  for  P(x,x2>ci)) 
in  an  acoustic  medium,  the  wave  number- frequency  transform  relating  the 
acoustic  pressure  to  the  panel  velocity  (see  Equation  (DIO),  Appendix  D)  is 


P(k,x  w)  =  Z(k,u)v(k,u)6ix2|/k  2  k2 


(Bll) 


idiere  the  radiation  impedance  Z(k,u)  is  given  by 

2  v  -1/2 


Z(k,o)  =  pqco 


(-  -J 
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Here  k  =  k  and  k  =  — . 

1  - 1  o  c 

o 

After  some  rearrangement,  we  obtain  from  the  above  definitions 
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PmW=-2Z  V(a,)  Jl  sqr(k)dk  (B13) 

1  q,r=l 


=  /  R  (w)  v  (w) 
mnqrv  J  qrv  ’ 


(B14) 


q,r=l 


where  dk  =  dk,dk„  and  S  (k)  is  a  shape  function  defined  by 
1  2  mn  -  v  1 
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(B15) 


S  fk) 
mn  - 


=  f  ¥  < 

I  mn 


(x)  e1-  '  5  dx 


and  Rjjjnqj.  the  coupling  coefficient  connecting  the  (m,n)  mode  and  the  (q,r) 
mode  is  defined  by 


R  (w)  = 

"nt*r  (27T) 


1  °° 
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o>)  S  (k)  S  (k)  dk 
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If  q  -*■  m  and  r  -*■  n,  then  the  nx>dal  coupling  coefficient  becomes  the  modal 
radiation  coefficient.  The  modal  coupling  coefficients  connect  the  vi¬ 
bration  of  one  plate  mode  with  that  of  other  plate  modes  because  of  plate- 
fluid  interaction.  The  modal  radiation  coefficients,  which  can  be  obtained 
as  special  cases  of  the  modal  coupling  coefficients,  are  a  measure  of  how 
efficiently  a  particular  mode  shape  resonates  when  no  other  modes  are 
excited.  The  real  parts  of  the  coefficients  are  associated  with  a 
radiation  damping  effect  on  the  plate  response.  The  imaginary  pa-  f ;  lead  to 
a  virtual  mass  to  be  added  to  the  structural  mass  of  the  plate,  thereby 
diminishing  the  modal  resonance  frequencies.  The  coupling  coefficients  can 
therefore  be  written 


R  =  S  +  iT 
mnqr  rnnqr  mnqr 


and  the  equations  of  motion  can  be  written 


(B17) 
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iojAJ  R  v  -iup  (o>)  (B18) 
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Approximate  values  of  T  =-  am  are  now  obtained  for  (1)  the 

mnqr  mnqr 

entire  frequency  range,  (2)  the  low  frequency  range,  and  (3)  the  high  fre¬ 
quency  range.  Moreover,  for  these  frequencies,  values  of  are  ob¬ 

tained  for  various  wave  number  domains,  i.e.,  for  edge  and  coimer  modes. 
Finally,  we  observe  that  the  acoustically  slow  edge  modes  are  the  major 
contributors  to  the  virtual  mass  (see  Figure  4). 
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ENTIRE  FREQUENCY  RANGE 

The  shape  functions  used  in  Equation  (B16)  are: 


f  ik-x  2k nk  „[M)Vikl£l-1][f-l)Vik3£3-l] 

=  V  (x)e  -  £  dx  =  -2-H - (B19) 

A  (k.  -k  2)(k/-k 

p  »p'-lm''vony 


Hence  the  integrand  of  Equation  (B16j  includes  terms  of  the  form 

!2[l-(-l)m  cos  kjJlj)  for  m,q 

both  odd  or  even:  (B20) 

0  otherwise 


(An  analogous  equation  holds  for  n  and  r.)  Each  mode  is  thus  coupled  to  at 
most  only  one-quarter  of  all  the  other  modes.  The  coupling  coefficients 
are  then  written 


R  =  S  +  iT 
mnqr  mnqr  mnqr 
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(B21) 
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We  consider  only  the  imaginary  part  of  the  integral  which  is  the  mass 
loading  coupling  term  ^mnqr •  This  part  of  the  integration  is  performed 
o\'er  all  values  |kj  >  kQ,  i.e.,  the  region  containing  acoustically  slow 
modes.  Inspection  of  the  integrand  indicates  that  the  largest  coefficients 
are  those  having  either  m=q  or  n=r,  or  both.  In  Equation  (B21)  we  let 
(ko2-kj2-k32) ^  -i(kj2+k3^-ko2)  ,*^2  The  region  of  integration  for 

Equation  (B21)  is  now  divided  into  three  regions  covering  acoustically  slow 
modes,  i.e.,  exterior  to  the  acoustically  fast  region  (see  Figure  4). 
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Edge  Modes 

It  is  clear  from  Figure  4  that  the  first  of  the  integrals  given  in 
the  right  aeaber  of  Equation  (B22)  are  the  dominant  ones  for  the  two  X-type 
edge  nodes.  We  first  perform  the  k^  integration.  Let 


[l-C-l)"  cos  k.i.] 

I  (kj  = _ Li— 

"q  1  o’)  oo 
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(B23) 


A  graph  of  this  function  is  plotted  in  Figure  5.  The  function  is  such  that 
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where  Inr(k_)  as  defined  by  Equation  (B23)  is 
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Ke  observe  that  there  is  no  contribution  from  the  signularities  at 

k,  =  k  jfk  or  k.  =  k  /k  for  n/r  since  in  that  case  the  term  in  brackets  in 
.>  n  r  o  r  n 

the  numerator  of  Equation  (B27)  becomes  zero  and  I^Qc^)  °-  Hence  the 
only  singularity  of  the  integrand  is  the  root  singularity. 

We  now  treat  the  terms  1  and  (-l)n  cos  k,£_  separately  as  Cauchy 

2*  *2  2  2  2  2 
principal  values.  We  make  the  approximation  k^  -k  ^-k^  and  k^  -k^  ^-k^ 

(since  there  is  little  contribution  from  k,~k  ,  k,  —k  we  consider 

3  n  3  r 

kn,kr  >  k,  from  which  the  approximation  follows).  We  write  for  the  first 
principal  value 
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for  k  »  k  -k  ,2k  »  k  +  k  .  For  coupling  between  resonant  edge  modes 
nom’nom  1 

only,  that  is  for  modes  close  together  in  wave  number  space. 
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The  integral  associated  with  the  (-1)  cos  k^l^  term  is  shown  in 
Reference  7  to  lead  to  an  asymptotic  result  due  to  the  square  root  singu¬ 
larity.  The  magnitude  of  this  result  can  be  ignored  in  comparison  with  the 

dominant  term  given  by  Equation  (B30).  Hence  for  the  edge  modes,  since 
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Finally  when  n=v  so  that  ^n=^r>  Equation  (B27)  becomes  for  acoustically  slow 
modes 
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by  virtue  of  the  analog  of  Equation  (B24)  for  I  n(kj) . 
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the  integral  in  Equation  (B27)  becomes 
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following  a  procedure  similar  to  that  used  in  evaluating  the  integral  for 
the  edge  modes.  Hence  Equation  (B27)  yields 
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Finally  when  n=r  so  that  k^=kr,  Equation  (B27)  becomes  for  acoustically 
slew  modes 
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Moreover,  inspection  shows  that  the  mass  coupling  is  negligibly  small 

for  two  acoustically  fast  modes. 

Retaining  only  the  dominant  terms.  Reference  7  summarizes  the  values 

of  T  as  follows: 
mnqr 

For  (m,n)  an  X-type  edge  mode  and  (q,r)  an  edge  mode. 
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For  large  inertial  coupling  between  modes,  it  is  necessary  that  two 
mode  numbers  be  the  same  so  that  the  modes  vibrate  in  the  same  shape  in  one 
direction.  This  is  symbolized  by  the  Kronecker  delta  functions. 


For  (m,n)  an  X-type  edge  mode  and  (q,r)  a  corner  mode. 
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For  (m,n)  an  X-type  edge  mode  and  (q,r)  an  acoustically  fast  mode. 
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For  (m,n)  and  (q,r)  comer  modes. 
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For  all  modes  in  K  space  for  which  |k|>  kQ  the  self-inertia  term  is 
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Equations  (B43)  and  (B44)  may  be  combined  to  give  the  following  result  valid 
for  all  corner  modes 
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is  of 


That  is,  T  is  the  same  for  all  acoustically  slow  modes  and  T 

mnmn  mnqr 

the  same  form  for  all  cases  treated,  irrespective  of  the  division  by 
radiation  characteristics  into  edge  and  comer  modes.  Equation  (B18)  is 
now  written  as 
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HIGH  FREQUENCY  LIMIT 

At  high  frequencies,  kQJ^ ,  »  it,  the  radiation  and  coupling 

characteristics  of  the  modes  are  not  the  same  for  all  modes  as  in  the  low 
frequency  case.  The  important  inertia  term  is  the  self-inertia  defined  by 
Equation  (B44)  so  that 
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The  high  frequency  analysis  was  restricted  to  frequencies  below  the 
acoustic  critical  frequency,  that  is,  the  resonant  modes  considered  all 
have  wave  speeds  on  the  plate  less  than  the  acoustic  wave  speed.  Hence  the 
case  of  resonantly  excited  acoustically  fast  modes  are  not  treated.  Since, 
however,  the  acoustic  critical  frequency  for  a  1/4  in.  steel  plate  in  water 
is  about  40,000  Hz,  the  restriction  is  oi:  no  great  practical  significance. 

For  still  greater  refinements  than  the  results  presented  here,  the 
reader  is  referred  to  Reference  7;  the  lefinements,  however,  do  not  sig¬ 
nificantly  alter  the  results  presented  here. 


APPENDIX  C 

THE  LEIBOVITZ  METHOD  I 


NOTATION 
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■npq 
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■npq 
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pq 


S(Y.-k.) 

1  y 


8.. 

ij 

Yx’Yy 


P 
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Area  of  plate 

Velocity  of  sound  in  fluid  Medium 

Driving  force  applied  at  coordinates  x=o,  y=o  (also  2=0)  i.e. , 
origin  of  coordinate  systca,  used  in  Appendix  A,  on  plate  sur¬ 
face 

Generalized  force  for  coupled  mnpq  nodes 

Integral  defined  by  Equation  (C2) 

Defined  by  Equation  (F25)  of  Appendix  B 

Acoustic  wave  number  equal  to  h/c 

Modal  wave  numbers  defined  by  Equation  (A3) 

T 


Surface  wave  nuaber  equal  to 


K1-2  *  k»2 


Wave  nuaber  in  the  direction  of  the  ordinate  (defined  as  in 
notation  for  Appendix  B) 

Half  length  and  half  width  of  plate,  respectively 

Equal  to  2£  and  21  ,  respectively 
x  y 

Mode  numbers 

Added  mass  per  unit  area 

Mass  of  coupled  nodes  *npq  per  unit  area 

Modal  mass  coupling  tern 

Displacement  amplitude  of  vibration  of  a  plate  for  the  pq 
mode 


Delta  function 
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1,  y.=k. 

1  j 

0,  Y.jfk- 

1  J 


Kronecker  delta  equal  to  1  for  i=j ;  equal  to  0  otherwise 

Wave  numbers  which  are  the  coordinates  in  Fourier  transform 
space 

Mass  density  of  fluid  medium 

Natural  circular  frequency  of  vibration 
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DESCRIPTION 


It  was  indicated  in  Appendix  A  that  Feit  and  Junger  avoided  the 
determination  of  the  exact  solution  of  Equation  (A17).  They  obtained  an 
approximate  solution  by  evaluating  the  integral  of  the  equivalent 
Equation  (A14)  for  the  high  wave  number  limit  and  u  <  u^.  In  the  present 
Appendix  we  attempt  to  obtain  approximate  solutions  directly  from 
Equation  (A17)  using  with  some  modifications  the  methods  and  results  of 
Davies  given  in  Appendix  B.  Thus  with  proper  interpretation  of  the  results 
as  applicable  to  various  modal  regions,  the  work  of  Appendices  A  and  B  can 
be  interrelated. 


DERIVATION 


For  a  single  node,  the  double  integral  of  Equation  (A17)  may  be 

* 

rewritten  in  the  following  form;  note  that  for  a  single  mode  the  summation 


J  is  dropped  and  wave  numbers  k  ,k  are  considered  as  known  quantities 

pq  p  q 

corresponding  to  any  particular  set  of  prescribed  mode  numbers  pq.  Thus 
Equation  (Cl)  represents  the  contribution  of  a  single  prescribed  pq  mode 
to  F  (see  footnote  to  sentence  above  Equation  (C29) ) . 
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The  problem  is  first  discussed  in  terms  of  the  results  for  various  modal 
regions.  Subsequently,  it  is  shown  that  for  the  even  modes,  the  general  in¬ 
tegral  expressions  for  the  virtual  mass  (or  mass  reactance)  obtained  from 
Equations  (A17)  and  (B21)  are  identical. 
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since  all  tents  in  the  integrand  are  even  functions  in  Y  and  Y  . 
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Ne  evaluate  I  for  the  following  cases: 

Case  I:  n/p,  n/q 

Ke  observe  that  for  n/p  there  is  no  contribution  to  the  integral  froa 
the  singularities  at  Y  =k ft  or  Y  =k  /k  because  in  that  case  the  ten  in 

X  m  p  X  p  ■ 


the  numerator  of  the  first  bracket  of  Equation  (C2)  bccoaes  zero  (see 
Equation  (A3)).  Evaluation  of  the  resultant  indctendnate  quantity,  using 
L"  Hospitals  rule,  yields  c  null  result,  i.e. , 


(  1  .  1  *  COS  2Vx\  *x  .  Sin  2knJ 

"k^U2-^2  k^-Y^  /  k  2-ka2  *■ 


lin 

Y  -L 

x  ■ 

k Jk 


P  ■ 

- - ?in  _  ^0,1,2— 


,  2  .  2 
k  -k 

P  ■ 


=  0  for  all  k 


A  similar  result  is  obtained  for  all  k  for  n/q.  Thus  we  conclude  that  the 

n  n 

singularities  at  Y  =k  /k  or  Y  =k  /k  and  Y  =k  /k  or  Y  ft  ft  nake  no  con- 
6  x  ■  p  x  p  m  y  n  q  y  q  n 

tribution  to  the  integral. 

Moreover  when  Yx=km/kp,  then  in  the  first  bracket  of  Equation  (C2), 

cos  2Y  £  =  cos  2k  £  =  cos  (2m+l)n  =-  1  for  all  m.  Since  cos  (2m+l)Tr  =- 
x  x  m  x 

(-1)  cos  nnr,  then 
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(1  ♦  cos  2Y  £  )  [1-C-l)"  cos  2k 1  ] 

- 5-5—  -»• - —50  for  Y  =  k  A  only  (C3) 

9  „  9  7  „  ?  „  7  J  „  7  7  *  *  P 

<VY>  v - V: »  <Yk.>  ^  - v} 

Si*ilar  relationships  obtain  for  Y  =k  A  ,  Y  =k  /k  ,  and  Y  =k  /k  .  For 

Y^/k^ ,  the  for*  of  the  left  *eaber  of  Equation  (C5)  differs  fro*  the  fora 

of  1  (k.)  C=I  (Y  )  here)  given  by  Equation  (B25)  of  Appendix  B  and  plotted 

■CJ  1  ®P 

in  Figure  5,  only  for  even  values  of  a  by  the  sign  in  the  numerator  pre¬ 
ceding  cos  Thus,  by  analogy  with  the  work  in  Appendix  B  and  fro*  an 

inspection  of  Figure  5,  it  is  apparent  that  the  integrated  contribution  is 
very  saall  over  the  range  (Ky^3*- 

The  foregoing  is  compatible  with  the  stateaent  aade  in  Appendix  B  to 
the  effect  that  the  largest  contributors  to  the  integral  occur  when  either/ 
or  a=p  and  n=q.  Appendixes  A  and  B  also  indicate  that  the  doninant  contri¬ 
bution  occurs  for  *=p  and  n=q,  i.e. ,  the  important  inertia  tern  is  the  self¬ 
inertia  ter*.  These  cases  are  discussed  next. 

Case  II:  ■= p,  n/q  or  n=q,  */p 

For  *=p,  n/q  we  use  the  Methods  of  Appendix  B.  For  this  case  we 
* 

have,  using  Equation  (B24) 

(1  +  cos  2Y  £  )  ri-C-lF  cos  2k JM  *21 

- —  = - - *  o(Y  -k  )  for 


(k  2-Y  2)(k  2-Y  2)  (Y  2-k  2)CY  2-k  2)  4ki 

^  ■  XJVP  X  '  v  X  ■■''•X  p 


Equation  (C2)  then  becomes 


Y  =k  =k 
x  ■  p 


n2 £  r  f  1  +  cos  2Y  £  “1  dY 

=  |  - ^ - y — 

4km  oLCk  2-V  2)0c  2-Y  2)  fY  2-(k2-k  2))1/2 
n  y  q  y  -l  y  n 


In  Equation  (B24),  let  l=2Z  and  k.  -*■  Y  ;  for  later  use,  we  note  that 
Z-=2l  .  l  x  i  x 

3  y 
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For  the  Z-izrpe  edge  nodes 


s2l  r 

4  f 

AV  A  J 


[1  +  cos  2Y  £  1 

y  y 


AV  j  2?  9  9  9  9  9  1/9 

*  k  [(k  -Y  j(k  ~-Y  z)]  [Y  -(k~-k  “)]*/z 

neglecting  as  in  Appendix  B  the  contribution  of  acoustically  fast  nodes 
occurring  in  the  region  0  <  Y^  <  kQ.  The  integral  in  Equation  (C6)  is 
siailar  in  fora  to  that  of  Equation  (B27).  Using  the  Cauchy  principal 
value  for  the  terns  1  and  cos  2T  ly  in  Equation.  (C6) ,  we  find  as  before 
that  the  integrated  contribution  of  the  cos  2Y  £  tern  can  be  ignored. 

y  y 

Siailar  results  are  obtained  for  n=q,  njfp.  Hence  the  results  obtained  for 

edge  cades  in  Appendix  B  can  be  used  directly. 

2  2 

Sinilarly  for  the  earner  nodes  when  k  «  k  ,  Equation  iC4)  is 

s 

written  as 


r ri  +  cos  2Y  £.  ~l 

Jo[ck  2-y  2)<k  2-y  2] 
L  n  y  q  y 


[Y  2.k  2],/2 
1  y  m  1 


The  integral  is  siailar  in  fora  to  the  left  aeaber  of  Equation  (B35). 

Hence  using  the  Cauchy  principal  value  for  the  tera  1  and  cos  2Yy£y  in 
Equation  (C7),  we  find  that,  as  before,  the  integrated  contribution  of  t£e 
cos  2Y  £  tera  can  be  ignored.  Siailar  results  are  obtained  for  n=q,  a^p. 
Hence  the  results  obtained  for  corner  codes  in  Appendix  B  can  be  used 
directly. 

Using  Equations  (B40)-(B43)  we  now  find  values  of  I,  F  ,  and  m 

Dul  DUl 

for  various  cases  of  nodal  coupling.  Substituting  Equation  (C2)  in 
Equation  (Cl),  we  have 

4CcVk  k  k  k 

m  m  p  q 

F  =  F(0,0)  +  -  K  •  I  (C8) 

mn  v  ’  J  _  pq 


Inspection  of  the  first  two  members  of  Equation  (C4)  indicates  that 
for  Yx=km=kp ,  k^k^,  the  integrand  of  Equation  (Cl)  and  of  Equation  (B21) 
times  -i  are  identical  if,  when  using  the  Cauchy  principal  values,  we 
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££*>'**  •  ft" Si-1 ^  *' 


neglect  the  cos  2TML  and  cos  k_£_  tens.  A  similar  identity  exists  for 
Yy=kn=k^,  k^/kp.  Hence  substituting  Equation  (C2)  in  Equation  (B21),  we 
get 


1  *  MpeTTTTT  =  *P  *  (2V  =  41  A 

*  n  p  q  v  3  3 

(2*)2  A 


where  we  have  considered  only  the  inaginary  part  of  the  integral  in 
Equation  (B21). 

Substitution  of  Equation  (C9)  in  Equation  (C8)  yields 


ck(-l)B+n+P^  A 


F«n  =  F(0»0)  ' 


If  T 

pq  »npq 


(CIO) 


where  T  is  given  by  Equations  (B40)-(B43).  (Note:  Davies  svnbol  q  •+•  p 
and  r  -*■  q  here;  ■,  p  stay  the  sane,  and  all  odd  node  niaabers  used  with 
respect  to  Davies  origin  for  the  plate  represent  even  nodes  with  respect 
to  the  Feit-Junger  origin  for  the  plate;  see  Appendix  D.) 

For  (*,n)  an  X-type  edge  node  and  (p,q)  an  edge  node. 


F  =  F(0,0)  - 
nm 


„  k2l  21 

’AP  V,  k 


k  k 

4oc  m  P 

k2£  - .1 

A  y  .  2t  2  W 

rnn  pn 


W  6  -(-l)B+q 

pq  ®p 


a  1 

i  pnj 


=  F(0,0)  +  C-l)n+q£  in 

TT  L  X 


K  5 
pq  <np 


(-l)®+q£,  k  k 
y  m  p 


(Cll) 


k  2k  2Fq  Pn 
mn  pn 


(C12) 


Hence 
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miU&i/AtV  •>  i *£ Hili  *  'tMixb  (lit  s  iff  a 


J 


For  (m,n)  and  (p,q)  corner  nodes. 


mn 


ckA  r-(-l)n+q4pck2£  k  k  4pck2£  k  k 

F  =  F(0,0)  -  _P  _ - If  « P - Y-Z- 

4  [  ,  .  t  2,.  2  «  ”P  if  A  k„  V  2 


.  .  2.  2 
A  k  k 
p  mn  aq 


W  6 

pq  nq 


l 


■jr  A  k 

pm  pn 


(C20) 


2 

;! 

!? 

ii 


=  F(0,0)  + 


2  2  r 

2pc  k 

if 


C-l)n+q£  k  k 
x  n  q 


.  2,  2 
k  k 
mn  aq 


K  S 

pq  nrp 


C-l)n+P£  k  k 

y»p 


w  s 


2k  2  W  nq 


mn  pn 


■] 


(C21) 


Hence 


a  =4 
mnpq  if  I 


C-l)n+qk  k 


Ik  h 

y  an  aq 


n  q 

- 6  + 

2  “P 


C-l)m4pk  k 


Ik  !t 
x  mn  pn 


k  “I 
m  p 

- 6 

2  nq 


(C22) 


Equations  (C13),  (C16) ,  (C19),  (C22)  are  identical  to  corresponding  results 
obtained  from  Appendix  B  and  presented  as  results  in  Table  1  as  Items  2-5 
for  the  Davies  method.  Having  obtained  the  results  for  the  added  mass  of 
the  coupled  modes,  i.e.,  the  coupled  inertia  terms,  we  now  consider  the 
added  mass  of  the  uncoupled  modes,  i.e.,  the  self-inertia  terms. 

Case  III:  m=p,  n=q 

At  m=p,  n=q  there  occurs  a  dominant  contribution  to  the  integral  in 
Equation  (Cl)  from  the  singularities  at  Yx=km=kp  and  Y^=kn=k^ .  For  this 
case 


(1  +  cos  2Y  i  )  [l-(-l)  cos  2k  l  ] 


x  x 


2  2  2 

(k  -Y/r 

m  x 


7r  •  2£ 


(k  2-Y  2)2 
^  m  x  J 


=  I  (Y  )  = 
mm  x 


6(Y  -k  )  (C23) 

v  x  a 


4k 


m 


% 

1 


(1  +  cos  2Y  £  )  [l-(-l)n  cos  2k  £  ]  n  •  2  £ 

: - ^  =  ^n^y)  «  — — ^  5CYy-kn)  (C24) 


(k  2-Y  V 

v  n  y  ' 


(k  2~y  2)2 

v  n  y  1 


4k 


Hence  Equation  (C2)  becomes 


I  = 


TT2  £  £ 

x  y 

2  2 
4k  ^k  Z 
m  n 


«  6(Y  -k  )  6(Y  -k  )  dY  dY 
■-  v  x  nr  v  y  n  xy 


JJ 

0  [Yx‘-Ck‘-Yy‘)] 


2  „  2,, 1/2 


(C25) 


2 

tT  £  £ 
x  y 


4k  2k  2[k  2-k2]1/2 

m  n  1  mn  1 


(C26) 


Substituting  Equation  (C26)  in  Equation  (Cl),  we  get 


F  =  F(0,0)  + 
mn  v  *  ’ 


pc2k2  £  £  W 
x  y  mn 

rv  2  .2,1/2 

*-Kmn  ~  ■* 


(C27) 


Equation  (C27)  is  identical  with  Equation  (A27) .  Hence  the  results 

given  by  Equations  (A28)-(A34)  apply  here.  We  note  that  since  the  mass 

loading  is  associated  with  the  imaginary  part  of  the  integral  for  which 
2  2  2 

Y^  +  Yy  >  k  (i.e.,  acoustically  slow  modes)  then  from  Equation  (A28) 


m  = 
mn 


(k  2-k2)1/2 
v  mn 


2>k2 

.  mn 

mn 


(C28) 


The  approximation  value  of  m^^  given  by  Equation  (C28)  agrees  with 

the  results  obtained  in  Equation  (B44)  since  T  =-  com  =-  kcm  .  The 

n  mnmn  mn  mn 

result  is  therefore  identical  to  the  corresponding  results  presented  in 
Table  1  as  Item  1  for  the  Davies  method. 
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PROOF  OF  IDENTITY  OF  FEIT-JUNGER  AND 
DAVIES  GENERAL  INTEGRAL  EXPRESSIONS  FOR 
THE  VIRTUAL  MASS  (FOR  EVEN  MODES) 

* 

The  double  integral  of  Equation  (A17)  may  be  rewritten  in  the  form 


F  =  F(0,0)  + 
mn 


i4pc2k2k  k  (-l)m+n 
m  nv 


A  W  k  k  (- 

pq  p  q 


00 


1  +  cos  2Y„£. 


2  2  2  2 
-Y  *)  (k  -Y  * 
m  x  J  K  p  x 


JL  “I 
X  X 

2  2 

-y  n  J 


[1  +  cos  2Y£  1 

(k  2-Y  2) (k  2-Y  2) J  [k 

vn  y  q  y 


dY  dY 
x  y 

2_y  2_y  2,1/2 
x  y  ■* 


(C29) 


In  Appendix  A  (Feit-Junger) ,  the  origin  is  taken  at  the  center  of 
the  plate  whereas  in  Appendix  B  (Davies)  the  origin  is  taken  at  a  comer 
of  the  plate.  As  shown  in  Appendix  D  and  Figure  11  (see  Appendix  D),  the 
modes  numbered  m  =  1,  3***  odd  with  respect  to  Davies  origin  represent  the 
even  modes  with  respect  to  Feit-Junger  origin.  By  the  Feit-Junger  stipu¬ 
lation,  Equation  (Al),  these  are  the  only  modes  to  be  considered;  similar 
relations  hold  for  n,p,q.  Hence  [-l]m+n  =  [-1]P+^  -*■  l  and 


1  +  cos  2Y  Z 


2  2 

(k  -Y  ^)  (k 
.  m  x  p 


2Y  Z  “I 
x  x 

2  2 
)  (k  -Y  ) 
p  x  J 


Feit-Junger 


”  l-(-l)  cos 

.(ko2-ki2)(k  2-ki2)J 


Davies 


"  i  +  cos  2Y  z  ~i  r  i-(-d  cos 

(kq2-Yy2)ieit-Junger  L'kn2-k52>  «r2-k32’  J 


Davies 


Here,  in  contrast  to  Equation  (Cl),  the  contributions  of  all  (p,q)  modes 


to  F  are  included, 
mn 
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-  '‘50*=“'- 


where  the  Feit-Junger  notation  m,n,p,q  m,n,q,r,  respectively,  and  2lx , 

2£y  -*■  S-j ,  &3,  respectively.  Also  Y^/Y  -»•  kj,  k3  =  components  of  k,  k  ->  kQj 

pC  *  PoCo*  Thus  the  imaginary  part  of  Equation  (C29)  rewritten  in  Davies 
notation  is:  co  » 


^for  the  imaginary  part,  ll  *  ii  md  we  let  (ko2-kj2-k32)^2 

o  |kj>k 


.  ,,  2  .  2  .  2,1/2^ 

-icki  +k3  -ko  )  ; 


(Imag)  Fmn  =  F(0,0)  - 


4p  c  2k  2k  k 

r  a  A  rr> 


i  c  k  k  k  __  rr  r  .  ,  ,,m  ,  „ 

O  o  O  m  n  y^w  kk  l-(-l)  cos  k  ^ 

q,r  |k|>ko  Lvki  "V }J 


'  l-( 


-n  cos  k3£3 


j  [kx2 


(C30) 


2  ,  2w,  2  .  2,  r.  2,2,  2,1/2 

n  -k3  3  (kr  “k3  U  -k0  ] 


In  Davies  notation  the  expression  for  the  virtual  mass  in 
Equation  (A18)  may  be  written 

(F  )  =  (Imag)  F  =  •  ■=£  Y*  W  co2m  =-  Y*  W  w  T 

mn  °  mn  2  2  qr  mnqr  4  ^  ^  qr  n 

Virtual  q,r  q,r 


(C31) 


where  A  =  JLJL,  T  =-  com  =-  kcm  and  where  now  T  includes 
p  13  mnqr  mnqr  mnqr  mnqr 

T  =-  wm  (=  -com  in  Equation  (A18))  and  q,r  includes  m,n  in  the 
mnmn  mnmn  mn  n  n 

summation. 

Comparing  Equation  (C31)  and  (excluding  the  driving  force  F(0,0)) 
Equation  (C30) ,  we  get 


S  'V 


64p  c  k  k  k  k  k 
o  o  m  n  q  r  o 

(2ir)2  A 


l-(-l)m  cos  Vi-iri-C-1)"  cos  k3£3 


l*l*o  K'V3  (k  2-ki2lJbkn2-k323  ^r'V3 


ri  2  .  2  ,  2,1/2 

[k2  +k3  -kQ  ] 


(C32) 


In  Equation  (B21)  let  ( kQ 2  -  kj2  -  k3V/2  *  +  -  k,2)-^ 

We  then  see  that  Equation  (C32)  is  identical  to  the  imaginary  part  of 
Equation  (B21)  which  has  the  same  integral  limits.  Hence  for  m,n,p,q  equal 
to  odd  numbers,  representing  Feit -Jungers  even  modes,  all  solutions  of 
Equation  (B21)  for  obtained  by  Davies  methods,  applicable  to  various 

modal  regions,  are  relevant  as  solutions  to  the  Feit-Junger  Equation  (A17). 
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APPENDIX  D 
LEIBOWITZ  METHOD  II 


NOTATION 


A 

c 

cfi 

c 

p 

E 

e 


^rad 

m 

_rad 

F 

m 


K(ks) 

k  ,k 
m’  n 


mn 


n 


k  ,k 
X1  x2 

kl,k2,lc3 

kl,k2,k3 

^l’^2 


Area  of  plate 

Velocity  of  sound  in  fluid  medium 
Compressional  wave  velocity  of  the  plate  equal  to 

Phase  velocity 


ps-r 

[ps(l-v2)J 


Young’s  modulus 

Equal  to  2.718;  base  for  natural  or  Naperian  system  of 
logarithms 

Radiation  force  for  the  mtk  mode 


th 


2  2 

^kg  -ko  in  Equation  (D38) 


Magnitude  of  the  radiation  force  of  the  m  mode;  equal  to 
-rad.  -iwt 
m  e 

Defined  by  equations  below  Equation  (T>36) 

Equal  to  ifT 
Equal  to 

Modal  wave  numbers  equal  to  and  respectively 

-  _ 

Nave  number  equal  to  ykm  +kfi 

Projection  of  kQ  on  a  normal  to  a  plane  lying  along 

Acoustic  wave  number  equal  to  w/c 

Wave  vector  equal  to  k  r 
n  o  o 

2  2  2 

Arbitrary  wave  number  equal  to  27i/X  =  w/c  ;  k  =  k  +  k- 

p  P  S  I  ^ 

Projection  of  kQ  on  a  plane 

Wave  numbers  in  x^-  and  x2-directions,  respectively 
Wave  numbers  along  the  x^-,  x^- ,  x3-directions,  respectively 
Vector  wave  numbers  in  the  x^,  x2-,  x^-directions ,  respectively 
Length  and  width  of  plate,  respectively 
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! 


L  ,L 
x  y 


M 


loading 


m 

mn 


m,n 


t 

VW 

V(ks) 


X1 ,X2,X3 


Half  length  and  half  width  of  plate,  respectively 
Total  added  (or  virtual)  mass 
Added  modal  mass  per  unit  area 


Mode  numbers 


Amplitude  of  acoustic  pressure  p 

Acoustic  pressure  in  fluid  medium,  i.e.,  in  half-space  x^  >  0 
Radiation  pressure  for  the  m*11  mode 

Sound  flux  or  source  strength  or  volume  velocity  of  j**1  source 
equal  to  the  product  of  the  velocity  of  the  source  and  the  sur¬ 
face  area  of  the  source 


Distance  from  origin  to  a  field  or  observation  point 

Distance  from  j**1  source  to  the  field  or  observation  point; 
equal  to  R-r^r^  which  is  the  projection  of  R^  on  R  (here 

rQ*r.  is  a  dot  producted  equal  to  the  projected  distance 

difference  R-R.) 

3 

Rest  position  of  the  plate 


Magnitude  of  wave  direction  vector,  and  wave  direction  vector 
respectively;  rQ  is  unit  vector  from  origin  to  field  or  ob¬ 
servation  point 

Magnitude  and  vector  for  distance  from  origin  to  jth  source 
lying  in  the  plane 

Time 


Instantaneous  modal  velocity  .of  a  point  on  the  plane 

Velocity  of  plate  at  its  rest  position  r;  equal  to  plate 
velocity  at  =  0 

Distribution  of  traveling  plane  waves  in  an  infinite  plane,  i.e., 
Fourier  distribution  of  the  velocity  V(r,t)  in  the  region  of  the 
plate 

Complex  amplitude  of  V(r,t);  modal  velocity 


Root  mean  square  of  |V(r)|;  see  Equation  (D22);  equal  to  |vj 

Abscissa  with  origin  at  the  corner  and  midpoint  of  the  plate, 
respectively 

Rectangular  coordinates,  x^^  lie  along  the  length  and  width  of 
the  plate  of  the  plate,  respectively;  is  normal  to  the  plate 
and  the  origin  lies  at  one  corner  of  the  plate 
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,rad 


r  (k  ) 

mv  s' 


vr> 


ll) 

* 


<>- 

r 


Radiation  impedance 

Modal  coupling  parameter  or  directivity 

Equal  to  (32/3x2)  +  (32/3y2)  +  (32/3z2) 
Wavelength  of  plate  equal  to  2ir /ks  =  2t!c^/co 

Poisson's  ratio 
Mass  density  of  fluid 

Mass  density  of  plate 

Normal  mode  function 

Natural  circular  frequency  of  vibration 
Denotes  complex  conjugate 
Denotes  magnitude 
Denotes  average  value  over  r 


DERIVATION 

Consider  a  siuply  supported  plate  in  an  infinite  plane  baffle 
immersed  in  a  fluid.  Plane  flexural  waves  form  a  wave  field  on  this  plate. 
We  treat  therefore  a  two-dimensional  problem,  in  which  radiation  in  the 

half-space  x_  >  o  is  of  interest.  The  instantaneous  modal  velocity  of  a 

^  -  8 
point  on  the  plate,  whose  rest  position  is  r  is  given  by  (see  Figure  6) 


(Dl) 


We  expand  V(r,t)  into  a  distribution  V(ks)  of  traveling  plane  waves 
in  an  infinite  plane  each  of  the  form9  (e  ia)t  is  tacitly  implied) 

V(r)  =  V(ks)el£s‘r  (D2) 

For  a  plate  vibrating  with  arbitrary  wave  number  kg  =  2tt/A^  =  u/c^ 
and  at  frequency  »,  it  seems  logical  to  make  the  following  formulation  for 
the  sound  pressure  in  the  half-space  x3  >  0  (see  Figure  6) 


„  i(k  •  r  +  k_x_) 
p(r,x_)  =  Pe  s  33' 


(D3) 


We  require  that  the  sound  pressure  represent  a  solution  of  the 

9 

Helmholtz  wave  equation 
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Figure  6a  -  Plan  View 


Figure  6b  -  Edge  View  Normal  to  Ks;  Vector  Relations  at  r 

Figure  6  -  Plan  View  and  Edge  View  Normal  to  K  of  Plate  in 
Infinite  Plane  Baffle  s 


1.  '^(r)  and  V(r,t)  arc  normal  to  plane;  ♦  V  is  out  of  plane,  -  V  is  into  plane 

2.  Xj ,  x^,  kj,  ^2*  ks  *n  planc,  x_,  arc  normal  to  the  plane 


V2?  +  k  2p  =  0 


(B4) 


Substituting  Equation  (D3)  in  Equation  (D4) ,  we  get 


.  2  ,2  .2 

k  +  k_  =  k 
s  3  o 


(D5) 


Me  also  require  that  the  noraal  exponent  of  the  sound  velocity  cal¬ 
culated  froa  Equation  (D3)  coincide  at  the  boundary  surface  x_=0  vith  the 
plate  velocity,  i.e., 

•  (D6) 

3 


Froa  Equations  (06)  and  (D2)  ve  obtain 


P  = 


_  -ik  •  r 
upoV(r)  e  s 


(D7) 


p  ck  V(k  ) 
o  o  s 


(D8) 


PocV(ks) 


4 


CD9) 


Substituting  for  P  from  Equation  (D9)  and  for  k^  from  Equation  (D5) 
into  Equation  (D3)  and  restoring  the  sinusoidal  temporal  variation  yields 
the  following  expression  for  the  sound  pressure  radiated  into  the  half¬ 
space  in  front  of  the  plate: 


,  pc  V(k  )  ilk  *r  + 

prad  =.  0  0  5  e  V  s 

one  wave  3’ 


(v?  *  |A?V  v) 


(»10) 
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Froa  Equation  (D2)  (with  tiae  factor  included),  the  sun  or  integral 

of  traveling  plane  waves  in  in  infinite  plane  with  amplitude  distribution 

V(ks)  provides  the  velocity  V(r,t)  in  the  region  of  the  plate  and  zero 

outside  the  plate.  The  distribution  V(k  )  *ust  therefore  satisfy  the 

condition  (note  dk  =  dk.dk,) 
s  1  o 


CO 

;,t)  =  f  V(ks)  J 


“^dk 


(Dll) 


Obviously  V(k^)  is  the  Fourier  distribution  of  the  function  V(r,t)  and  we 

know,  therefore,  that  such  a  distribution  exists  and  that  it  can  be  calcu- 

*  _ 

lated  fro*  the  equation  (note  dr  =  dx^dx,) 


V(ks) 


=  -4  f  V(f,t)  ,-rts  ■ i 

a *r  J. 


■“’af 


(D12) 


The  corresponding  total  sound  pressure  obtained  fro*  the  super¬ 
position  of  all  traveling  waves  of  the  fora  represented  by  Equation  (DIO) 


rad,- 
P  (r 


n  f  V<**)  ifl 

.x„t)  =  o  c  I  - e  \ 

00  J  /~~2  v 


k  •  r  +  yin  2-k  2  x_-tot)ji; 
k  s  f  o  3  3  /  dks  (D13) 


Ke  are  interested  in  the  pressure  on  the  surface  of  the  plane. 
Hence  setting  x,=0 


00 

Tad,-  .  I 

p  (r,x3=o,t)  =PocoJ 


V(V  ci(k.  •  r  -ait)dks 


(D14) 


No  contribution  to  V(k  )  is  rendered  by  the  integral  for  the  region  ex¬ 
ternal  to  the  plate  area. 


I 


Now  substituting  Equations  (Di)  in  Equation  (D12),  we  get 

V(k 


)  -  f  V  ?  (fje^s  *  r 

(27!)2  J  *" 


dr 


^  f  V 
p»r  J,„  " 


—  -ik  •  r  — 
(r)e  s  dr 


(D15) 


(D16) 


lfe  now  prove  that  the  integral  is  related  to  a  quantity  [F^* (k^jj/2 

where  T  (k  )  is  called  the  nodal  coupling  parameter  or  directivity.  The 

■  ^  _  2 

functional  forns  and  average  values  of  |i‘B(ks)|  described  in  some  detail 

in  Reference  8  are  useful  in  making  approximate  computations. 

PROOF 

The  pressure  radiated  into  the  half-space  x„  >  0  from  an  array  of 

9  * 

sources  on  an  infinite  rigid  plane  is  (see  Figure  7) 


m  **o  V  Qj  -ik  P. 
p(R)  =  ~xr  2-  —  -  - 


2tt  —  e  °  J 
j  R. 

3 


(D17) 


Projecting  R^  on  R,  we  obtain  R.«fR-r0*?.  where  ro*r-,  the  projected  dis¬ 
tance  difference,  is  2  dot  product  (see  Figure  8).  Since  1/R_.  = 


r  *r.\  R  \  R  / 


1/R  then 


iwp 


°  e 


ik  R  V  n  ik  *r. 
o  Z-  Qj  e  0  3 


(D18) 


where  k  =  k  r  is  a  wave  vector  =  (wave  number  k  )  •  (wave  direction  r  ) . 
o  o  o  o  o 

If  kg  =  projection  of  kQ  on  the  plane  and  k^  =  projection  of  kQ  on  the 

normal  to  the  plane,  i.e.,  k  lies  along  x_,  then  k  *r.  -  k  *r.+k  *r.=k  »r. 

v  ’  ’  n  &  V  oj  njsjsj 

since  ^n*?j  =  0  (see  Figure  6).  Hence  Equation  (D18)  becomes 


P(R>r0)  = 


xwp 


2ttR 


o  -ik  R 
e 


(D19) 


1 

i 


1 


P  (FIELD  OR  OBSERVATION 
POMT  OR  RECEIVER) 


i 

t 


i 
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Figure  7  -  Array  of  Point  Sources  in  an  Infinite  Rigid  Plane 


Figure  8  -  Coordinate  System  Showing  Vector  Relationships 


Notes: 

1.  Xj,  x^,  r.  lie  in  plane  of  plate;  x^  is  normal  to  plane, 

rQ  and  kQ  are  unit  and  acoustic  spaccvectors ,  respectively, 

ro>  icQ  and  the  distances  (not  vectors)  R  and  from  the 

origin  and  source  to  the  fluid  point,  respectively,  are  at 
arbitrary  angles  to  the  plane  depending  upon  the  position 
of  the  field  point  in  space. 

2.  R.  =  R  (1  -  r  *  r./R  +  - ) 

)  o  J 

ft:  R  -  r  •  r.  for  R,  R.  »  Ir.l  or  for  |r.|  sin  <f>  small, 

°  }  3  3*  J 

i.e.,  applicable  to  plate  area  A;  which  encloses  all  sources,  for 

which  Vl*  +  t22  «  R,  R.. 


•  'J 


2 


£ 


f 


5 


* 

I 

% 


£ 

Z 


■V 


% 

i 

£ 

fj 

I 

I 

■> 

f 

V 5 
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v  |v  1  r  *(k  ) 

V(k  )  =  '  2~ 

M  v 

Now  substituting  Equation  (D25)  in  Equation  (D14),  we  get 


-p  c  V  |V  J 

rn  m  I  m  * 


(D25) 


Pi.  au  f  .  %  u  ui  w  I  ui 

(r,t)  - - =■ -  I - - 

”  Sir2  V  *  J  / TJ 


|v  I  r  r  *  ei(Vr 

m1  m‘  |  m _ .. 

v  *  J  /  7i r 


(D26) 


■  v 

The  radiation  force  for  the  mth  mode  is8  ^noting  that  p*ad  =  P*ad  e 

cd  ■  cd  ■  U c  -  j »r  c 


E - =  F —  e 

ra  m 


rad  . 

p  ip  (r)dr 
rm  Ym 


Substituting  Equation  (D26)  in  Equation  (D27), 
,  p  cV  |V  I  °C  T  *dk 

to  H  n  m 1  Tn  *  I  m  cr 


.  p  cv  v  r  r  *dk  r  -c  - 

Frad  =  °,  “  °  _5_i_  I  V  (f)elks'r, 

■  s»2v  *  J  m  J  ■ 


(D28) 


Using  Equation  (D23) 


8it  \  V  * 
\  m 


i  /p  cv  jv  |\  /  |v  j\  r  r  *r 

*2  V*  '2Vm  J  f— 


— dk 


(D29) 


and  since  V  *V  =  jv  Z| ,  T  *r  =  Ir  I  .  then  we  find  as  a  basic  working  ex- 
mm  m'  mm'm' 

pression 


crad  „  \v  1 2 

zrad=_js_=V  U - d£ 

m  V  2  ]  / - r  s 

m  167T  -»  /  k  2 


(D30) 


We  write  this  equation  as  the  sum  of  real  and  imaginary  terms: 


k0  r_k 

r  |T  | 2dk  ip  c  r°  |T  |2dk 
I'm1  s  .  Ko  |  1  m1  s 


r  jv  dks ,  ipoc  r 
J  I  7~2  16tt2  J 


|r  |2dk 
1  m'  s 


-k  /,  s 
o  ,/l - = 


(D51) 


The  first  term  in  the  right  member  represents  the  radiation  damping  whereas 
the  second  and  third  terms  represent  the  fluid  loading  or  added  mass  terms 
=  i^ioading*  We  consider  the  fluid  loading  problem  only.  (Note  that 
each  integral  in  Equation  (D31)  represents  a  double  integral  since  dks  = 
dkjdk,^.)  Since  the  integrands  are  even  functions,  the  terms  within  the 
brackets  may  be  combined  so  that  the  fluid  loading  term  may  be  written: 


ioading  -  47T2o 


|T  (k  )|2 
1  m  s' 


(D32) 


00 

=  j_!s  f 
4*2  ko  J 


(D33) 


The  product  of  the  left  and  right  members,  respectively,  of 
Equations  (D23)  and  (D24)  yields 


|rj  =4  J  t  ( 

1  m  m 

A 


J  Y Weiks'r  dr  J  Y  (?)e'iks*3 


(034) 


The  modal  shape  functions  for  the  normal  modes  of  a  simply  supported 
* 

panel  are: 


miTXj  mrx_ 

T  (r)  -  2  sin  — g — ■  ■  sin  — 5 —  =  2  sin  k  x.  sin  k  x0 
mv  y  ml  n2 


(D35) 


Hence 


|rj2 .  u  j-e^i 


sin  k  x,dx. 
m  1  1 


j*kx  X2  sin  k  x.dx  • 
2  n  2  2 


^1  ^2 

f  e  ikx,  X1  sin  k  x,dx,  f  e  *kx„  X2  sin  k  xdx  (D36) 

I  1  mill  2  n  2  2 


^1  ^2 

Defining  I±  =j-  J  e1^  X1  sin  k^dXj,  I2  =  ~  f  elkx2  X2  sin  knx2dx2 


t  J 


I1*!2  =  16a2  IlI2Il*I2*  =  16A2|I1|2|I2|2 


(D37) 


Substituting  Equation  (D37)  in  Equation  (D33) ,  we  obtain 


loading 


4p/  f  !  Ij.12!  I2|2  - 

-7-15^* 


(D38) 


i  1/k  2  -  k  2 

k0  y s  ° 


where  dk  = 


d^dk.,  k  2=k  2  +  k  2 
1  2  s  1  2 


* 

Later  in  this  Appendix  it  is  shown  that  this  normal  mode  representation 
includes  the  representation  used  in  Appendix  A  as  a  special  case;  see 
Equation  (Al). 
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'/Kk 


We  consider  kg  »  kQ  and  for  this  region  the  quantity  1 
=  K(kg)  is  a  smooth  well-behaved  function.  We  can  then  use  the  analytica. 
approximation  (presented  in  Appendix  II  of  Reference  8)  for  the  weighted 


integral  given  in  Equation  (D38) .  For  the  range  of  k  including  k  and 

*  s  ro 

k^  (see  Equation  All. 10  of  Reference  8),  the  added  mass  per  unit  area  for 
radiation  into  the  half  space  x^  >  0  is  (see  Figures  9  and  10) : 


or 


m 

mn 


M 


loading 


A 


(D39) 


m  = 
mn 


if  k  »  k 


mn 


-  k 


mn 


mn 


(D40) 


The  results  given  by  Equation  (D40)  are  in  agreement  with  results 
presented  by  Equations  (A27a) ,  (B47) ,  and  (C28) . 

Finally,  it  is  shown  that  the  normal  mode  representation  for  a 
simply  supported  plate.  Equation  (D35),  is  more  general  than  the  correspond¬ 
ing  representation  used  by  Feit-Junger,  Equation  (Al).  To  see  that 
Equation  (D35)  includes  Equation  (Al)  as  a  special  case,  consider  the  sin 
k  x  factor  in  Equation  (D35);  similar  results  are  obtainable  for  the  sin 
kny  factor  (let  k^  kn,  x  y).  We  translate  the  origin  of  the  abscissas 
from  0  to  O'  (the  midpoint  along  the  plate  length)  as  shown  in  Figure  11, 
by  use  of  the  equation: 

x'  =  x  -  -j  (D41) 


k.>k 
1  m 


f 


K(kjn)|l|2dk1 


-k  <k.<k 
m  1  m 


V 


ir 


K(km)  and  similarly  for  the  integral  for 
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Figure  9  -  Sketch  of  |l|  for  Large  Mode  Numbers  m  (Applicable  to  small 
mode  lumbers  if  lower  limit  of  integration  of  Equation  AII-10  of 
Reference  8  is  set  to  zero.  See  Figure  10  and  page  216  of  Reference  8). 
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The  new  boundaries  are: 


x=0,  x*  =-  4-  =-  L 
2  x 

x=£,  x'  =+A=+  l 
2  x 

where  the  notation  *  for  the  new  plate  boundaries  has  been  introduced 
for  compatibility  with  the  Feit-Junger  notation.  Thus: 

sin  kfflx=sin  kffl  (x'+|)=sin  k^ (x '  +Lx)  =sin  1^— -  cos  ~~  *  cos  1^X-  sin  ~  (D42) 


where  k  = 
m 


nnr 

£ 


m  =  0,1,2,3**** 


x 

It  is  clear  from  Figure  11  that  m  =  1,2,3***  represents  odd  modes 
with  respect  to  the  origin  at  x=0.  However  ia  =  1,3***  odd  and  m  =  2 ,4*  •  * 
even  represent  even  and  odd  modes,  respectively,  with  respect  to  the  origin 
at  O'.  In  accordance  with  the  Feit-Junger  postulation.  Equation  (Al),  we 
wish  to  retain  the  even  modes  only,  i.e.,  m  =  1,3***  odd. 

Equation  (D42)  shows  that  for  the  even  modes,  m  =  1,3* ••  odd,  the 
first  term  in  the  right  member  is  zero  so  that  as  a  special  case 


.  ,  .  nnrx'  .  nnr 

sm  k  x  -*■  cos  —  sxn  -z- 
m  zL  z 


m  odd 


-*■  cos 


(2m  2L  (-1)ID  =  t-1)m,cos  kmX* 
x 


(D43) 

m'  =  0,1,2***  (D44) 


where  m  =  2m'  +  1,  an  odd  number  and  k'  = 
’  m 


(2m'  +  1)tt 
2L 


Thus  Equation  (D44) ,  which  contains  even  modes  only  in  the  axial 

system  with  origin  O’,  represents  only  half  of  the  modes  of  sin  kmx  in 

the  axial  system  with  origin  0.  Moreover  in  the  axial  system  with  origin 

O',  if  we  take  the  mirror  image  of  (i.e.,  reflect)  the  even  modes  m=3,  7, 

11,  etc.  about  the  x'  axis,  then  all  even  modes  will  have  a  positive  value 

in  ® 

at  x*  =0.  Mathematically  this  is  accomplished  by  letting  (-1)  -*■  1.  Thus 

for  even  modes  with  respect  to  the  origin  O'  and  positive  values  of  sin  kmx 


at  O' 


sin  k  x 
m 


cos  - - zr—} - =  cos  k'x'  which  is  identical  to  the 

2L  m 

x 


corresponding  factor  in  Equation  (Al). 
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APPENDIX  E 

BOLT,  BERAXEK  AND  NEWMAN  METHOD 


i 


DERIVATION 


Several  researchers  for  Bolt  Beranek  and  Neman  have  presented 
Mathematical  relationships  for  the  added  mass  and  submerged  natural  fre¬ 
quency  of  a  vibrating  finite  rectangular  plate  radiating  into  half-space. 
The  fornulations  are  briefly  considered. 

Reference  10  gives  the  following  relationship  between  the  fire- 

* 

quency  and  wavelength  of  a  vibrating  subaerged  structural  panel 


-  2/  pf  \  ,.4  2  2 

(  ppV7 "  p  C£  ; 


k  >k_ 
P  f 


(El) 


4  2  2  2 

From  Equation  (iv.S.16)  of  Reference  8,  k  k  c»  =  a  ,  so  that 

T)  X  Mil 


U _ =  U 

am  an 


P 

r,  Pf  t1/2 

Id  : 


k  >  kr 
P  £ 


(E2) 


Equation  (E2)  is  identical  in  for*  to  Equation  (A34) . 

References  11  and  12  give  the  following  relationships  as  precise 
for  waves  on  a  large  flat  submerged  plate. 


(p_  A/M  \  '  /  P~  \ 

1+-V-^ 1  +~M~  J 

h  ’  V  liV' 


-1/2 


-1/2 


Apf  /  Pf  ' 

where  M  =  Ap  h  and  M  , ,  ,  =  - - so  that  M  =  Ap  h  II  +  — -v -v 

o  s  added  k,  mn  s  y  Ps\” 


(E3) 


Hence 


P  "1/2 

“mn  =  “  L1  *  ^l]  ;  kh 


>  k. 


(E4) 


Equation  (E4)  is  identical  in  form  to  Equation  (A34) . 


In  Equation  (11)  of  Reference  10,  let  co  -*■  comn  and  in  Equation  (iv.5.16) 


of  Reference  8,  let  w  -»•  w  . 

’  m  mn 
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Wavelength 

Mass  density  of  water 

Density  (mass  per  unit  volume)  of  plate 

Coefficient  (reactance)  depending  on  aspect  ratio  a/b 

Reactive  (added  mass)  component  of  impedance  for  the  m 
mode  of  vibration 

Reactive  (added  mass)  component  of  impedance  of  a 
rectangular  piston 

Circular  frequency  for  the  mnt^1  mode  of  vibration 


Note:  In  Equations  (F5) ,  (F6),  (F7),  and  (F8)  the  subscript  m  -*■  ran  for  con¬ 


sistency  with  notation  in  the  previous  appendixes. 
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S5&a*atti 


DERIVATION 

References  13  and  14  present  the  following  methods  to  account  for 
the  added  mass  of  a  rectangular  plate  vibrating  in  water. 

Method  1  (see  Reference  13) 

Reference  13  quotes  Reference  15  as  the  source  of  the  following 

equation  for  the  natural  frequency  of  the  first  mode  of  a  simply  supported 

* 

rectangular  plate  vibrating  in  water. 


^11  ^  water 


(F, ,) 

11  vacuum 


|/i+1'X  ^p7.85h 


(FI) 


The  function  x(a/b)  shown  in  Figure  12  has  been  derived  in  Reference  16 
for  a  plate  which  is  clamped  on  the  edges  y  =  o,  b  and  simply  supported  on 
the  edges  x  =  0,  a.  This  equation  can  be  used  to  obtain  the  order  of 
magnitude  of  the  correction  due  to  the  added  mass  for  both  simply  supported 
and  clamped  plates.  However  the  reader  is  referred  to  Reference  15  for  a 
more  accurate  analysis  of  added  mass. 


Method  2  (see  References  13  and  14) 

References  13  and  14  present  the  following  analysis  for  the 
determination  of  the  added  mass  per  unit  area. 

Divide  the  plate  into  equal  elemental  areas  and  let  [P  (x,y)l  be 

th  m  ^6 

the  average  pressure  on  the  y  elemental  area  A  due  to  modal  vibration 

■ffi  I 

of  the  6U  area  Ag  (see  Figure  13).  Then 


Since  (F  )  = - — - -a-^— n-  (see  Appendix  A  or  Reference  15),  then 

11  water  Vl+(Ma/Mp) 

the  present  author  deduces  that  the  added  mass  per  unit  area 

Ma  =  £i'x(j~)  7 " ^5h J  mp  where  mp  is  the  structural  mass  per  unit  area. 
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Masrrraafiga 


i.  .-v 


[P  (x,y)]  „  = 
1  nr  J  Y6 


(Z  )  [w(x,y,t)] 

Yg  m  Jm 


(F2) 


The  total  pressure  on  the  Ytb  area  due  to  the  vibration  of  all 
other  elements  will  be 


(Z  J 

K  -  1  -23  KJ 


m 


6im 


(F3) 


In  these  equations,  the  mutual  radiation  impedance  which  is  a  function 

Yo 

of  the  frequency  and  the  distance  between  the  elements  may  be  written 


<V»  =  PocoV(V»  *  i(V„l 


(F4) 


The  reactive  or  added  mass  component  (Xy^m  is  a  function  of  the 

aspect  ratio  a/b  and  the  nondimensional  parameter  nb/A  =  <o  ' b/2co  (note 

is  denoted  by  P  in  References  13  and  14).  The  impedance  of  a  rectangular 

piston  is  computed  in  Reference  17.  For  relatively  low  frequencies 

irb/A  <  1,  the  reactance  (added  mass  component)  of  the  piston  can  be  written 

(for  consistency  with  previous  results  in  Appendixes  A-F  exclusively,  we 

let  w  -»■  w  ) : 
m  mn 


.  w  b 

Y  =  £(•£■)  Zk  -  f  rZ)  _ _ . 

XP  V  X  V  2c 


(F5) 


* 

In  the  analysis  presented  in  References  13  and  14,  leading  to  the  com¬ 
putation  of  the  plate  natural  frequencies,  the  sums  involving  8  and 

(Equation  (F3))  are  quite  complex  and  require  considerable  computation. 

Therefore,  it  was  assumed  that  the  average  pressure  on  the  y  area  can  be 
computed  by  considering  the  entire  plate  to  act  as  a  rectangular  piston 
with  a  deflection  equal  to  the  average  of  wm  over  the  plate  (see 

Reference  13  for  additional  details).  This  simplified  the  mathematical 
analysis. 


I 
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>i 


where  f(a/b)  is  a  function  of  the  aspect  ratio  and  can  be  obtained 

directly  from  Figure  14.  The  circular  natural  frequency  of  the  plate  in 
13* 

water  is  then 


(w  1 

■  mn  water  on 
one  side 


(w  ) 

mn  vacuum 


[  P  b  a  A?. 
P 


i,3=m,n 


The  A.^  and  B.^  depend  on  the  beam  functions  used  to  represent  the  mode 

shapes  and  therefore  depend  on  the  boundary  conditions  of  the  plate.  The 

values  of  and  B_  for  several  modes  of  plates  which  are  clamped  or 

simply  supported  on  all  edges  are  shown  in  Table  2. 

Plates  with  combinations  of  boundary  conditions  can  also  be  solved 

so  long  as  the  mode  shape  can  be  approximated  by  a  product  of  beam 

functions.  For  example,  consider  a  uniform  plate  which  is  clamped  at  the 

longer  edges  (0,b)  and  simply  supported  at  the  shorter  edges  (0,a).  (This 

13 

is  the  case  given  in  Reference  16.)  Then  for  the  first  mode 


A?..  =  (0.8309  x  0.6366) 2  =  0.2798 


B. .  =  1  x  0.5  =  0.5 


(a)  ) 

mn  water  on 
one  side 


(to  ) 
mn  vac 


"p  b~| 

_2_J  0.2798  f  (g) 


In  accordance  with  the  previous  footnote,  we  deduce  that  the  added  mass 

2 


fp  b  A.  .“I 

2m  V  B.  . 

L  d  nJ 


per  unit  area  f(jj)  m  for  water  on  one  side.  For  water  on 

a  L  p  ij-»  p 

both  sides,  we  double  this  value  in  Equations  (F6)  and  (F7)  and  (F8) .  The 
symbol  y  used  in  Reference  13  to  denote  the  plate  mass  per  unit  area  has 
been  replaced  here  by  the  symbol  m  . 
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For  a  steel  plate  vibrating  in  water 


(w  ) 

mn  vac 


(F8) 


Equation  (F8)  is  of  the  same  form  as  Equation  (FI)  originally  found 
in  Reference  16.  For  most  of  the  practical  cases,  the  correction  of  the 
frequency  due  to  added  mass  given  by  Equations  (D8)  and  (Dl)  will  be  of 
the  same  order  of  magnitude. 


APPENDIX  G 
LEIBOWITZ  METHOD  III 


NOTATION 


The  notation  included  for  Appendix  F  also  applies  to  this  appendix. 

> 

DESCRIPTION 

The  approximate  expression  presented  is  convenient  for  the  compu¬ 
tation  of  the  added  mass  and  (low)  natural  frequencies  of  a  vibrating 
rectangular  plate.  The  development  of  this  expression  is  based  on  the 
work  of  Reference  14  discussed  in  Appendix  F. 

DERIVATION 


Ignoring  plate  damping  and  considering  clamped  or  simply  supported 
* 

plates,  the  natural  frequency  of  a  vibrating  plate  in  water  is  given  by 
Equation  (F6) .  Following  Reference  14,  we  take  as  a  low  frequency 
approximation 


1/2 


Xp  =  0.48k  [ab]'  =  0.48 


to  „  0) 

mn  ,  ,  ,1/2  -  mn  ▼  n  .  . 
-  [ab]  for - <  1 

o  o 


(Gl) 


Substituting  Equation  (Gl)  in  Equation  (F5) ,  we  get 

f(g)  =0.96  (|)1/2 

which  is  an  approximate  fit  to  Figure  14. 


(C-2) 


? 


J. 


P 

% 


Reference  14  also  permits  a  solution  of  the  frequencies  for  the  more 
general  case  of  a  rotationally  constrained  beam  including  plate  radiation 
and  structural  damping  in  the  frequency  equation. 
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* 

Substituting  Equation  (G2)  in  Equation  (F6) ,  we  get 


(u  ) 

v  mn  vacuum 


(o)  1 

v  mnwater  on 
one  side 


y 


i,j=m,n 


0.48  Po (ab) 1/2  A2. 

m  5T7 

P 


(G3) 


where  A. .  and  B. .  are  values  given  in  Table  4  for  clamped  and  simply 
13  xj 

supported  plates. 


*  0.48  p  (ab)1/2  A?. 

The  added  mass  per  unit  area  M  =  - - -  ^  m  for  water  on 

r  am  e.  .  p 

P  il 

one  side.  For  water  on  both  sides,  this  quantity  is  doubled  in 
Equation  (G3) . 


'-'v 
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APPENDIX  H 


tsasss’e- : 


NOTATION 


PROCEDURE  FOR  MODIFYING  THE  MAESTRELLO  PROGRAM  TO 
INCLUDE  EFFECTS  OF  FLUID  LOADING  (OPTION  3) 


A. 

l 

A. . ,  B. . 
13 


a,  b 


a  ,  a 
mn  mn 


D’  °d 


E 

e 


Constants 

Quantities  which  depend  on  the  beam  functions  used 
to  represent  the  mode  shapes  (see  Table  4) 

Dimensions  of  plate  (see  footnote  to  w  of 

aluminum  plate  in  this  Appendix  for  more  precise 
description) 

Plate  modal  damping  in  air  and  water,  respectively 
Speed  of  sound  in  fluid 

Compressional  wave  velocity  of  the  plate  equal  to 

,1/2 


~  E  T 
.P.Cl-v2). 


P.Cl-v 

- 

Specific  damping  energy  at  any  stress  a  and  at  peak 
stress  respectively  in  a  part  under  nonuniform 

stress  (0  <  D  <  Dd) 

Young's  modulus 

Equal  to  2.718;  base  for  natural  or  Naperian  system 
of  logarithms 

Equals  =  ” 


mn 


K. 

1 

k 

k 


Frequency 

Plate  natural  frequency 

Function  of  aspect  ratio 

Panel  thickness 
Mode  numbers 
Constants 

Acoustic  wave  number  equal  to  co/c 

7  2 

Surface  wave  number  equal  to  [(lMT/a)'  +  (mr/b)  ] 
for  Method  I;  for  Method  II,  a  and  b  may  or  may  not 
be  interchanged  (see  footnote  to  00  of  aluminum 
plate  in  this  Appendix) 
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M,  m 

P 


M* 


M  ,  m 
a  mn 


m,  n 

W  ,  W 
p  mn 


K* 


P(«) 

T 

p 

q 


u,  u 

9  00 

u 


Plate  mass  per  unit  area  in  air 

Total  mass  per  unit  area  of  plate  in  water;  equal 
to  sum  of  plate  and  added  masses  per  unit  area 

Added  mass  (or  apparent  mass  or  virtual  mass)  per 
unit  area 

Mode  numbers 

Plate  weight  per  unit  area  in  air  and  added  weight 
per  unit  area,  respectively 

Total  weight  per  unit  area  of  plate  in  water;  equal 
to  sum  of  plate  and  added  weight  per  unit  area 

Power  spectrum 

Mean  square  turbulence  pressure 

Equal  to  k/ks  (see  footnote  to  q  of  aluminum  plate 
in  this  Appendix) 

Free-stream  velocity 

Convection  velocity 


(x,y),  (x',y') 
a 


3 


6* 

6  ,  S 

mn  mn 


n 


n 


s 


e 


V 


Points  on  the  panel  at  which  displacements  are 
measured 

Equal  to  1  for  fluid  loading  on  one  side  of  plate 
only;  equal  to  2  for  fluid  loading  on  both  sides  of 
plate;  dimensionless  damping  energy  integral  (see 
section  on  steel  plate  parameters  in  this  Appendix) 

Total  damping  coefficient  of  plate  in  a  fluid;  di¬ 
mensionless  strain  energy  integral  (see  section  on 
steel  plate  parameters  in  this  Appendix) 

Critical  damping 

Boundary  layer  displacement  thickness 

Total  damping  ratio  of  plate  in  air  and  water, 
respectively 

Material  loss  factor;  equals  y  -  y1,  lateral 
partial  separation 

Loss  factor  of  a  specimen  or  part 

Eddy  lifetime  for  steady  convection,  i.e.,  time  in 
which  value  of  correlation  coefficient  obtained 
from  envelope  of  correlation  maxima  (maxima- 
maximorum)  drops  to  1/e 

Kinematic  viscosity  of  fluid  near  wall 


K 


Equals  x  -  x’,  longitudinal  partial  separation 
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K 


P,  P 


a,  o. 


T 

T 

W 

<t>  (x,y),  <P  (x',y‘) 

Mnn  mn  1 

w 

0)  ,  (0 
mn  mn 


Radius  of  gyration 

Miss  density  of  fluid  mediiun 

Mass  density  of  plate 

Amplitude  of  reversed  stress  and  maximum  value  of 
stress  in  a  part  (0  <  a  <  a^) 

Equals  t  -  t*,  time  delay 

Local  wall  shear  stress 

Plate  eigenfunctions 

Circular  frequency  equal  to  27rf 

Plate  modal  frequency  in  air  and  water, 
respectively 
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The  Maestrello  program  for  determining  the  turbulence-induced  vi¬ 
bration  and  radiation  of  plates  is  modified  and  then  extended  to  include 
the  effects  of  fluid  loading.  The  modified  program  is  used  to  obtain  the 
vibratory  response  for  a  fluid-loaded  simply  supported  rectangular 

aluminum  isotropic  plate  and  for  a  fluid-loaded  simply  supported  rectangular 

* 

steel  isotropic  plate.  The  effects  of  hydrophone  size  and  boundary  layer 
thickness  are  excluded  for  simplicity  and  other  plate  boundary  conditions 
on  the  response  are  not  considered.  However,  the  user  can  correct  for 
these  effects  and  treat  various  plate  boundaries  as  well  as  plate  curvature 
by  including  Option  1  (Reference  2)  and  Option  2  (Reference  3)  in  the  com¬ 
putations  (see  DISCUSSION  AND  EVALUATION) .  The  methods  for  determining  the 
input  data  used  in  the  computations  are  also  described. 

MATHEMATICAL  ANALYSIS 

The  results  of  the  analytical  study  for  the  fluid  loading  of 
rectangular  plates  presented  in  Appendixes  A-G  are  summarized  in  Table  1. 

Of  these  results,  only  the  uncoupled  (dominant)  modes  (see  DISCUSSION  AND 
EVALUATION)  are  considered  here  for  the  modification  of  the  original 
Maestrello  program  which  excludes  fluid  loading.  With  this  practical 
restriction.  Table  1  shows  that  there  are  then  two  basic  methods  for  compu¬ 
ting  the  virtual  or  added  mass  and  the  associated  natural  frequencies  of  a 
fluid  loaded  plate.  The  equations  corresponding  to  the  methods  designated 
as  Methods  1  and  2  below  are  now  summarized  and  modified  for  purposes  of 
practical  computation. 

Method  I :  Based  on  Analyses  in  Appendixes  A-E 


(added  mass  per  unit  area)  =  - x — •  =  -  -  - -■  "  1  . 


&S  ”k  >  k  [1-£)2]1/2  V1-^' 

S  K 

S 


where  0  <  y  =  q  <  1 . 


That  is,  the  excitation  function  which  represents  measured  turbulence 
data  is  not  corrected  for  hydrophone  size  in  the  present  computations. 
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Hence  in  the  original  Maestrello  program  when  fluid  loading  is  in¬ 


cluded. 


M  -*■  M*  =  m  +m  =m  |  1  + 
p  mn  p 


1 


ap 


k  m  (1  -  q 

sp  M 


!)1/2J 


(H2) 


also 


0)  =  w 

mn  mn 


r  ap  1  1/2  r,  ap  l"1/2 

1  *  2  .2.1/2  .  “un  1  *  r  ,,  .2-11/2 

L  Ps(ks-k)  hj  L  pkhl-(|-)  J 

L  '  s  '  - 


(H3) 


=  u 


ap 


”  1  Psksh[l  -  q2]1/2 


-1/2 

] 


In  summary,  when  fluid  loading  is  included,  the  original  Maestrello 

program*  is  modified  in  accordance  with  Method  I  by  letting  M  M*  and 

a)  -  u  where  M*  and  w  are  given  by  Equations  (H2)  and  (H3) , 
mn  mn  mn  &  J  n 

respectively. 


Method  2:  Based  on  Analyses  in  Appendixes  F  and  G 


M 


V  A2 

[ap  b  A.  .  -i 

— —  f  (r-)  U  ;  i,j=m,n 

2m  ij  J  P 


(H4) 


Values  for  f(a/b)  are  obtained  from  Figure  14.  For  practical  computation  a 

finite  number  of  values  of  f(a/b)  versus  a/b  can  be  tabulated  and  stored  in 
* 

the  computer.  The  numbers  of  such  values  should  be  sufficient  to  allow 
the  computer  to  calculate  required  intermediate  values  with  an  acceptable 
accuracy  by  means  of  linear  interpolation  of  the  stored  values.  For  the 


* 

In  this  method  of  computation,  b  always  represents  the  longest  side  of 
the  plate.  Thus  b  may  lie  along  either  x  or  y  and  in  the  direction  of  or 
orthogonal  to  the  flow. 
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In  summary,  when  fluid  loading  is  included,  the  original  Maestrello 
program'*'  is  modified  in  accordance  with  Method  2  by  letting  M  M*  and 
where  M'  and  comn  are  given  respectively  by  Equations  (H5)  and 
(H6)  in  general  and  by  Equations  (H8)  and  (H9)  for  low  frequencies. 


METHOD  FOR  DETERMINING  INPUT  DATA 


The  methods  for  determining  the  input  data  required  for  computer 
calculations  are  now  described.  Generally  the  methods  are  similar  to  those 
given  in  Appendix  B2  of  Reference  1.  Eitner  the  data  are  arbitrarily  pre¬ 
scribed  by  the  user,  i.e.,  the  values  are  chosen  to  represent  the  range  of 
interest,  or  the  selections  may  correspond  to  experimental  or  analytically 
determined  values  for  a  parameter.  However,  the  form  and  types  of  basic 
(raw)  data  available  for  use  in  computing  some  of  the  input  data  as  well  as 
some  of  the  features  of  the  response  computations  for  the  aluminum  and 
steel  plates  differ  somewhat.  In  consequence,  different  methods  are  used 
to  evaluate  certain  of  the  input  data  for  these  plates.  The  logical 
presentation  is  then  to  furnish  for  each  plate  a  description  of  the 
particular  methods  used  for  determining  the  data.  Tabulations  of  the  actual 
computed  input  data  used  in  response  calculations  are  also  given  for  each 
plate. 

The  following  input  data  are  furnished  to  the  computer. 

Flow  data:  U  ,  t  ,  6*  =  FU  =  FU  ,  w,  p^  9,  A. ,  K.  where  i  =  1,2,3  and  a,  q 

C  W  XI 

Panel  data:  a,  b,  h,  Smn,  5jn,  E,  M,  M',  4>mn,  u>mn,  wmn,  Z,  n,  t,  m,  n,  x, 
y,  x’,  y' 

Aluminum  Plate 

The  values  for  the  input  data  are  tabulated  in  Tables  2a  and  2b. 

The  methods  used  in  determining  these  values  are: 


Parameter 

Description 

A.,  K. 
i’  i 

Prescribed  constants  used  in  Equation  (B7)  of 

Reference  1 

a,  b,  h 

Prescribed  quantities 

E 

A  prescribed  quantity 

M  =  m 

P 

A  prescribed  quantity 

M' 

Using  Method  1  only  for  the  aluminum  plate,  M’  is  com¬ 
puted  by  means  of  Equation  (H2) 

m,  n 

Prescribed  data 

T 

P 

Equals  P(oj)  dw  where  P(w)  is  obtained  from 

•'o 

Equation  (B7)  of  Reference  1.  This  quantity  can  also 

I 
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aeasured  directly ._  For  the  present  calculation  it  is 

2 

convenient  to  set  p  =  1  so  that  the  autocorrelation 
function  for  the  turbulence  pressures  and  the  cross 
correlation  function  for  the  panel  displacements  given 
by  Equations  (B8)  and  (B52)  respectively  of  Reference  1 
can  be  regarded  as  normalized  to  this  value 

2 

(We  can  also  compute  p  .  Thus  for  subsonic  speeds 
1  2 

Jacobs  shows  that  p  *3.1  T  and  from  Reference  20, 

w  -1/5 

for  a  smooth  plate  T  /q  ~  0.060  (R  )  =  0.060 

-1/5  "  2  x 

(U^x/v'j  where  q  =  1/2  p  is  the  free  stream  dy¬ 
namic  pressure.) 

U  A  parameter  whose  values  are  prescribed  by  the  user, 

c  Uc  =  0.8  U^;  for  the  present  problem,  Uc  =  8,  16,  32, 

64  ft/sec 

x,y,x',y'  Prescribed  points;  the  cross  correlation  of  the  displace¬ 

ments  are  computed  for  these  points 
*  -1  /5 

6*=(FU)=(FUoo)  Equals  0.37x/8  (U^x/v)  '  (see  Equations  21.6  and  21.8 

of  Reference  20),  Tables  1.1  and  2.1  of  that  reference 
give  values  of  v  in  air  and  water.  Using  this  equation, 
values  of  6*  for  a  given  fluid  can  be  prescribed  over  a 
Range  of  U^.  For  an  isolated  plate,  the  equation  shows 

that  for  a  given  value  of  x,  6*max  occurs  for  Uot  equal 

to  the  minimum  prescribed  value,  i.e.,  U  =  (UJ  .  ,  and 

for  v  =  v  .  .  For  all  U  >  (U  )  .  the  corresponding 
air  00  00  mm  r 

values  of  6*  <  6*  for  both  v  =  v  .  and  v  =  v  .  . 

max  air  water 

However,  we  can  also  treat  the  displacement  thickness  of 

a  plate  in  a  structure,  acting  as  a  baffle.  The  plate 

can  then  be  moved  sufficiently  rearward  of  the  edge  of 

the  structure  at  zero  incidence  so  that  for  a  given 

position  along  the  plate  and  all  Um  >  (U  )  .  ,  6*  = 

00  mm 

constant  >  <5*  for  both  v  .  and  v  ^  .  In  effect  we 

max  air  water 

are  considering  series  of  plates  of  identical  geometry 

so  located  in  a  structure  that  for  all  U,  >  (UJ  .  the 

displacement  thickness  at  a  fixed  position  on  one  of  the 
plates  in  the  series  is  always  equal  to  6*  =  constant 


Due  to  a  typographic  error,  this  equation  was  incorrectly  written  as 
0,37/8  (U l^x/v)  in  Appendix  B2  of  Reference  1. 
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for  v  =  v  .  and  v  =  v  _  .  To  reduce  the  aisount  of 

air  water 

computation  for  the  aheximen  plate  only  (i.e.,  to  elimi¬ 
nate  5*  as  a  parameter)  we  select  for  any  of  the  pre¬ 
scribed  values  of  U  =  U  /O.S  and  for  v  =  v  .  and 

c  air 

v  =  vKater  the  plate  whose  position  is  such  as  to  yield 

a  displacement  thickness  6*  =  0.033  >6*  .It  is 
r  max 

assumed  that  on  a  given  plate  6*  is  sensibly  constant  in 

the  direction  of  flow  along  the  plate  length 

Total  damping  ratio  in  air  including  the  damping  in  the 
plate  as  well  as  radiation  damping.  In  the  present 

problem  the  radiation  damping  in  air  is  considered  to  be 

* 

much  smaller  than  the  material  damping.  According  to 
Jacobs  the  damping  used  for  the  aluminum  plate  in  air 
was  experimentally  determined  from  an  analysis  of  band- 
widths  between  half-power  points  on  response  resonances 
from  segmented,  unpublished,  high-resolution  deflection 
power  spectral  density  measurements  (see  page  5  of 
Reference  21,  pages  24-25  and  34  of  Reference  22,  and 
page  273  of  Reference  1) 

A  plot  of  the  measured  results  given  in  References  21 


and  22  indicates  that  6 

n 

a  =6  u  12  —  15ir 
mn  mn  mn 


=  15/fmn  approximately.  Hence 


Total  damping  ratio  of  panel  in  water  including  the 
damping  in  the  plate,  the  damping  associated  with  fluid 
loading  and  the  radiation  damping.  However ,  in  the 
present  problem  the  contribution  of  the  radiation  damping 

** 

to  o  is  excluded.  From  page  41  of  Reference  !  we 
rrm 

deduce  that 

a  5  W 
3  mn  mn  . 

a  =-Hir=  - xn  air 

mn  2M 


5  w 

mn  mn 


in  water 


In  any  event  the  experimentally  measured  damping  will  include  the  con¬ 
tributions  of  all  damping  mechanisms  influencing  the  measurement  results. 
The  mathematical  representation  of  damping  in  the  equations  of  plate  motion 
include  these  results  as  a  viscous  damping  term. 

** 

We  are  interested  here  in  the  radiation  damping  of  a  mode  in  the  range 
kp  >  k.  For  heavy  fluid  loading  the  corner  mode  contribution  to  radiation 

damping  is  applicable  at  all  frequencies  whereas  the  edge  mode  contribution 
is  significantly  applicable  at  high  frequencies.7  Simple  equations  for 
determining  these  contributions  are  given  in  Reference  7  and  will  appear  in 
a  companion  report. 
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The  first  equation  represents  the  viscous  damping  of  the 

plate  in  air  and  agrees  with  Equations  (A7)  and  (B18)  of 

* 

Reference  1  if  we  ignore  the  hysterestic  and  acoustic 
radiation  terms  respectively  in  these  equations.  The 
second  equation  which  includes  the  contribution  of  both 
viscous  damping  and  fluid  loading  (also  called  added  or 
virtual  mass  or  liquid  coupling)  agrees  with  Equation  (A51) 


of  Reference  1  if  the  hysteretic  and  radiation  damping 
term  in  this  equation  is  excluded.  It  is  of  interest  to 
note  that  when  radiation  damping  is  not  considered  the 
increased  mass  associated  with  fluid  loading  reduces  the 
modal  damping  term  in  water  (see  Equation  (AS1)  of 
Reference  1)  so  that  a__  <  a _ .  As  shown  in  Tables  3a 


mn 

.2 


mn 


and  3b,  =  0.6  lb/ft  for  the  aluminum  plate.  From 

the  ratio  of  the  above  equation  we  get 


0  Corresponds  to  the  time  in  which  the  value  of  the 

measured  correlation  coefficient  of  the  fluctuating 
pressures  at  the  wall,  obtained  from  the  envelope  of  the 
correlation  maxima,  drops  to  1/e.  Plots  of  0  versus  Mach 
number  for  broad-  and  narrow-band  frequencies  are  given 
by  Maestrello  in  Reference  23,  Figure  5.  The  Maestrello 
narrow  band  measurements  of  eddy  lifetime  for  fre¬ 
quencies  centered  at  1200  Hz  were  extrapolated  to  zero 


* 

The  user  interested  in  including  a  hysteretic  damping  term  as  a  separate 

contribution  to  a  and  a  can  easily  determine  this  term  by  use  of 
mn  mn  J  J 

Equations  (A7)  and  (A8)  and  the  relevant  term  in  Equation  (A51)  of 

Reference  1,  substituting  M*  for  M  and  u>  for  a)  when  fluid  loading  is 

mn  mn 

considered;  see  also  Equations  (13)  through  (15)  of  Reference  13  for  similar 

results  in  terms  of  both  a  and  6  .  However  to  simplify  the  present  com- 

mn  mn  r  r 

putations  and  in  accordance  with  the  Maestrello  procedure  discussed  in 
Appendix  B  of  Reference  i,  explicit  division  of  the  damping  into  its  hy¬ 
steretic  and  viscous  components  is  not  made  here. 


— V 
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K,  n,  r 


T 


W 


<j>  (x',y') 

mnv  1 


0) 

(0 

mn 


Mach  number  at  which  0  =  2.25  x  10  5  sec.  For  all  U 

c 

under  consideration,  (Mach  number)water  <  (Mach  num- 

ber) air  =  Uoo/Cair  5  64/(0. 8)  (1129)  <  0.08.  Since  the 

extrapolated  curve  shows  that  0(Mach  numbers  <  0.08) ~ 
0(Mach  number  equal  zero)  then  approximately 

_3 

0  =  2.25  x  10  sec  for  all  U  considered. 

c 

Prescribed  data 


Determination  of  this  quantity  is  based  on  the  law  of 
the  wall  which  is  further  discussed  on  page  62  of 
Reference  1.  The  Maestrello  measurements  presented 

in  Figure  1  of  Reference  23  indicate  that  /pVr 

=  2.9  to  3.5  for  M  =  0.35  to  M  =  0.75.  In  reasonable 
agreement  Jacobs  finds  that  the  results  of  various 

investigators  yield  VP^/T  =  3.1  as  an  average  value 

for  all  subsonic  Mach  numbers  (see  pages  301  and 
302  of  Reference  1) .  Using  the  latter  for  the 


n 

present  calculations  we  have  t  =  yp/3.1  =  1/3.1^ 
2  w 

0.323  lb/ft  ;  see  also  description  of  computation 
for  p2 . 


Data  required  for  the  computer  program  are  calcu¬ 
lated  by  the  digital  computer  for  a  range  of  pre¬ 
scribed  values  of  m,  n,  x,  y,  x',  y*. 

Prescribed  in  Equation  (B7)  of  Reference  1  to  obtain 

P(«) 

For  a  plate  of  given  geometry,  boundary  conditions, 
and  structural  properties,  this  quantity  can  be  com¬ 
puted  by  the  methods  of  Option  2  (Reference  3) . 

In  the  present  problem  values  for  comn  were  computed 

by  use  of  the  IVarburton  program^  for  a  simply 

supported  plate  and  substantiated  by  means  of  the 

* 

simple  frequency  expression 

■  «* [(r/  *  (r)2] 


"k 

In  Method  I,  a  and  b  lie  along  x  and  y  respectively  and  are  therefore 
identified  with  m  and  n  respectively.  In  Method  II,  b  is  always  the  longer 
side  so  that  in  general  a  and  b  may  lie  along  either  x  and  y  or  y  and  x 
respectively  and  are  correspondingly  identified  with  either  m  and  n  or  n 
and  m  respectively.  If,  in  particular,  a  and  b  are  identified  with  n  and  m, 

2  2 

then  for  Method  II  (only)  the  equation  w  =  kc^  [(mir/b)  +  (nu/a)  ]  should 
be  used. 
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to  Using  Method  I  only  for  the  aluminum  plate  to  is 

'  “  computed  by  means  of  Equation  (H3)  1101 

a  Equal  to  1  for  fluid  loading  on  one  side  of  the  plate 

only;  equal  to  2  for  fluid  loading  on  both  sides  of 
the  plate 


0) 


* 

See  footnote  on  previous  page. 


122 


Steel  Plate 


The  values  foi-  the  input  data  are  tabulated  in  Tables  3a-3f.  The 
methods  used  in  determining  these  values  are: 


Parameter 

A.,  K. 

1  1 

a,  b,  h 
E 

M  =  m 
P 

M* 


m,  n 


Description 

Prescribed  constants  used  in  Equation  (B7)  of 
Reference  1 

Prescribed  quantities 

A  prescribed  quantity 

A  prescribed  quantity 

Using  Method  I,  M*  is  computed  by  means  of 
Equation  (H2);  Using  Method  II,  M'  is  computed  by 
means  of  Equations  (H5)  and/or  (H8) 

Prescribed  data 


Equal  to 


r 


P(a))du  where  P(w)  is  obtained  from 


Equation  (B7)  of  Reference  1.  This  quantity  can  also 
be  measured  directly.  For  the  present  calculation  it 

~2 

is  convenient  to  set  p  =  1  so  that  the  autocor¬ 
relation  function  for  the  turbulence  pressures  and 
the  cross  correlation  function  for  the  panel  dis¬ 
placements  given  by  Equations  (B8)  and  (B52) 
respectively  of  Reference  1  can  be  regarded  as 
normalized  to  this  value 


x,  y,  x’,  y' 


(We  can  also  compute  p  .  Thus  for  subsonic  speeds 
1  ~2 

Jacobs  shows  that  p  ~  3. 1  T  and  from  Reference  23 

W  -1/5 

for  a  smooth  plate,  t  /q  ~  0.06  (R  )  '  = 

-1/3  w  x  2 

0.060  (U^x/v)  '  where  q  =  1/2  p  U^  is  the  free- 

stream  dynamic  pressure.) 

A  parameter  whose  values  are  prescribed  by  the  user; 
Uc  =  0.8  U.  For  the  present  problem  =  8,  16,  32, 

64  ft/sec 


Prescribed  points;  the  cross  correlation  of  the  dis¬ 
placements  are  computed  for  these  points 


123 


I 


6*=(FU)S  (FUJ 


The  displacement  thickness  is  determined  by  use  of  a 
computer  program  available  at  the  Center  entitled 
XG  75,  Fortran  IV  by  R.W.  Brown.  The  program  uses  a 
relation  developed  by  Mangier,  i.e.,  the  Mangier  inte¬ 
gral  transformation,  which  permits  reduction  of  the 

calculation  of  axially  symmetrical  boundary  layers  on 

* 

arbitrary  bodies  of  revolution  to  that  in  two- 
dimensional  flow.  The  Mangier  method  relates  the 
distance  along  the  axis  of  a  body  of  revolution  to  the 
distance  along  a  flat  plate  at  which  the  boundary 
layer  thickness  is  identical.  The  boundary  layer  dis¬ 
placement  thickness  is  then  calculated  using  an  ex¬ 
pression  which  is  a  function  of  the  flat  plate  distance 
and  flat  plate  Reynolds  number  (using  kinematic  vis¬ 
cosity  of  air  at  90  deg  F  and  of  water  at  39  deg  F) , 

24 

developed  by  Granville.  The  derivation  of  the  ex¬ 
pression  is  based  on  similarity  arguments.  The 
undetermined  coefficients  in  the  expression  are 
evaluated  by  use  of  experimental  data.  The  required 
input  data  for  the  program  are  the  axial  distance  (in 
feet)  and  the  radius  of  the  body  (in  feet)  at  that 
position  as  well  as  the  free  stream  velocity  (in  knots). 
The  limitations  of  the  program  are  those  due  to  the 
assumptions  required  for  the  Mangier  transformation 
and  to  the  fit  of  the  available  data  which  Granville 
used  in  his  theory.  The  assumptions  involved  are 
considered  to  have  a  greater  bearing  than  the  data 
fit  on  the  accuracy  of  computation.  In  particular, 
we  observe  that  Mangier  assumes  a  two-dimensional 
flow  which  does  not  explicitly  account  for  the  local 
pressure  gradient.  The  saving  feature,  with  regard 
to  this  omission,  is  the  fact  that  the  boundary  layer 
growth  is  sufficiently  slow  so  that  the  imprecision 
in  the  calculation  is  considered  to  be  approximately 

10  percent.  Figure  15  gives  the  results  of  the  compu- 

* 

tation  for  5*  of  the  actual  structure  in  water  and 
in  air  obtained  by  use  of  the  Brown  computer  program. 

Total  damping  ratio  in  air  including  the  damping  of 
the  plate  as  well  as  radiation  damping.  In  the 
present  problem  the  radiation  damping  in  air  is  con- 

sidered  to  be  much  smaller  than  the  material  damping. 


The  actual  structure  under  consideration  is  a  cylinder  which  is  part  of 
a  cigar-shaped  body  of  revolution.  Our  computations  are  made  above  the  ring 
frequency  which  allows  the  cylinder  to  be  treated  approximately  as  a  plate. 
The  supports  at  the  cylinder  or  equivalent  plate  are  essentially  simple. 


See  second  footnote  for  <5  of  aluminum  plate. 
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8*  IN  INCHES 


In  Reference  25,  Figure  36.9  used  in  conjunction  with 
Figure  36.15  shows  that  the  loss  factor  q  of  various 
steels  at  low  and  intermediate  reversed  stress  values 
and  with  cyclic  stress  values  well  below  the  fatigue 
limit  range  approximately  from  a  maximum  of  about 
0.005  to  a  minimum  of  0.001.  Analytically,  we  can 
also  obtain  a  value  of  this  order  of  magnitude  by  use 
of  Equation  (36.13)  of  this  reference,  namely. 


n  *  JL  A 

S  TT  „  2 
ad 


—  In  this  equation,  E  for 


steel  is  known  and  Dd  the  specific  damping  energy 
associated  with  the  peak  stress  level  reached  anywhere 
during  the  vibration  (i.e.,  the  value  of  D  correspond¬ 
ing  to  a  =  a^)  can  be  obtained  from  published  curves 

such  as  Figure  36.15  and  Figure  36.17  of  Reference  25 
or  tabulated  data.  For  materials  with  uniform  stress 
distribution  a/B  =  1  whereas  for  various  specimens 
with  variable  stress  distribution  a/B  is  given  by 
Figure  36.9  of  Reference  25.  For  a  rectangular  beam 
with  either  a  constant,  linear,  or  quadratic  moment 
distribution,  at  low  and  intermediate  stress, 
a/8  ^0.66.  Additional  details  on  this  procedure  as 
well  as  a  sample  calculation  are  given  in  Reference  26. 

27-34 

Experimental  data  for  the  loss  factors  of  iso¬ 
tropic  flat  plates  of  steel  and  other  materials,  with- 

* 

out  substantial  stress  concentration,  over  the  range 
of  frequencies  100  to  1000  Hz,  also  vary  approximately 
from  about  0.005  to  0.001.  Considering  the  scatter 
of  these  data,  a  reasonably  approximate  empirical  ex¬ 
pression  for  the  loss  factor  (over  this  range  of  fre¬ 


quencies)  is  ri  =  0.5/f 


Hence,  B/8C  =  0/2  =  0.25/fmn 


=6  /2  =  a  /w  . 

mn  mn  mn 

a  =  0.5ir. 


Therefore,  6  =  0.5/f 


The  damping  of  the  panel  modes  is  a  combination  of 
viscous  and  structural  damping  where  the  viscous  damp¬ 
ing  force  is  proportional  to  velocity  and  the 


Thus  the  data  refer  to  the  loss  factor  for  the  bare  (free)  plate  or  to 
a  simple  supported  plate  if  the  corresponding  stress  concentration 
associated  with  the  support  is  not  large.  The  value  of  q  for  the  particular 
boundary  condition  under  consideration  should  be  used  in  computations.  For 

35 

example,  from  data  on  clamped  plates,  Bies  determined  that  n  =  77/f.  This 


1  21-22 

result  is  close  to  that  used  by  Jacobs  ’  for  a  clamped  aluminum  plate. 

7 

Davies,  in  his  investigation  of  plate  response  and  radiation,  also  assumed 
that  n  =  B/f  where  8  is  a  constant. 
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I 


hysteretic  damping  force  is  proportional  to  dis¬ 
placement.  Measurements  have  shown  that  the  damping 
of  panels  is  small  so  that  the  damping  can  be 
represented  adequately  in  most  cases  by  considering 
it  to  be  entirely  viscous. 


Total  damping  ratio  of  panel  in  water  including  the 
damping  associated  with  fluid  loading  and  the 
radiation  damping.  However 3  in  the  present  problem 
the  contribution  of  the  rad.iation  damping  to  5  is 

excluded.  Following  the  analysis  and  argument  made 
for  the  aluminum  plate  and  noting  that  W  =  20.4  lb/ 

ft^  for  the  steel  plate  as  shown  in  Table  3,  then 


(  20.4  \ 

l  20.4  +  W  ) 
\  mn  / 


/  20.4  \ 

1 

/  20.4  \ 

120.4  +  W  ) 
\  mn  / 

2  f 

mr-> 

1  20.4  +  IV  / 
\  mn  / 

£>,  n,  t 


Corresponds  to  the  time  in  which  the  value  of  the 
measured  correlation  coefficient  of  the  fluctuating 
pressures  at  the  wall,  obtained  from  the  envelope  of 
correlation  maxima,  drops  to  1/e,  Plots  of  0  versus 
Mach  number  for  broad-  and  narrow-band  frequencies 
are  given  by  Maestrello  in  Reference  23,  Figure  5. 

The  Maestrello  narrow-band  measurements  of  eddy  life¬ 
time  for  frequencies  centered  at  1200  Hz  were  ex¬ 
trapolated  to  zero  Mach  number  at  which  0  =  2.25  x 

_3 

10  sec.  For  all  Uc  under  consideration,  (Mach 

number)  .  <  (Mach  number)  .  =  U  /c  .  =  U  /0.8c  . 

water  v  air  03  air  c  air 

<  64/ (0.8) (1129)  <  0.08.  Since  the  extrapolated  curve 

shows  that  0(Mach  numbers  <  0.08)  ^0(Mach  number 

equals  zero),  then  approximately  0  =  2.25  x  10  3  sec 

for  all  U  considered, 
c 

Prescribed  data 

Determination  of  this  quantity  is  based  on  the  law  of 
the  wall  which  is  further  discussed  on  page  62  of 
Reference  1.  The  Maestrello  measurements  presented  in 


See  first  footnote  for  6  of  aluminum  plate. 

mn  r 
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vs 


Figure  1  of  Reference  23  indicate  that  Vp  /x  =2.9 

to  3.5  for  Mach  number  =  0.35  to  Mach  number  =  0.75. 

In  reasonable  agreement,  Jacobs  finds  that  the  results 

/“ 2 

of  various  investigators  yield  V  p  /x  =  3.1  as  an 

average  value  for  all  subsonic  Mach  numbers  (see 
pages  301  and  302  of  Reference  1).  Using  the  latter 


for  the  present  calculations  we  have  x  =  y  p  /3.1  = 

2  w 

1/3.1  ss  0.323  lb/ft;  see  also  description  of  p 

4>  (x ,y) ,  <*>  (x',y*)  Data  required  for  the  computer  program  are  calculated 

by  the  digital  computer  for  a  range  of  prescribed 
values  of  m,  n,  x,  y,  x',  y' 

u  Prescribed  in  Equation  (B7)  of  Reference  1  to  obtain 

P(u>) 

oj  For  a  plate  of  given  geometry,  boundary  conditions 

11,11  and  structural  properties,  this  quantity  can  be  com¬ 

puted  by  the  methods  of  Option  2  (Reference  3) .  In 
the  present  problem  values  for  to  were  computed  by 

use  of  the  Warburton  program^  for  a  simply  supported 
plate  and  substantiated  by  means  of  the  simple  fre¬ 


quency  expression 


id  =  KC„ 
mn  i 


[(F)2  <¥)2] 


Using  Method  I,  co^  is  computed  by  means  of 

Equation  (H3);  Using  Method  II,  tomn  is  computed  by 

means  of  Equations  (H6)  and/or  (H9) 

Equal  to  1  for  fluid  loading  on  one  side  of  the  plate 
only;  equal  to  2  for  fluid  loading  on  both  sides  of 
the  plate 


Equal  to  k/ks  = 


See  footnote  for  a)  defined  for  the  aluminum  plate, 
mn 
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COMPUTER  PROGRAMS 


Four  computer  programs  associated  with  the  mean  square  displacement 

computations  performed  in  this  report  for  a  simply  supported  fluid-loaded 

* 

plate  subject  to  turbulence  excitation  are  described  below. 

Program  A':  Modification  of  the  Maestrello  subprogram  A  (designated 

TURAD) ,  described  in  Reference  1 ,  for  computation  of  the  mean 

square  displacement.  The  modification  which  incorporates 

certain  corrections,  additions,  deletions,  and  improvements 

in  efficiency  in  conjunction  with  its  running  on  the  CDC  6700 

supersedes  subprogram  A  and  should  be  used  henceforth.  Details 

of  the  modification  are  discussed  below. 

Program  B’ :  Computes  the  fluid-loaded  natural  (or  modal)  frequencies 

as  well  as  k  .  k  ?  added  weight  w  =  m  g,  and  total  weight 
s  s  mn  mn 

w'  =  M'g  =  (m  +  m  )g.  In  the  program  w  =  m  g  WP,  W  = 
p  mn  p  p  mn 

mmn®  W  =  M'g  FW,  and  wmn  -*•  BOMEGA.  This  program 

is  designated  as  Option  3. 

3 

Program  C':  Warburton  program  for  computing  the  plate  natural  frequencies 
f  and  to  in  air  for  both  simply  supported  and  clamped- 
clamped  boundary  conditions. 

Program  D' :  Computes  the  displacement  thickness  for  an  arbitrary  body  of 
revolution.  The  displacement  thickness  is  used  to  calculate 
the  variable  FUCSQ  in  Program  Af.  Program  D’  was  developed  by 
Mr.  R.W.  Brown  (see  previous  section  of  this  Appendix). 


((MTURAD)) 


Computes  turbulent 
joundary  layer  thickness. 
Symbol  is  5* 


Computes  wmn. 

Fluid  loadinq  routine 


Computes  wmn  by  the 

Warburton  Method. 
Air  loading  routine 


As  previously  discussed,  computations  for  clamped-clamped  boundaries 
can  also  be  performed  with  the  same  program.  The  program  can  also  compute 
the  modal  acoustic-power  radiation  of  a  plate  in  a  reverberant  field. 

These  computations  are  not  made  here. 
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As  shown  in  the  preceding  chart,  when  fluid  loading  is  omitted 
Programs  C'  and  D'  generate  data  for  use  in  Program  A’.  The  methods  for 
determining  additional  imput  data  have  been  described  in  the  previous 
section. 

When  fluid  loading  is  included.  Program  C*  generates  data  for 
Program  B'  and  then  Programs  B’  and  D'  generate  data  for  use  in  Program  A’. 
Determination  of  additional  data  inputs  have  previously  been  described. 

The  computer  and  running  times  of  a  calculation  associated  with  each 
program  are  given  in  the  following  chart. 


Program 


Computer 
(at  NSRDC) 

Time 

Cost  for  Total 
Computation 

CDC  6700 

Using  the  Simpson  rule  of  integration 
the  computer  running  time  is  approxi¬ 
mately  5  min  for  obtaining  a  curve 

of  Y2(u>)  versus  U  for  each  set  of  mode 
c 

numbers  (m,n)  and  four  convection 
velocities  Uc 

$75.00 

IBM  7090 

1.0  min  for  (m,n)  ranging  from  1,1 
to  10,10,  i.e.,  100  values  of  d>mn 

4.00 

IBM  7090 

1.3  min  for  (m,n)  ranging  from  1,1 
to  10,10,  i.e.,  100  values  of  u>mn 

5.50 

IBM  7090 

1.1  min  for  198  values  of  axial 
distance,  flat  plate  distance,  local 
flat  plate  Reynolds  number,  and  6  =  86* 

4.60 

Program  A' 

To  make  the  original  Maestrello  subprogram  A,  designated  TURAD,  a 
more  efficient  program  and  to  enable  it  to  run  on  the  CDC  6700,  the 
following  changes  have  been  made  in  the  coding  of  the  original  program  (the 
modified  program  A'  is  designated  MTURAD) : 

1.  Reduction  of  the  four-dimensional  array  IXYZ  to  a  three- 
dimensional  array 

*IXYZ  (M,N,IT,L2) - ►  IXYZ  (N,1,L2) 

2.  Increase  in  the  dimension  of  the  following  variables: 

*IXYZ  (10,1,50,3) - ►IXYZ  (10,1,3) 

IYZ  (20,1) - ►IYZ  (20,10) 
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G3  (11,1) - ►  G3  (101,10) 

YY  (11) - -  YY  (101) 

TEMPI  (11) - ►TEMPI  (101) 

3.  Rearrangement  of  certain  data  cards 
*FUCSQ 

FM2 

4.  Deletion  of  some  original  data 
*NP  -  no  longer  needed 

PARAM  -  no  longer  needed  (see  below) 

THETA  -  no  longer  a  dimensional  variable 

5.  Correction  in  coding  of  Simpson's  rule.  Original  coding  is: 

*525  G2  (J,M)  =  G2  (M,M)*2 

should  read 

*525  G2  (J,M)  =  G2  (J,M)*2 

6.  Elimination  of  double  precision  functions  due  to  no 
declaration  of  such  at  the  beginning  of  the  program.  Such  occurrences  were 
on  cards  numbered  2310  and  2330. 

2 

7.  Introduction  of  a  tolerance  as  an  option  for  Y  (to)  for  any  plate 
material.  An  exceeded  tolerance  produces  an  error  message  to  be  printed  in 
the  output.  If  use  of  the  option  is  not  desired,  the  user  reads  in  the 
value  0  for  the  variable  IOPT.  If  the  option  is  used  then  any  number  other 
than  0  is  read  in  for  the  variable  IOPT. 

8.  Rearrangement  of  lines  of  coding  due  to  rearrangement  of  data. 
The  lines  are  on  cards  numbered  0620-0670,  0750  and  0760  in  the  original 
listing. 
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For  completeness  we  include  an  explanation  of  the  procedure  by  which 
0  and  amn  were  obtained  by  the  computer  program  in  Reference  1.  For  the 
present  program,  the  method  given  in  Section  2  of  this  Appendix  for  the 
determination  of  0  and  a^  has  superseded  this  procedure.  However,  if  the 
user  finds  the  earlier  procedure  useful  in  computation,  it  can  be  rein¬ 
troduced  into  the  present  program. 

From  page  7  of  Reference  44  and  from  Section  3  of  DISCUSSION  AND 
EVALUATION  (Equation  7b) ,  the  following  condition  obtains  ac  hydrodynamic 
coincidence 


w  0  =  k  U  0  =—-  U  0 
mn  tan  c  Ac 
mn 


We  define  PARAM  =  2U0/A  .  Hence  ui  0  =  ir  •  PARAM  and 

mn  mn 


7T  ♦  PARAM 


it  •  PARAM 


PARAM 


*  ■’  te'  •  ©■]  ■  •  [(;)’  •  ©1 


Now  Uc  is  the  convection  velocity  along  x  (i.e.,  flow)  direction  only. 
Therefore  at  hydrodynamic  coincidence,  the  trace  wave  speed  of  the  free 
plate  bending  wave  is  matched  to  the  convection  speed  in  the  flow  direction, 
yielding  a  greatly  increased  response.  Thus, 


U 

c  k 


=  kcj"[(t) 


K  C„  Tf  — 

a  a 


so  that  (neglecting  the  quantity  (n/b)  in  the  expression  for  0) 


PARAM  _  PARAM  •  a 


U  2L 
c  a 


Maestrello  (Figure  21  of  Reference  45  or  Figure  18  of  Reference  46) 
has  shown  experimentally  that  the  maximum  vibratory  response  (i.e.  maximum 
mean  square  displacement)  occurs  at  hydrodynamic  coincidence  where  the  tur¬ 
bulence  and  plate  modal  wave  numbers  are  matched  for  a  constant  frequency. 
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Maestrello  has  also  performed  computations  which  show  that  the 
condition  ir  ♦  PARAM  =  w  0=1  represents  the  aerodynamic  coincidence  con¬ 
dition  for  whicli  maximum  excitation  of  the  panel  occurs;  see  Figure  7  of 

Reference  1.  This  figure  also  gives  results  for  w  0  <  1  and  w  0  >  1 

&  &  mn  mn  * 

corresponding  to  the  conditions  below  and  above  coincidence,  respectively. 

Thus,  selection  of  the  value  PARAM  for  a  given  plate,  value  of  Uc 
and  mode  number  m  will  yield  the  corresponding  value  of  0.  For  PARAM  =  1/ir, 
the  value  of  0  corresponds  to  the  coincidence  condition.  For  PARAM  <  1/ir, 
the  value  of  0  corresponds  to  the  condition  below  and  above  coincidence, 
respectively. 

In  subprogram  A,  Maestrello  chooses  the  magnitude  of  amn  according 

to  the  co0  region  of  the  curve  that  interests  him.  For  example.  Figure  7  of 

2  -1 
Reference  1  plots  Y  (w)  against  w0.  In  the  region  w0  =  10  to  1.0,  he 


uses  a  /10;  however  for  u)6  =  1.0  to  40,  he  uses  a  . 
mn  mn 

methods  for  determining  a^  are  given  in  Reference  1. 


The  Maestrello 


For  0  «  1  the  modal  mean  square  displacement  is  inversely  pro¬ 
portional  to  the  total  damping  and  at  the  peak  is  inversely  proportional  to 
the  square  of  the  damping.  In  this  region  the  effect  of  damping  on 
~2 

Y  («)  is  greatest.  For  tomn  0  >  1  coincidence  is  not  possible.  In  this 

__ 

region  the  effect  of  damping  on  Y  (w)  is  much  smaller,  see  References  45  and 
46. 
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TABLE  6 


Input  Data,  Computer  Listing,  Flow  Chart,  and  Column  Headings 
for  Input  Forms  on  Data  Cards  for  Updated  Maestrello  Pro- 
gran  A’  Used  to  Compute  Mean  Square  Displacement  of  Plate 
with  and  without  Fluid  Loading 

TABLE  6a 

Input  Required  for  Program  A'  (MTURAD) 


(Units  are  given  in  foot-pound-seconds) 


Data 

Description 

Format 

Symbol  Used 
in  Program 

Unit 

Flow  Characteristies  (Program  A' -an  updated  version  of  Subprogram  A  given 
in  Reference  1) 

Uc 

Broadband  convection 
velocity 

F10.0 

UC(I) 

ft/ sec 

7 

Mean-square  wall- 
pressure  fluctuations, 
which  vary  with  Uc* 

F10.0 

PB2*DPB2(I) 

(lb/ft2)2 

(FUc)2 

/6*Uc\ 

Quantity  l-g— ) 

squared  where: 

6*  =  boundary  layer 
displacement 
thickness 

U  =  free  stream  velocity 

F10.0 

FUCSQ 

ft2 

K.j  ,K2»I<2 

Universal  constants: 

K1  =  0.470 

K2  =  3.0 

K3  -  14.0 

F10.0 

AK 

A-|  ,^2  >^2 

Universal  constants: 

A]  =  1.6 

A2  =  7.2 

A3  =  12.0 

F10.0 

AN 

2  2 

(*PB2  would  represent  a  unique  value  of  p  if  p  were  independent  of  U  .  It 

2  c 

enters  the  program  (i.e.,  data  cards)  once  only.  Since  p  actually  varies  with 

U  ,  a  correction  factor  DPB2(I)  is  entered  with  every  value  of  U  .  Thus 

p2  as  a  function  of  Uc  is  accounted  for  by  the  quantity  PB2*  DPB2(I). 
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Table  6a  (Continued) 


Description 

Format 

Symbol  Used 
in  Program 

Unit 

Data 


Plate  Characteristics 


Panel  thickness 


F10.0 


11 


ft 


*  (W)2 


Square  of  total 
weight  per  unit  area 


F10.0 


FW2 


lb/ ft 


a,b 


Lengths  of  panel 
sides 


F10.0 


ZUP.YUP 


ft 


**  6,6 

mn’  -mn 


Total  damping  ratio 


F10.0 


DAMP 


U)  ,  W 

mn  mn 


Modal  frequencies 
of  the  plate 


F10.0 


OMEGA 


rad/sec 


Additional  Quantities 


110 

mow 

no 

MUP 

110 

DM 

no 

MSTEPS 

no 

NLOH 

no 

no 

no 

NSTEPS 

Range  of  plate 
mode  numbers  for 
which  calculations 
are  desired 


First  m  mode  number, 
last  m  mode  number, 
interval  between  m 
mode  number,  total 
number  of  m's 
m  £  20,  MSTEPS  £  20. 
Same  information  as 
■previously  described, 
with  respect  to  n 
mode  numbers  n  t  10, 
NSTEPS  <  10- 
To  run  program  A‘ 
for  m  >  20,  n  >  10  the 
size  of  the  dimen¬ 
sioned  variables 
IXYZ,IYZ,IY,G2,G3,YY, 
TEMP .OMEGA ,DAMP , FA , 

FC, EIGEN  and  FUDGE 
must  be  examined  and 
increased  accord¬ 
ingly 


Time  delay 


HO.O 


TAU 


sec 


Number  of  values 
of  U.  to  be 

calculated 


110 


KUC 


V*o 


Coordinates  of  a 
point  on  plate  at 
which  mean  square 
displacement  and 
acoustic  power  are 
calculated 


F10.0 


X0,Y0 


ft 


Total  weight 
area 


6_ 


mn 


W  =  M' •  g  where  M'  =  total  mass  per  unit  area.  We  stress 
that  the  user  submits  inertial  data  in  terms  of  weight  per 
unit  area  and  the  computer  program  converts  these  data  to 
mass  per  unit  area. 

Total  damping  ratio  of  plate  in  water;  6mn  =  total  damping 
ratio  of  plate  in  air,  u>mn  =  modal  frequencies  of  plate  in 
water;  u>mn  =  modal  frequencies  of  plate  in  air. 
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Table  6a  (Continued) 


Data 

Description 

Format 

Symbol  Used 
in  Program 

Unit 

Any  point  on  plate 

different  from  x„,y^ 
o  ■7o 

F10.0 

X0P,Y0P 

ft 

Calculated  Output 


4>, 


m,n 


Value  of  eigen¬ 
functions  of  mean 
square  displacement. 

A  value  of  EIGEN  is 
computed  for  each 
mode  (m,n)  with 
three  values  of  total 


damping;  1/10  a, 

10am  „ 
m,n  m,n 


m,n’ 


E16.8 


EIGEN 


'a  ,a 
m,n’  m,n 


Values  of  total 
damping  associated 
with  each  mode 
(m,n) 


E16.8 


FA(m,n,l)  for 
computation; 
A(m,n,DAMP) 
in  output 


1/sec 


Vol 


m,n 


]fz  ? 


Volume  under  each 
eigenfunction 


E16.8 


VOL 


I(m,n) 


Triple  integral  of 
Equation  (B52)  of 
Reference  1 ; 
integral  of  cross 
correlation 


E16.8 


IXYZ 

I(m,n) 


in. 


y2M 


Mean  square  displace¬ 
ment  for  values  of 

1/10  am  ;  am  n; 
m,n  m,n 


E16.8 


ANS 


10a 


m,n 


in. 


m,n 

am,n 


Total  damping  for  each  mode  of  fluid  loaded  plate. 
Total  damping  for  each  mode  of  plate  in  air. 


uuu 


€£SjClu 


r'XWOKSSirf- 


TABLE  6b 

Computer  Listing  for  Updated  Maestrello  Program  A* 


FTN(T) 

tea. 


PROGRAM  MTURAD! INPUT ♦OUTPUT *TAPE5*INPUT »TAPE6»01JTPUT) 

REORGANISED  PROGRAM  FOR  COMPUTING  TRIPLE  0  10 

INTEGRAL  ?0 

USING  SIMPSONS  RULE  30 

COMMENT  TRANS  =  0.0  MEANS  NO  TRANSDUCER  EFFECTS  INCLUDED 
COMMENT  THIS  RUM  OF  MTURAD  HAS  SIMPLY-SUPPORTED  ******* 

C*#*«***#**FREOUENCIES  W//UT  TRANSDUCER  *)**********«********♦******#**** 

DIMENSION  VOL 1 20 #10 )  AO 

DIMENSION  F3 (3 ) »AN5(3 )  50 

DIMENSION  IXYZ(10,1,3)»IYZ(20#10) ,IY(10).G?(401#20) ,G3(101»10) • 

1YY ( 101 ) «TEMP1(101 ) #G5(4000) 

DIMENSION  AK(3)»AN(3)»TITLE(7) 

DIMENSION  OMEGA (20# 10 ) »DAMP ( 20#10)  0100 

DIMENSION  FA(20#10#3)#FC(20#10»3)#cIGEN(20#10#3) » FUDGE (20#10#3)  0110 

DIMENSION  UC(20) #OPB2(20)  0120 

REAL  IXYZ.IYZ.IY  0130 

INTEGER  DM«DN  OlAO 

READ (5«13 )  TAU#TRANS 
13  FORMAT (2F10. 01 

WRITE(6#30l)  TAU  0170 

801  FORMAT! i HI, /5H  TAU=E15.61  0180 


READ(5,1212>  I OPT 

1212  FORMAT ( 15 1 
WRITE(6#2213)  I OPT 

1213  FORMAT (1HO»5HIOPT=I5) 

READ(5»103)H»XO*YO»XOP#YOP#PB2#AK  #AN#(TITLE(I1#I=1#71 
103  FORMAT (12(F10.P/)#(7A10) J 

WRITE(6#20l)  (TITLE ( I ) » 1=1 #7) 

201  F0RMATI7A10) 

WRITE! 6«2n3)PB2  »AK#AN#XO#YO# 


1XOP.YOP  250 

203  FORMAT! 1X»8H  RHO-BAR  9H  SQUARFO=E16.8 # 

24H0K1=E16.8#4H  K2»E16.8,4H  X3*E16.8/  0290 

34H  A1=E16.8.4H  A2=E16.8,4H  A3=E16.8/  0300 

44H  X0=E16.8#4H  YO=E16.8.5H  XOP=E16.8»5H  Y0P=E16.8)  0310 

READ (5# 102 1  KUC»(UC( I) »DPR2! II »I  =  1#KUC)  0320 

102  FORMAT! I10/(2F10,o) )  0330 

WRITE(6#204)  (UC( 1 1 #DPB2 ( 1 1 ♦ I  =  l.KUC)  0340 

204  FORMAT ! 1H07X2HUC13X4HDPB2/! 1H  2E16.8))  0350 

99  READ ( 5# 1 )  ZUP»YUP#H»MLOW.MUP#DM,MSTEPS#NLOW#NUP*DN#NSTFPS  0360 

1  FORMAT(3(F10. 0/). (4110))  0370 

WRI  TE(  6,2  )MLOW»MUP#DM«MSTFPS#NLO'i,#NUP  #DN»NSTEPS  0380 

2  FORMAT ( 9H0M  FROM  I5»4H  TOI5#7H  PM=Is,  0390 

1 12H  A  TOTAL  OF  I5#7H  STFPS/9H  N  FROM  15#  0400 

24H  T0I5«7H  DN=I5»12H  A  TOTAL  0FI5.7H  STEPS)  0410 

WRITE! 6,202 1  ZUP»YUP*H  0420 

202  FORMAT ( 3H  A=E16.3#5H  8=E16.8,5H  H=E16.8)  0430 

READ (5,101)  t ( OMEGA (M»N ) ,M=MLOW,MUP#D") ,N*NLOW,NUP#DN) 

101  FORMAT <F10.0)  450 

IF! Nl)P  .GT.  10)  WRITE! 6 ,888 ) 

888  FORMAT! 1H0,78HDIMENSI0NS  FOR  G3»IYY#TFMP1  ARE  EXCEFDFD.  CHFCK  FOUA 
1TIOM  FOR  KtJP  ON  CAPD  1360,) 

WRITE! 6,2n5 ) ( ( OMEGA (M»N) #M»MLOW,muP»Dm) ,  0460 

1N=NL0W,NUP,DN)  470 

205  FORMAT ( 7H00MEGA=/ ( 8E16.8 ) )  0480 

IF!  TRANS. EO#ri.O)  GO  TO  17 

READ ! 5 ,1477)  R #DEL 
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Table  6b  (Continued) 


l 


i 

i 


1477  FORMAT (2F10.0) 

00  1453  J  =  1*3 
1453  AK(J)  =  A KIJ)  *  SQRT(R/DEL) 

17  00  7  N=1 *NUP»DN 

00  7  M=l,fluPtOM 
7  DAMP(M,N)=.01 

READ  (5*104)  ((DAMP(M»N)*M«MLOW*MUP»DM)*N=NLOW#NUP*DN) 

104  FORMAT (F10.0) 

WR I T E  ( 6  *  206 )  ( ( DAMP  ( M  #  N )  * M=M|_OV*Mt )P  *  DM )  * 

1N=NL0W*NUP»DN) 

206  FORMAT (6H00AMP=/ (8F16*8 ) ) 

Nl=8 

READ ( 5*207) THETA 
207  FORMAT ( F10«5 ) 

A=ZUP 

B=vuP 

8  B7  «  AN ( 1 )  /  AK (1 ) 

B8=AN(2)/AK(2) 

B9=AN( 3 )/AK (3 ) 

C1=4.0/(A*R) 

02=3.14159265 
03=1 */C2 
C4=C2»C2 
06=04/ (A*B) 

DO  10  M=MLOW  «  MUP ♦ DM 

DO  10  N«NLOW*NUP»DN 

XM=M 

XN=N 

X0=A/2. 

Y0=B/2. 

IF(M-M/2*2.NE.0)G0  TO  43 
X0=A/2.-A/(XM*2.) 

43  IF(N-N/2*2.NE.O)GO  TO  44 
YO=B/2.-B/(XN*2.) 

44  XOP=XO 
YOP=YO 

DAMP ( M. N ) “DAMP (M«N)/ 10* 

DO  45  L=1 *3 

FA(M,N.L)=t,MP(M.N)/2.«0MEGA(M»N) 

FA(M»N»L)«FA(M.N*L)*0*5 

FC(M.N*L)«OMEGA(M,N)/FACM*N»L) 

FUDGE (M»N*L)=XM*XN*0MEGA(M*N)*(FA(M»N*L)**2  +OMFGA(M*N ) **2  ) 

EIGEN(M»N.LJ=Cl/FUDGE(M»N*U«SIN<XM*C2*XO/A)*SIN(XN*C2»YO/B)* 
IS  I N ( XM#C2*X0P/A ) *S IN ( XN*C2* YOP/B ) 

45  DAMP(M*N)=DAMP<MjN)*10. 

10  CONTINUE 

WRIT£(6.1003) 

1003  FORMAT ( 33H0X0=X0P* Y0=Y0P *AND  THEY  VARY  WITH 
117H  THE  MODE  NUMBERS) 

DO  46  L=1  *3 

WRITE(6»3) ( (E!GEN(M*N*L) *M=MLOW*MUP*DM) *N»NLOW*NUP»DN) 

3  FORMAT ( 7H0EIGEN=/(8E16*8 ) ) 

46  CONTINUE 

DO  47  L“1 *3 

WRITE(6*48)( (FA(M*N*L) *M=MLOW*MUP*DM)*N«NLOW*NUP*DN) 

48  FORMATS 13H0A(M*N*DAMP) */ (8E16*8) ) 

47  CONTINUE 

CALL  VOLUM ( A  *  B  *MLCW  *MUP ♦ DM  *  NLOW  *  NUP  *  DN  * VOL ) 

IC0UNT=1 

READ(5»225)  FUCSO 


138 


490 

500 

510 

530 

0540 

550 

0560 

570 


600 

610 

690 

700 

720 

730 

740 

770 

780 

790 

800 

810 

820 

830 

0840 

850 

0860 

870 

380 

890 

0900 

910 

0920 

0930 

0940 

0950 

0960 

0970 

0980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 


'A. 


■1 


Table  6b  (Continued) 


FORMAT { F10.0 ) 

00  777  KU=l»tCUC 
WRITE(6»303)UClKU) 

FORMAT (4H  UC=*E16.8) 

Bl=AK<l)**2*FUCSO 

B2=AK<2)**2*FUCS0 

B3=AK(3)**2*FUCSQ 

B4=AN(  1  )*AK<1  )*FtlCS0 

B  5 = AN ( 2 1 *AK ( 2 ) *FUCSQ 

B6=AN(3)*AK<3)*FUCSQ 

DO  778  M*MLOW *  MUP  *  DM 

KC0UNT=1 

R£AD(5»222)FW2 

FORMATIFIO.2) 

C5=  < A*B*PB2 ) / <2 .*C4*FW2* ( B7+B8+B9 ) ) 
C5=C5*32.2*32«2*12»*12. 

CON5T=C5*DPB2 (KU) 

WR I TE 1 6  *  304 )  CONST 

FORMAT  (12X.7H  C0NST=»F16.8) 

XM=M 

WRITE(6*780  J  THETA 

FORMAT ( 1H0 ,6HTHETA=  *£12 .5  J 

IF(THETA  .GT.  100)  GO  TO  999 

XLOW=0« 

DX=1./(20.*OMEGA(M*NUP  ) /2./3. 14159265 ) 
IUP*5./0X*THETA+1. 

IF( I UP  .GT.  3599)  IUP=3599 

IF( IUP-IUP/2*2)  5°1»500,501 

IUPsIUP+1 

IUP  =  IUP  +400 

0X=THFTA/40. 

ZLOW=-A 
JUP=20*M+1 
ZJIJP=JUP-1 
DZ=2  »*A/ZJUP 
YLOW=0. 

KUP=10*NUP+1 
YKUP=KUP-1 
DY=B/YKUP 
DO  11  1=1 » IUP 
G5(  I )=1. 

I F (  I.E0.40UG0  TO  11 
IFII.NE.l.ANn. I.NE.IUP)  GO  TO  510 
GO  TO  511 
G5 ( I ) =G5 t I )*2. 

IF ( I-I/2*2.EQ.O>  GO  TO  513 

GO  TO  11 

G5(  I  )=G5( I )*2« 

CONTINUE 

Z=ZLOW 

00  21  J=1»JUP 
XM=M 

D1=XM*C2*Z/A 

02=COS(D1) 

G2 ( J  »M ) =02+1 . / ( XM*C2 ) * ( S I N ( APS ( OJ ) ) -APS <01 ) 
1*02) 

IF( J.NE. l.ANO. J.NF« JUP )  GO  TO  520 
GO  TO  521 

G2(J*M)=G2(J*M)*2« 
l  IF ( J-J/2*2.EO.O)GO  TO  525 
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Table  6b  (Continued) 


GO  TO  21 

1600 

525 

G2(J»M)=G2(J*M)«2. 

21 

Z=Z+DZ 

1620 

Y=YLOW 

1630 

DO  31  K*1  »ICUP 

1640 

DO  30  NsNLOH*NUP»DN 

1650 

XN=N 

1660 

D3=XN»C2*Y/6 

1670 

D4=C0S(D3> 

1680 

G3(K»N)*D4+1./<XN*C2)*(SINID3)-D3*D4» 

1690 

IFOC.NE.l.AND.K.NF.KUP)  GO  TO  530 

1700 

GO  TO  531 

1710 

530 

G3(KtN)=G3(K»N)«2* 

1720 

531 

1FCK-K/2*2.EO»0)  GO  TO  533 

1730 

GO  TO  30 

1740 

533 

G3(K»N)*G3(X*N)*2« 

1750 

30 

CONTINUE 

176'i 

31 

Y=Y+DY 

1770 

Y=YLOM 

1780 

DO  39  K^1»KUP 

1790 

YY(Kl=Y*Y 

1800 

39 

Y=Y+DY 

1810 

DO  40  Ll=l»3 

1820 

DO  40  N-NLOH#NUP»DN 

1830 

40 

IXYZ(N.l»Ll)-0. 

KAPPA=0 

1860 

X=XLOW 

1870 

DO  160  I*1.IUP 

1880 

IFITAU  .EO.O.O)  GO  TO  630 

1890 

El  *  UCCKU»*(TAU-X) 

1900 

GO  TO  632 

1910 

630 

El  *  UC(KUJ*X 

1920 

632 

CONTINUE 

1930 

DO  50  N*NLOH#NUP»DN 

1940 

50 

IYZ(M.NJ»0. 

1950 

Z-ZLOH 

1960 

DO  120  j«l.JUP 

1970 

E2=fZ-El)*(Z-El) 

1980 

E4=81+E2 

1990 

E5=B2+E2 

2000 

E6=B3+E2 

2010 

DO  60  K»1«KUP 

2020 

60 

TEMPI lK)=B4/( E4+YY  <K I J  +B5/ { E5+YY (K) ) +B6/ ( E6+YY ( K )) 

2030 

DO  70  N*NLOW»NUP*DN 

2040 

70 

IYIN)=0. 

2050 

DO  90  K-lfKUP 

2060 

DO  90  N-NLOW»NUP«DN 

2070 

90 

IY(N)=IY(N)+G30C*N)*TEMP1<KJ 

2080 

DO  100  N*NLOW.NUP»DN 

2090 

100 

IYINJ=IY(N)«DY/3» 

2100 

DO  110  N*NLOW»NUP#DN 

2110 

TEMP=IY(N) 

2120 

110 

IYZ(M»NJ=IYZ<M«N>+G2<J»M)*TEMP 

2130 

120 

Z=Z+DZ 

2140 

DO  130  N«NLOW»NUPtDN 

2150 

130 

IYZ(MtN)*IYZ(MtN)*DZ/3» 

2160 

1001 

CONTINUE 

2170 

DO  140  N*NLOW»NUP#DN 

2180 

TEMP=IYZ(M»N) 

2190 

XM=M 

2200 

Table  6b  (Continued) 


XN=N  2210 

F1=XM*XN«C6  2220 

F2=OMEGA!M»N)*ABS!X)  2230 

DO  140  Ll*l»3  2240 

F3!L1)=FA(M#N.L1)*ABS!X)  2230 

IF!F3!Ll)»GT»5«)G0  TO  140  2260 

F6=F1*EXP ! -F3 1  LI ))* ! S IN I F2 ) +FC! M*  N  * LI ) *COS ( F2 1 )  2270 

IFIX  »GT«  5.*THETA)  GO  TO  137 

IFITAU.E0.0.0)  GO  TO  634  2300 

G6=EXP ! -ABS ( ( TAU-X ) /THFTA) I 

GO  TO  636  2320 

634  G6=EXP!-ABS!X/THETA)) 

636  CONTINUE  2340 

G1=2.*F6*G6*G5 ( 1 )  2350 

IXY2(N*1  .LI)  =  IXYZIN.l  »L1)  +  Gl*TEMP*DX/3. 

137  CONTINUE  2370 

140  CONTINUE  2380 

IF!KAPPA.NE.O)GO  TO  160  2390 

I F C I .NE.401 )GO  TO  160  2400 

KAPPA=1  2410 

DX=1./(20.«OMEGA(M»NUP  J/2./3. 14159265) 

GO  TO  1001  2430 

160  X=X+DX  2440 

DO  601  N=NLOW.NUP.DN  2450 

IF! THETA  .GT.  -100)  GO  TO  999 

WRITE (6. 141)  M.N  2470 

141  FORMAT (3H1M=15»3H  N=I5//45X6HI (M.N) 38X»18HHM.N)«FIGEN*C0MST//  2480 

12X5HTHETA4X11H0MEGA*THETA6X11HDAMP=1./10.5X7HDAMP=1.9X8HDAMP=10.  2490 

2fiXllHDAMP=l./10.5X7HDAMP=1.9X8HDAMP*10.) 

DO  602  L2=l .3  2520 

602  ANS1L21=IXYZ!  N.l  .L2)«EIGEN!M.N»L2)*CONST 
T=THETA*OMEGA(M*N) 

WRITE! 6 .142) THETA  »T.(IXYZ(  N.l  »L2) »L2»1.3 ). ( ANSCL2 ) »L2»1 .3)  2550 

142  FORMAT ( 1X.F9.6.E14.6.6E16.8)  2570 

IFdOPT  .EO.  0)  GO  TO  600 
IF( N  .EO.  NLOW )  TH0LD=ANS1 2 ) 

IF!  N  .EQ.  NLOW)  GO  TO  600 

IF( ANS( 2 )  .EQ.  T H0LD/10.  .OR.  ANS(2)  .GT.  THOLD/IO.)  GO  TO  600 
WRITE (6 .998) 


998  FORMAT! 1H0.72HM0DE  IS  MORE  THAN  1/10  BELOW  VALUE  OF  Y**2  OF  OMFGA 
1AT  CORRESPONDING  UC.) 

GO  TO  778 

600  KCOUNT =KCOUNT +1 

IF! KCOUNT  .GT.  NSTEPS)  GO  TO  778 
READ! 5.222 )FW2 

C5= !  A»B«DPf'2 )  / 1 2.*C4*FW2  * !  B7+B8+B9 )  ) 

C5=C5*32.2**2*12.**2 

CONST=C5*DPB2!KU) 

601  CONTINUE 

778  CONTINUE 

779  I COUNT * ICOUNT +1 

IF! I COUNT  .GT.  KUC)  GO  TO  777 
READ  15.225)  FUCSO 


777  CONTINUE  2600 

999  STOP  2610 

END  2620 

SUBROUTINE  VOLUM! A.B.MLOW.MUP.DM.NLOW.NtlP.DN.VOL)  2640 

INTEGER  DM.DN  2650 

DIMENSION  VOL ! 20.10 )  2660 

PI=3. 14159265  2670 


V 


Table  6b  (Continued) 


00  10  N*NLOW»NUP»DN 
XN=N 

00  10  M»MLOW»MUP*DM 
XM=M 

VOL(-M»N)*0. 

IF(N-N/2*2.EC.O}GO  TO  10 
IF(M-M/2*2.EQ*0)GO  TO  10 
GAMMAN® I 2.*XN+1 . J  *PI /2 • 

GAMMAMs  1 2«*XM+1 . >  »PI/2 • 

XKN«S IN  <  G AMMAN /2 • J /S I NH ( 6AMMAN /2  •  > 

XKM«SIN ( GAMMAM/2 * I /SINH ( GAMM AM/2*  > 

VOL ( M  »N ) *16.*A«B/GAMMAM/GAMMAN/ ( 1 *+XKM  J  *144* 

1/ ( 1«+XKN)*S IN (GAMMAM/2* )*SlNI GAMMAN/2*  J 
CONTINUE 

WRITEC6.20H (V0L(M.N)»M«ML0W»MUP#0M)*N®NL0W»NUP»0N) 
FORMAT C 28 HOVOLUME  UNDER  EIGENFUNCTION*// 

K8E16.8)) 

RETURN 

END 

FUNCTION  SINH  CX» 

SINH«0.5*<EXP(XJ-EXP(-XJ ) 

RETURN 

END 


2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 


TABLE  6c 

Flow  Chart  for  Updated  Maestrello  Program  A* 


READ  AND  WRITE 
INPUT 
DATA 


CALCULATE 

VARIABLES  USED 
INSIDE  LOOPS 


CALCULATE 

nhf 7 
3  Av 

2*J»JX  r- 

I  "» 


DOUBLE  DO  LOOP  ON  10 
m  and  n  Vary 
.  .VO  -  A /S.  >0  -  H/2 


IS  EVEN 


TO  -  A/2  -  A/(2m) 


IS  EVEN 


TO  -  H/2  -  R/2n 


DO  on  IS 
Sot  up 
for  .1 


io*MO 


CALCULATE  for  ouch  damp 

r-. 


45  CONTINUE 


10  CONTINUE 


Valuos  juat 
Calculated 


SUBROUTINE  VOLUM 


Calculates  1  per  »,n 


(2m*  l\ir  (2 n*  1\» 

- )~  K,n  l  — ~ — )- 

2/2  \  2  /2 

n-1  r~ 

sin<2m«l)-  sin(2n4l) 

l+l)|(2«+  l)j  U  - -  - 

sinh(2»i4  1)-  sinh(2n4l) 

2J  _ 


sinh  (2i»4  1)~ 
2 


CALL 

SUBR 

VOLUM 


WRITE  answers 
RETURN 


For  convenience,  fluid-loaded  parameters  (represented  by  barred  quantities)  are  not 

distinguished  from  corresponding  parameters  for  the  plate  in  air.  Thus,  for  the  fluid-loaded 

plate  we  must  think  of  certain  parameter  symbols  being  converted  as  follows:  -*■  to  , 

a  -*-a  ,6  -»8  ,  i.e.  to  is  really  to  in  the  water  case,  etc. 

mn  mn  mn  mn’  mn  '  mn 


Table  6c  (Continued) 


Set  counter  for  readinq 
FUSCQ 


Set  counter  far  readinp  FW2 


Set  upper  and  lower  limits  for  x,y,z 
integrals;  Determine  number  of  uteps 
for  each  integral  and  the  step  size 
for  each 


Table  6c  (Continued) 


y  Set  variable  of  integration  =  lower  limit  on  y 
Calculate  function  of  y  values 


Set  up  and  initialize  to  zero  the 
array  for  summing  on  outside 
integral 


Table  6c  (Continued) 


9 


Table  6c  (Continued) 


Table  6d  (Continued) 


AN(1)  10 


AN(2)  1C  80 


KUC  10 


Z'JP  10 


NUMBER  OF  CARDS  NEEDED  FOR  U^,  DPB2  ARRAYS  IS  EQUAL  TO  KUC 


YUP  10 


MLOW  10  MUP  20  DM  30  MSTEPS  40 


Table  6d  (Continued) 


NLOW  10  NUP  20  ON  30  NSTEPS40 


0MEGA(1 ,1)10 


0MEGA(2,1)10 


WITH  THE  ORDER  OF  CARDS  (1,1),(2,1),...(M,1),(1,2),...(M,2)...,(1,N),..(M,N)  TO  COMPLETE 
THE  OMEGA(M,N)  ARRAY 

R  10  DEL  20  80 


DAMP(1 ,1 )  10 


WITH  THE  ORDER  OF  CARDS  (1,1),(2,1),...(M,1),(1,2),..(M,2)...,(1,N),..(M,N)  TO  COMPLETE 
DAMP(M.N)  ARRAY 

THETA  10  80 


FUCSQ(1)10 


FW2(1)  10 


FW2(2)10 


THE  NUMBERING  OF  FW2  IS  DETERMINED  BY  THE  NUMBER  OF  (M,N)  MODES 
FUCSQ(2)10 


THE  NUMBER  OF  FUCSQ  IS  DETERMINED  BY  USER. 
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Program  B' 


A  description  of  Program  B'  (Option  3)  is  given  in  Tables  7a-7d. 


TABLE  7 

Input  Data,  Computer  Listing,  Flow  Chart,  and  Column  Headings 
for  Input  Forms  on  Data  Cards  for  Program  B’  (Option  3) 
Used  to  Compute  Added  and  Total  Weight  per  Unit  Area  and 
Modal  Frequencies  of  Fluid-Loaded  Plate 

TABLE  7a 

Input  Required  for  Program  B' 


Variable 

Name 

Format 

Description 

Unit 

WP 

F10.8 

Weight  per  unit  area  of  plate 

lb/ft2 

RHOF 

F10.4 

Weight  density  of  fluid 

lb/ft3 

RHOW 

F10.4 

Weight  density  of  plate 

IBS! 

H 

F10.4 

Panel  thickness 

ft 

RL1 

F10.4 

Length  of  panel  side,  x-direction 

ft 

RL2 

F10.4 

Length  of  panel  side,  y-direction 

For  a  =  1,  fluid  loadinq  is  on  one  side 

of  plate 

ALPHA 

F10.4 

a  =  2,  fluid  loadinq  is  on  both  sides 

- 

of  plate 

MLOW 

no 

Lower  limit  of  M  mode  m  <  20 

MUP 

no 

Upper  limit  of  M  mode 

NLOW 

no 

Lower  limit  of  N  mode  n  <  10 

NUP 

no 

Upper  limit  of  N  mode 

OMEGA (m,n) 

F10.2 

Array  for  modal  frequencies  of  plate  in 
air 

rad/sec 
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n  r>  r> 


TABLE  7b 

Computer  Listing  for  Program  B' 


SIBJOR  MAP»FIOCS 

$ EXECUTE  IB JOB 

SIBFTC  OMEGA 

DIMENSION  OMEGA ( 20 » 10) *POMEGA( 20*10 ) »ROMEGA( 20*10) *W(20»10) 
READ(5*15)WP 
15  FORMAT (F10. 2) 

READ (5 *20)  RHOF*RHOW*H*RL1»RL2»ALPHA 
20  FORMAT (6F10. 4) 

READ(5»1) MLOW*  MUP  *NLOW »  NUP 
1  F0RMAT(4I10) 

READ ( 5 ♦  10 )  ( (OMEGA (M*N )  *M=MLOW*MUP )  #N=MLOW*NtlP ) 

10  FORMAT  (Fi0«. 4) 

METHOD  1 

WRITE(6*40) 

40  FORMAT(lHl.lOX#8HMETHOD  1) 

0=0. 

DO  244  L=l»4 
DO  25  M=MLOW»MUP 
DO  25  N=NLOW»NUP 

FKS2=( (FLOAT (M )*3»1416) /RLl )**2+( (FL0AT(N)*3.1416) /RL2)**2 
FKS=5QRT ( FKS2 ) 

W ( M  *  N )  =  ( ALPHA*RH0F ) / ( FKS*SQRT ( 1 .-0**2 ) 

FW=WP* ( 1 .+ ( ( ALPHA*RH0F ) / ( FKS*WP*SQRT ( 1 • -0**2 ) ) J ) 
POMEGA(M*N)=OMEGA(M»N) 

BOMEGA  (M»N)  =POMEGA  ( M * N )  *S0RT ( 1 .  /  ( 1  •  +  ( ( A(_PHA*RHOF )  /  ( RH0W*FKS*H*SORT 
1(1. -0**2)))) ) 

WRITE(6»30)M»N*Q 

30  FORMAT ( 1H0  »2HM=  » 12  »2X*2HN=* I 2  »2X»2HQ=  *F5.3 ) 

WR I TE ( 6 ♦ 3  5 ) FKS2 . FKS ♦ FW ♦ B  OMEG A  <  M  »  N ) »  W ( M  *  N ) 

35  FORMAT ( 1H0»5HFKS2=  »F12»4»2X  »4HFKS=»F12*4»  2X  *3HFW=*F12.4  *2X*7HR0MEG 

1A=»F12.4»2X#2HW=»F12.4 ) 

25  CONTINUE 
0=0+.3 

244  CONTINUE 
STOP 
END 


TABLE  7d 


Column  Headings  for  Input  Forms 


j  Data  Cards 


RL2  50  ALPHA  60 


Program  C* 


This  program  yields  a  solution  for  clamped  and  simply  supported 
plates  corresponding  to  the  Warburton  Method  as  described  in  Appendix  1  of 
Reference  3.  Warburton  treats  the  frequency  parameter  subscripts  m,n  as 
the  number  of  nodal  points  along  the  plate  length  and  width  respectively. 
However,  most  other  authors  treat  m,n  as  the  mode  numbers  along  these  di¬ 
mensions  (or  define  it  for  the  opposite  dimensions).  Thus  (m  =  2, 
n  =  3 Warburton  means  the  2  mode  containing  2  nodes  along  x  and  3  along 
y  whereas  (m  =  2,  n  =  3)  ^  ®eans  the  2,  3  mode  containing  either  3  nodes 
along  x  and  4  along  y  or  4  nodes  along  x  and  3  along  y  depending  on  the 
definition  of  m,n  with  respect  to  the  x,y  coordinates.  To  avoid  confusion 
and  for  compatibility  with  most  investigators,  the  program  assigns  the  modal 
(not  nodal)  meaning  to  m,n  for  all  computations. 

In  computing  the  simply-supported  plate  frequencies  by  the  Warburton 
Method  the  value  of  SPEC  must  be  1.0.  In  computing  the  clamped-clamped 
piate  frequencies  by  this  method  any  value  of  SPEC  other  than  1.0  is  used. 

In  all  computations,  the  frequency  f  (Hertz)  is  obtained  as  the  pro¬ 
duct  of  the  frequency  parameter  (or  )  and  a  factor.  For  the 

Warburton  computations  the  factors  are  expressed  as 

hTT_  /  E 
a2  y  48  Pm(l-a2) 

where  the  mass  density  Om  =  RHOW/G 
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TABLE  8 


Input  Data,  Computer  Listing,  Flow  Chart,  and  Column 
Headings  for  Input  Forms  on  Data  Cards  for  Program  C’ 
(Warburton  Program)  Used  to  Compute  Natural  Frequencies 
of  Simply  Supported  and  Clamped-Clamped  Plate  in  Vacuo 

TABLE  8a 

Input  Required  for  Program  C* 


Variable  Name 

Format 

Description 

Unit 

NCASE 

15 

Number  of  cases 

— 

M 

15 

Modes  in  x-direction 
m  <■  20 

— 

N 

15 

Modes  in  y-di recti on 
n  s  10 

— 

A 

F12.6 

Plate  dimensions, 
length  in  x  direction 

in. 

B 

F12.6 

Plate  dimensions, 
length  in  y  direction 

in. 

H 

F12.6 

Plate  thickness 

in. 

E 

E16.8 

Young’s  modulus 

lb/in.2 

SIGMA 

F12.6 

Poisson's  ratio 

— 

RHOW 

F12.6 

Weight  density  of  plate 

lb/in.3 

G 

F12.6 

Acceleration  due  to 
gravity 

in. /sec 

SPEC 

F10.0 

Option  for  obtaining  fre¬ 
quencies  of  either  simply 
supported  or  clamped- 
clamped  plate 

! 

TABLE  8b 

Computer  Listing  for  Program  C' 


SEXECUTE  I 8 JOS 

SIBJOB  MAP*FIOCS 

SIBFTC  WARS _ _ _  _  _  _ 

COMMON  M»N,A»e»H»E#SlGMA«RHOM»PI*G 
C  K  -  MOOES  IN  X  DIRECTION 

C  N  •  MOOES  IN  Y  DIRECTION 

C  A  -  LENGTH  IN  X  DIRECTION 
C  B  -  LENGTH  IN  Y  DIRECTION 
C  H  -  PLATE  THICKNESS 

C  E  -  YOUNGS  MODULUS 

C  SIGMA  -  POISSONS  RATIO 
C  RHOW  -  PLATE  DENSITY 
C  G  -  ACCELERATION  DUE  TO  GRAVITY 

f  WHHHHHHHWWUHHHHHHHHHHHHHHHHHHHHIWHtWHHHHHHHHHHHHMHHHMHIWHMHHHHHHHHW 

PI*3#H1M27 

READ I 5 *2 INCASE 

00  500  L*1»NCASE 

READ*5*2)  M  «N 

READ (5 *3 I  A»B*H 

READ(5»4)E.SIGMA*RH0W*G 

2  FORMAT I 2151 

3  FORMATi3Fl2*61 

4  FORMAT ( £1&«8«3F12  *6) 

RH0M-RH0W/6 
10  CALL  WARS 
500  CONTINUE 
STOP 
END 

SIBFTC  WARBER 

SUBROUTINE  WARB 

REAL  LAMBDA • JX • J Y  »K  tKP 

DIMENSION  OMEGAC20.10) 

DIMENSION  FREQ<25*10)  •  GXUOO)  »HX  C100)  .JXI 100)  #GY(  1001  *HY(  100)  • 

1  JY(IOO) 

COMMON  M«N*A#8»H«E*SIGMA*RH0M*PI»G 
READ! 5 •111) SPEC 
111  FORMAT  C  F10.0 ) 

A2«A*A 

B2»B*8 

A4*A2*A2 

B4«B2*B2 

MP1-M+1 

NP1-N+1 

IFtSPEC.EQ.  1.0)  GO  TO  510 
Gxm«i. 

HXtl)»l. 

JXtl)-l. 

GY ( 1 1*1. 

HY{ 1 1*1. 

JY(1)*1» 

GX{ 21*1.506 

HX( 21*1.248 

JX( 21*1.248 

GY( 21*1.506 

HY( 21*1.248 

JY( 21*1.248 

DO  100  M1«3»MP1 

GX<  Ml)  “FLOAT  (Ml  )***5 

HX ( Ml  1“ ( (FLOAT (Ml )“.5) **2 1*1 l»“2./( (FL0AT<M1)-.5)*PI 1 1 


I 


4 


I 
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Table  8b  (Continued) 


JX(M1)»HX(M1) 

100  CONTINUE 

DO  150  N1*3*NP1 
GY*Nl)*FLOAT*Nl)-*5 

HY*Nlf***FLOAT*Nl)-.5)»*2)**l«-2«/' *FL0AT*N1)-.5)*PI)) 

JY(Nl)*HY(Nl) 

150  CONTINUE 
GO  TO  590 

510  DO  500  Ml  *  1»MP1 

GX(M1)  *  FLO/.T *  Ml )  -  1*0 
MX  (Ml )  »  GXtM)  **2 
500  JX(M1)  ■  HX(M1) 

DO  550  N1  *  1 *MP1 
GY(Nl)  »  FLOAT (N1)*»1#0 
HYCN1I  «  GY‘.N1)**2 
550  JYIN1I  »  HYtNII 
590  WRITE*6»20)A*B*H#E*SIGMA»RHOM 

20  FORMAT* lHlf.3H  A*#F7*2t3H  B*»F7.2.3H  H*«F7.4#3H  E*#E11.4*7H  SIGMA** 

1  F7.2.6H  RMOM**E11.4) 

WRlTE(6*19) 

19  FORHAT * //23X*  22H  WARBURTON  FREQUENCIES) 

IW  *  1 

DO  400  N2'*2*NP1 

N21-N2-1 

WRITE*6*21)N21 

21  FORMAT I 3F  N**I2) 

WRITE*6»22) 

22  FORMAT  *  9?*  *  1HM*  1 5X  »6HLAMBDA  *  1 6X  * 5H  FREO) 

DO  300  M2*2*MP1 

M21*M2-1 

XLAMSQ*CX{M2)*GX*M2)*GX*M2)*GX*M2)+*GY*N2)«GY*N2)*GY(N2)*GY*N2) 

1  *A4)/B<+12»*A2/B2)**SZGMA*HX*M2)*HY(N2)+*1.~SIGMA>*JX*M2)*JY(N2)) 
LAMBDA*  SORT ( XLAMSQ ) 

FREQ (M2  »N2 ) * ( ( LAMBDA*H»P I)/A2) *SQRT ( E  / ( 48 . *RHOM* ( 1 #-S IGMA**2 ) ) ) 

WRITE*  £ »23)M21 »LAMBDA»FREQ* M2»N2 ) 

0MEGA(f»2*N2)  *  2.  *  PI  *  FREQ(M2»N2) 

WRITE*(‘.*30)  0MEGA*M2»N2 ) »  IW 

23  FORMAT J5X»I5*5X»E15*8f5X»El5. 8) 

WRITE* 3*30)  OMEGA* M2 «N2 ) #  IW 

30  FORMAT  * 1X*F15.4«65X» 15) 

IW  *  IW  +  1 
300  CONTINUE 
400  CONTIAUE 
RETURA 
END 
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Table  8c  (Continued) 


SUBROUTINE  WARB 


161 


The  input  data  and  results  are  labeled  and  printed  out  for  each 
value  of  NCASE.  The  mode  numbers  (m,n),  nondimensional  frequency  A,  and 
final  frequency  f  (Hertz)  are  given. 

Program  D' 

This  program  calculates  turbulent  boundary  layer  thickness  for  an 
arbitrary  body  of  revolution. 


TABLE  9 

Input  Data,  Computer  Listing,  Flow  Chart,  and  Column 
Headings  for  Input  Forms  on  Data  Cards  for  Program  D' 
(Brown  Program)  Used  to  Compute  Turbulent  Boundary  Layer 
Thickness  for  an  Arbitrary  Body  of  Revolution 

TABLE  9a 

Input  Required  for  Program  D' 


Variable  Name 


Format 


F12.8 

F12.8 


6F10.5 

6F10.5 


ft/sec  6F12.8 


Description 

Number  of  pairs  of  data  points 
describing  the  body  of  revolution 

Length  of  the  body  of  revolution 

Constant  to  which  all  of  the 
data  points  are  normalized 

Axial  distance  along  the  body 

Radial  distance  of  the  body 
from  the  axis 

Number  of  speeds 

Speed 


on  on 


TABLE  9b 

Computer  Listing  for  Program  D' 

c  mangler  integrals 

C  R.  WILLIAM  BROWN*  INITIALS  CWCt  NSROC  CODE  942*  PHONE  227-1252 

C  THIS  IS  A  PROGRAM  FOR  THE  CALCULATION  OF  THE  TURBULENT  BOUNDARY 

C  LAYER  THICKNESS*  DELTA*  IT  USES  A  RELATION  AFTER  MANGLER  WHICH 
C  RELATES  THE  DISTANCE  ALONG  THE  AXIS  OF  A  BODY  OF  REVOLUTION  TO 

C  THE  DISTANCE  ALONG  A  FLAT  PLATE  AT  WHICH  THE  BOUNDARY  LAYER  THICK- 

C  NESS  IS  IDENTICAL*  THE  BOUNDARY  LAYER  THICKNESS  IS  THEN  CALCULAT- 

C  ED.USING  A  RELATION  DUE  TO  GRANVILLE  (DTMB  REPORT  NO*  1340)*  THIS 

C  EXPRESSION  IS  BASED  ON  FLAT  PLATE  DISTANCE  AND  FLAT  PLATE  REYNOLDS 

C  NUMBER* 

C 

c 

c  688  COORDINATES*  (REFERENCE-GENERATING  FUNCTION— BRIAN  BOWERS) 

C 

C 

ODIMENSION  Xf 500) *  Y<500>*  XZ(300)»  YZ(SCO)*  XZR(500)*  YZR(500)*  YZ 
1R2I500) »  XBARC500)*  RXBAR(500#15) *  DELTAI300*15) »  Utl5) 

C  READ  IN  THE  CONSTANTS  N»  R*  AND  Z* 

C  N  -  THE  NUMBER  OF  PAIRS  OF  DATA  POINTS  DESCRIBING  THE  BODY  OF 

C  REVOLUTION* 

C  R  *  THE  LENGTH  OF  THE  BODY  OF  REVOLUTION. 

C  Z  *  THE  CONSTANT  TO  WHICH  ALL  OF  THE  DATA  POINTS  ARE  NORMALIZED* 

C  READ  IN  THE  DATA  POINTS  DESCRIBING  THE  BODY  OF  REVOLUTION*  X(U 

C  »  THE  AXIAL  DISTANCE  ALONG  THE  BODY.  Y(I)  •  THE  RADIAL  DISTANCE 
C  OF  THE  BODY  FROM  THE  AXIS* 

READ(5*1)  N*  R*  Z 

1  FORMAT  (IS*  2F12.8) 

READC5.2)  ( X ( I ) •  I  *  1*  N) 

READ t 5*2)  (Y(I),  I  ■  1*  Nl 

2  FORMAT (6F10.5) 

C  READING  IN  THE  NUMBER  OF  SPEEDS*  L* 

READ  (5*  18)  L 

18  FORMAT  (15) 

C  READING  IN  THE  SPEED*  U(K)*  IN  FEET  PER  SECOND* 

READ  (5*  19)  (U(K) •  K  *  1*  L) 

19  FORMAT  (6F12.8) 

SETTING  X(l)  AND  Y(l)  EQUAL  TO  ZERO*  AND  SHIFTING  ALL  OF  THE  DATA 
POINTS*  IF  NECESSARY. 

IF  (X( 1) )  15.  8.  9 
9  M  *  N  ♦  1 
TEMPI  «  X(l) 

DO  10  J  “  2*  M 
TEMP2  =  X(J) 

X(J)  =  TEMPI 
TEMPI  =  TEMP2 

10  CONTINUE 
8  IF  (Y(l) )  15*12*  11 

11  M  =  N  +  1 
TEMPI  *  Y(l) 

DO  13  J  *  2*  M 
TEMP2  =  Y(J) 

Y(J)  =  TEMPI 
TEMPI  =  TEMP2 

13  CONTINUE 
X ( 1 )  =  0.0 
Y( 1 )  =  0.0 
N  =  N  +  1 

CONVERTING  THE  DATA  POINTS  WHICH  WERE  RE*D  IN  INTO  ACTUAL  BODY  CO¬ 
ORDINATES. 

12  DO  3  I  =  1,  N 
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r>  r» 


Table  9b  (Continued) 


xzm  =  xm/z 

YZ(I)  *  Y(I)/Z 
XZR ( I )  *  R*XZ( I ) 

YZR(I)  *  R»YZ(I) 

YZR2CI)  »  ( YZR( I ) )**2 
IF! 1-2)  6.  7.  7 

6  XBAR(I)  =  0. 

GO  TO  3 

C  CALCULATING  THE  DISTANCE  ALONG  A  FLAT  PLATE#  XRAk. 

7  XBAR(I)  «  ( 1 • /YZR2 ( I ) ) *S I MPUN ( XZR  *  YZR2»  I) 

3  CONTINUE 

WRITEC6#  25) 

250FORMATC99H1  INPUT  X  VALUES  INPUT  Y  VALUES  AXIAL  DISTAN 

ICE  BODY  RADIUS  FLAT  PLATE  DISTANCE//) 

WRITE(6*5)  (X( I ) •  Y( I ) #  XZR ( I )  #  YZR(I)#  XRAR(I)»  I  *  1#  N) 

5  FORMAT (5F20. 8) 

DO  26  K  *  li  L 
RXBAR! 1*K)  =  0.0 
DELTA! 1»K)  =  0.0 
26  CONTINUE 

DO  17  K  *  1*  L 

CALCULATING  THE  LOCAL  FLAT  PLATE  REYNOLDS  NUMBER*  USING  A  KINEMA¬ 
TIC  VISCOUSITY  OF  WATER  AT  39  DFGREES  F. 

DO  16  I  «  2#  N 

RXBAR! I»K)  =  U(K)*XBAR(I)/0#00001684 
C  CALCULATING  THE  TURBULENT  BOUNDARY  LAYER  THICKNESS*  DELTA# 
IF(RXBAR(I*K) )16. 30*31 

30  DELTA! I. K)  *  0.0 
GO  TO  16 

31  DELTA! I #K)  *  0.0598*XBAR!I)/!ALOG10{RXBAR(I#K) )  -  3.170) 

16  CONTINUE 

17  CONTINUE 

DO  20  K  »  1*  L 
WRITEI6.  24)  U!K) 

24  FORMAT ! 12H1 VELOCITY  ■  .F12.8.  14H  FEET  /  SECOND///) 

WRITEI6*  21) 

210FORMAT ( 116H  AXIAL  DISTANCE  FLAT  PLATE  DISTANCE  LOCAL  FLAT 

1  PLATE  REYNOLDS  NUMBER  BOUNDARY  LAYER  THICKNESS#  DELTA//) 

WRITE!6»22)  (XZR!  I )  ♦  XBARUl#  RXBAR(I#K)#  DELTA!  I*K)»  I  *  1#  N) 

22  F0RMATIF17.8*  F20.8*  18XE15.8#  28XF12.8) 

20  CONTINUE 
15  STOP 
END 
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TABLE  9c 

Flow  Chart  for  Program  D* 


TABLE  9d 


Column  Heads  for  Input  Forms  on  Data  Cards  for  Program  D' 

N  5  R  17  Z  29 


NUMBER  OF  VALUES  ON  CARD  IS  EQUAL  TO  N 
Y(l)  10  Y(2)  20  Y(3)  30  Y(4)  40  Y(5)  50  Y(6)  60 


NUMBER  OF  VALUES  ON  CARD  IS  EQUAL  TO  N 
L  5 


U(l)  12  U(2)  24  U(3)  36  U(4)  48  U(5)  60  U(6)  72 

NUMBER  OF  VALUES  ON  CARD  IS  EQUAL  TO  L 
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COMPUTER  RUNS 


Results  obtained  from  the  computer  programs  for  the  input  data  of 
Tables  2  and  3  are  presented  in  Figures  1  and  2  respectively.  The  figures 
show  computer  runs  for  the  normalized  modc-1  mean  square  displacement  of  the 
turbulence  excited  simply  supported  aluminum  and  steel  plates  with  fluid 
loading  effects  included. 


APPENDIX  I 

HYDROSTATIC  PRESSURE  EFFECTS  ON  NATURAL  FREQUENCIES 


a 

b 


E 

G 

h 

K 


k  ,k 
nr  n 

M 

m,n 


H 

R,R 
’  a 


Plate  length 
Plate  width 

Equal  to  - 


Eh' 


12(1  -  v) 

Young's  modulus 

Constant  dependent  on  k  /k  (see  Reference  37) 

m  n 

Plate  thickness 

Constant  depending  on  aspect  ratio  of  plate;  for  the  panel 
dimensions  cited  in  Appc-lx  I,  K  =  0.0018  (see  Reference  42) 

2 

Effective  wavelength  equivalencing  factor,  dependent  on  oha  /D  (see 
References  36  and  37) 

Mode  numbers  in  the  x-  and  y-directions  respectively 

Mass  per  unit  area  of  plate 
Mode  numbers 
Hydrostatic  pressure 

Radii  of  curvature 


W(x)  Deflection  shape  of  panel  to  loading  by  uniform  pressure 
Wq  Deflection  at  the  center  of  the  plate  due  to  hydrostatic  pressure  P^ 

x,y  Plate  coordinates 

y  Equal  to  1.5 

\  Wavelength 

v  Poisson's  ratio,  generally  taken  to  be  equal  to  0.3 

°x>°y  Stresses  in  directions  associated  with  m  and  n  respectively 

Circular  natural  frequency  of  vibration 


DESCRIPTION 


A  clamped  panel  below  the  water  surface  is  subject  to  essentially 
uniformly  distributed  pressure  which  deflects  the  panel  and  creates  bend¬ 
ing  stresses.  The  natural  frequencies  of  singly  curved  stressed  plates 
are  given  approximately  by^  ^ 


,  ■  or  [. 


2  2 
(k  +  k  ) 
m  n  J 


a  hk 
x  m 


a  hk 
y  n 

D 


where  k  =  - — 
m  K^a 


K  =  — 

*  n  Kb* 


Equation  (II)  is  based  upon  wavelength  equivalencing  which  effectively 

reduces  the  dimensions  of  the  clamped  plate  to  those  of  the  equivalent 

simply  supported  plate  in  each  mode.  For  an  unstressed  flat  plate  this 

equation  yields  values  of  natural  frequencies  which  agree  within  ±10  per- 

3  40 

cent  of  those  calculated  by  the  Warburton  method.  *  In  this  equation,  G 
and  are  considered  to  be  determinable  quantities  (see  Notation) .  Hence 
it  remains  to  determine  the  stress  and  curvature  produced  by  the  hydrostatic 
pressure . 

The  deflected  shape  of  the  panel  due  to  loading  by  a  uniform 

40  41 

pressure  is  the  same  as  the  fundamental  mode  shape.  It  is  given  by  * 


W  / 
=  T7T33  (C 


Ytrx  A  ,  Yttx1 

cos - +  0.133  cosh  - 

,  a  a  , 


42  4 

where  W  =  KP  b  /D 
o  H 


Since 


2  2 

yVw 


dl  1.133  a  { 
i..2  I 


cos  -0.133  cosh 
a2  a  a 


From  Equation  (13)  then  we  can  find  R  =  R  _  the  radius  of  curvature  at 
the  center  of  the  plate  in  the  direction  of  a.  Similarly  we  can  find  the 
radius  of  curvature  of  the  center  of  the  plate  in  the  direction  of  b.  For 
either  case  the  effect  of  curvature,  or  the  mean  radii  Ra  and  R^,  may  be 
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3g  » — 

estimated  by  using  R  =  Rc y2.  The  effect  of  the  curvature  upon  is  then 

found  by  substituting  the  estimated  mean  curvature  in  Equation  (II). 

Similarly,  the  mean  stresses  a  and  a  are  estimated  by  letting 


(a  ) 

^  mo  v* 


(a  ) 

v  mo  v  \ 


'  max'x  „  '  max'y  ,  .  .  ...  . 

a  -  - 1  a  =  - i.  where  (a  ).  is  a  maximum  bending  stress  along 


max  1 


i  derived  from  Reference  42.  The  effect  of  the  stresses  on  w  is  then 

n 

found  by  substituting  the  estimated  values  of  a  and  a  in  Equation  (II). 

x  y 

For  a  30  in.  x  24  in.  x  3/16  in.  clamped  plate  in  a  horizontal 

plane  subject  to  a  uniformly  distributed  hydrostatic  pressure  of  1.5  psi 

the  effects  of  stress  and  curvature  were  found  in  Reference  39.  The  results 

show  that  the  natural  frequency  is  more  affected  by  stress  than  curvature. 

Moreover,  while  stress  and  curvature  due  to  hydrostatic  loading  caused  the 

* 

natural  frequencies  of  the  panel  to  increase,  the  virtual  mass  effect 
(i.e.,  fluid  loading)  caused  a  more  significant  reduction. 


If  hydrostatic  loading  is  considered  in  computations  then  the  plate  (in 
an  infinite  baffle)  cannot  be  considered  to  be  submerged  in  an  infinite 
water  medium.  The  effect  of  the  proximity  of  the  water  surface  on  the 
virtual  mass  must  be  treated.  However,  for  a  plate  located  more  than  A/6 
below  the  waterline  the  presence  of  the  free  surface  will  have  no  signifi¬ 
cant  effect  on  the  fluid  loading.  This  is  clearly  a  frequency-dependent 
43 

criterion.  If  the  plate  is  not  in  a  horizontal  plane  then  the  variation 
of  pressure  with  depth  necessitates  an  integral  formulation  for  the  hydro¬ 
static  pressure.  The  problem  is  then  to  determine  the  center  of  pressure. 
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